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TOLKMITT’S BACKWATER AND DROPDOWN CURVE TABLES 


SYNOPSIS 


Tables and formulas are presented for the computation backwater and 
dropdown curves for channels parabolic section. These are followed 
cautions the selection the coefficients and when used connec- 
tion with the computation water-surface curves. 


INTRODUCTION 


The integration equations for varied flow natural and artificial chan- 
nels continues interest hydraulic engineers. This evidenced 
the many papers published the Transactions and Proceedings the ASCE, 
together with their applications the determination backwater and drop- 
down curves; excellent list references and investigations was given 
Ven Chow 1955.2 Several papers have been published these subjects 
since 1955, among which that Joe Lara and Kenneth 

Most the tables available for the computation water-surface curves 
have been developed for channels rectangular, trapezoidal, and circular 
cross section.2 The formulas and tables Tolkmitt4 were developed for 


open until October 1960. extend the closing date one month, 
written request must filed with the Executive Secretary, This paper part 
the copyrighted Journal the Hydraulics Division, Proceedings the American 
Society Civil Engineers, Vol. 86, No. May, 

Engr., Cons, Engr., Grand Junction, Colo. 

“Integrating the Equation Gradually Varied Flow,” Ven Chow, Proceed- 
ings, ASCE, 838, November, 1955. 

“Two Methods Compute Water Surface Profiles,” Joe Lara and Kenneth 

Schroder, Proceedings, ASCE, Part Vol. 85, No, April, 1959. 


“Wasserbaukunst,” (River Engineering), Tolkmitt, Ernst and Sohn, Berlin, 
1948, 
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channels parabolic cross section and the writer’s knowledge are only 
found the German texts that author. 


FORMULAS AND TABLES 


1932, Boris discussed Tolkmitt’s formulas only briefly. 
Bakhmeteff’s formulas and tables are intended applicable channels 
any type, although many cases they are not easily applied. The tables 
Bakhmeteff have been recomputed and extended Chow,2 and their use illus- 
trated for rectangular, trapezoidal, and circular channels. Since appears 
that Tolkmitt’s tables may fill gap the material available aid the 
computation backwater curves, andalthough admittedly the gap not seri- 
ous one, they are submitted herewith. The tables were furnished some years 
ago van der Veen, then member the Chihli River Commission. 


also supplied the following notes with the formulas for the backwater and drop- 
down curves: 


Backwater curve caused the raising the water level certain 
point river. (For example, weir, dam sluice gates.) 


Tolkmitt (Wasserbaukunst, River Engineering) gives for the calculation 
the distance over which certain raising the water level cer- 
tain point river affects the water level above that point, the follow- 


ing formula: 
d+z 


which distance which raising water level raising water 
level causing the backwater curve; original depth section the struc- 
ture; function; slope fall per (unit) length; coefficient formula 


Backwater curve caused the sudden lowering the water level 
certain point river. (For example, cuttings, the removal rapids, 
the opening weirs, dykes, etc.) 


Tolkmitt (Wasserbaukunst, River Engineering) gives for the caclulation 
the distance over which certain lowering the water level cer- 


tain point river affects the water level; above that point, the follow- 
ing formula; 


which distance which the lowering the water level the 


“Hydraulics Open Channels,” Boris Bakhmeteff, McGraw-Hill Book Co., 
Ed, 
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TOLKMITT'S CURVE TABLES 


lowering the water level causing the backwater curve; original mean 
depth the river; function; slope fall per (unit) length; coeffi- 
cient formula CVR and 


The tables were entirely recomputed and enlarged the addition twenty 
items each table. The formulas were modified@ order facilitate the 
computation, and the new computations were carried out six decimals, then 
rounded four places correspond with the usual tables Bresse’s func- 
tions for rectangular channels. The results are given Tables and and 


TABLE 


0.6196 

0.6390 

0.6574 1.1637 
0.6749 
0.6917 
0.7078 
0.7382 
0.7665 
0.7931 
0.8183 1,2622 
0.8422 
0.8651 
0.8869 1.3022 
0.9080 
0.9479 
0.9670 1,3539 
0.9854 
1,0209 
8.9996 
1,0712 
14,9999 
20.0000 
1.1186 1.7404 


ogico 


2.3756 
2.5809 
2.6829 
2.7847 
3.4922 
3.9948 
5.9985 
6.9990 


PS 


the appropriate backwater and dropdown sketches are presented Figs. land 
Table Eqs. and were used with (d+ z)/d and Table Eqs. 
and were used with z)/d. Table the functions are founded 
decimals, those values that have horizontal bar over the the fourth 
decimal place should have the number the third decimal place increased. 
The use tables the computation backwater and dropdown curves has 
been fully illustrated and explained elsewhere and, therefore, there dis- 
cussion this procedure included herein. However, some suggestions and 
cautions the application formulas this type work are offered. 


The formulas were rearranged for more convenient recomputation the tables 
shown below each table. 


1.000 1.075 
1.001 -0.5069 1.080 
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The Manning formula has come into very general use, replacing the Chezy- 
Kutter formula very large extent. The Chezy-Kutter formula, however, 
was the one general use the time the Bresse and Tolkmitt tables were 


published. use the coefficient with the Manning formula, Chesley 
Posey® has proposed the formula: 


parabolic, would the depth the vertex the theoretical curve. The 
other symbols have the usual connotation. 


TABLE 


1.8038 
1.7701 
1.7135 
1.6893 
1.6672 
1.6468 
1.6103 
1.5496 
1.4685 
1.4312 


deriving these formulas for the water-surface curve, the coefficient 
was assumed constant, but this not always the case. fact, most 
natural channels will vary considerably, except some carefully selected 
reaches. Hence, average values are usually taken. Another assumption was 
that the top width, the section was wide that the hydraulic radius, 


Chesley Posey, Engineering Hydraulics, Hunter Rouse, John Wiley and 
Sons, Inc., New 
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TOLKMITT'S CURVE TABLES 


could assumed equal two-thirds the central depth, But again, 
wide and deep sections, this will not the case. 

channel parabolic form, the top width the water section, 
twice the center depth, then the hydraulic mean depth, will 0.45 
16, will 0.66 For greater values this ratio, the hydraulic mean 
depth may taken two-thirds the central depth the parabola without 
appreciable error. Let T/y equal the formula Then for any par- 
abolic cross section with the ratio top width depth equal greater 
than two, the value the coefficient can computed the formula: 


The average error computed Eq.7 dozen trials was less than 


with maximum error equal practice, when working with 


data obtained from complete hydrographic surveys, Eq. may become mat- 
ter academic interest only, since areas and wetted perimeters channel 
cross sections are usually obtained from plots the sections. 


SELECTING AND 


backwater curves are computed for widely different volumes 
flow given stream, for low water, average and flood flows, and the 


FIG, FIG, 


stream bed subject scour, often difficult, not impossible, de- 
termine the characteristic depth (that is, average central depth 
rectangular parabolic section) during high stages flow. Experience has 
shown that for every foot rise during floods, there may 
scour, even more. such cases, the exercise judgment based 
experience necessary the selection coefficients and factors 
used any formula. The selection the coefficient roughness, case 
point. Experience indicates that the value this important coefficient may 
change much from 0.030 0.025 0.020 with large change the hy- 
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draulic radius from low flow flood conditions channel.7,8 quote only 
one authoritative statement, 


“For example, the Manning coefficient roughness changes great deal 


practically any channel when the variation depth extends con- 
siderable range.” 


avoid the necessity making this adjustment, some engineers have used 
exponential formulas various forms, this means, was intended 
make the coefficient roughness purely descriptive, was originally in- 
tended be. The more common practice present, however, use the 
Manning formula and adjust this coefficient one’s best 

Additional material for those who may wish study some 
the more recent papers and discussions onthe subject backwater curves. 


CONCLUSIONS 


While not new, Tolkmitt’s tables are made available American engineers 
slightly enlarged form, and the use the Manning velocity formula rec- 
ommended with care selecting the necessary coefficients. 
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HOOD INLET FOR CLOSED CONDUIT 


Fred Blaisdell,! ASCE 


SYNOPSIS 


The hood inlet for closed conduit spillways formed cutting circular 
pipe angle and placing the pipe that the crown projects beyond the in- 


vert. This forms hood over the pipe entrance. 

The use the hood inlet causes the pipe tofill with little submergence 
the inlet crown. This filling occurs even though the pipe laid slope 
that hydraulically steep. Smooth pipes having slopes greater thanabout 1.5%, 
and corrugated pipes having slopes greater than about are hydraulically 
steep. Therefore, many highway culverts and other drainage structures fall 
into this steep category and will not fill unless the entrance receives special 
attention. 

The capacity and performance the spillway for variations the hood in- 
let length, the conduit slope, the wall thickness, and the approach conditions 
are described. The great effect vortices the spillway capacity shown 
and anti-vortex devices are developed. Scour the vicinity the hood inlet 


determined for various sizes stone and equations for the scour-hole di- 
mensions are presented. 


Note.—Discussion open until October 1960. extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE. This paper part 
the copyrighted Journal the Hydraulics Division, Proceedings the American So- 
ciety Civil Engineers, Vol. 86, No. May, 1960, 

Presented the August 1957 ASCE Hydr. Div. Conference Cambridge, Mass. 


proj. Supervisor, U.S. Dept. Agric., Research Service, St. Anthony Falls 
Hydr. Minneapolis, Minn. 
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INTRODUCTION 


Renewed interest the hydraulic performance culverts evidenced 
the number studies reported since 1950. One reason for this that the hy- 
draulics culverts andother types closed conduit spillways not simple 
was once thought. Another reasonis because large numbers these struc- 
tures are being built each year and small savings improvements each 
structure result large total savings. 

Among the recent investigations study the hood inlet which was dis- 
briefly tested Malcom Karr and Leslie Clayton, ASCE.2 
This appeared simple, economical, and easily installed inlet that would 
particularly useful the United States Soil Conservation Service (SCS) 
engineers cooperating conservation districts structures are 
installed the landowner. The agricultural conservation, flood prevention, 
and watershed protection programs the United States Department Agri- 
culture also would have use for this type inlet. Accordingly, the request 
thorough study the hood inlet was initiated the Soiland Water 
Conservation Research Divisionof the Agricultural Research Service (ARS) 
the St. Anthony Falls Hydraulic Laboratory, Minneapolis, Minn. 


PREVIOUS WORK 


This review previous work will largely confined closed conduit spill- 
ways used the SCS its programs soil and water conservation, flood 
detention, and watershed protection with the recognition that this restriction 
scope necessarily excludes noteworthy pertinent information. 

The earliest studies closed conduit spillways were those conducted 
drop inlet spillways under the direction L.H. Kessler, ASCE, 1933 and 
1934.3 few field tests were made during the 1930’s plugging the drop in- 
let, waiting for runoff fill the reservoir, and removing the The bar- 
rels some these spillways were slope steeper than that the hy- 
draulic friction grade line, yet the barrel flowed completely full—a phenomenon 
that previously had been thought impossible when the outlet was unsubmerged. 
result these field tests and other field observations full pipe flows, 
theSCS recognizeda lack knowledge the hydraulics closed conduit spill- 
ways and need for additional laboratory investigations. They, cooperation 
with the Minnesota Agricultural Experiment Station and the St. Anthony Falls 
Hydraulic Laboratory the University Minnesota, established research 
staff the St. Anthony Falls Hydraulic Laboratory 1940. 

Experiments the hydraulics closed conduit spillways were begun 
1941 and have been conducted, with interruptions, since that time the SCS 
and the ARS which research activities were transferredon January 


and Clayton, Bulletin No. 35, Oregon State College, Corvallis, Oreg., Engrg. Ex- 
periment June, 1954, 

“Experimental Investigation the Hydraulics Drop Inlets and Spillways for 
Erosion Control Structures,” Kessler, Bulletin No. 80, Univ. Wisconsin, 
Madison, Wis., Engrg. Experiment Sta., 1934. 

“Hydraulics Closed Conduit Spillways, Part VIII, Miscellaneous Laboratory 
Tests,” Fred Blaisdell, St. Anthony Falls Hydr. Lab., Tech. Paper No. 19, Series 
Minneapolis, Minn., March, 1958, 20. 


: 

4 

1954. During the early and middle 1940’s, the theory closed conduit spill- 
“Model Studies Inlet Designs for Pipe Culverts Steep Grades,” Karr 


way hydraulics was developed, relationships were testedand 
test methods were standardized, and methods analysis were de- 
veloped. Refinements were, and continue be, made experience indicates 
their desirability. 

The earlier studies the drop inlet must barrel diameters 
the barrel entrance For many installations this drop inlet 
height much greater than economically desirable use. has been 
found possible reduce the drop inlet height two barrel diameters 
the barrel entrance But even this low drop inlet under 
some circumstances, the search continued for lower drop inlet or, prefer- 
ably, its complete elimination. 

The investigation the hood inlet Karr and indicated that 
would probably fill the need the SCS for simple closed conduit spillway in- 
let which would insure that the barrel would fill even were laid steep 
slope. The studies reported this paper can considered the develop- 
ment the hood inlet originated Karr and Clayton. 

shown Figs. and the hood inlet formed cutting pipe 
angle and laying the pipe that the longer side the crown. This forms 
hood over the pipe entrance. The name was selected ARS and SCS con- 
ferees being descriptive the inlet. 


THEORY 


hydraulic laws beyond those elementary hydraulics are involved 
computing the flow through closed conduit spillways. However, their appli- 
cation the hood inlet will reviewed for the case spillway having 
steep barrel slope. 

The definition conduit slope order. The slope defined steep 
the slope the conduit steeper than the slope the friction grade line. 
The pressure willordinarily sub-atmospheric throughout most the conduit 
length the outlet unsubmerged. contrast, conduit with flat slope has 
slope less than that the friction grade line. The pressures the conduit 
will then greater than atmospheric, except the crown near the outlet when 
the outlet unsubmerged. The friction grade line slope represents the rate 
loss energy due pipe friction. the hydraulic grade lineslope when 
the losses are due pipe friction alone, losses due bends, enlargements, 
and other sources being excluded. 

Weir Flow.—The conduit entrance weir the conduit slope steep 
and the head less than the pipe diameter The head-discharge relationship 
expressed the familiar equation 


“Hydraulics Closed Conduit Spillways, Part Theory and its Application,” 
Fred Blaisdell, St. Anthony Falls Lab., Tech. Paper 12, Series Minn- 
eapolis, Minn., January, 1952, Revised February, 1958, 11. 

“Hydraulics Closed Conduit Spillways, Parts through VII, Results Tests 
Several Forms the Spillway,” Fred Blaisdell, St. Anthony Falls Hydr. Lab., 
Tech, Paper 18, Series Minneapolis, Minn., March, 1958, 24. 
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FIG, 1.—HOOD INLET AND ANTI-VORTEX DEVICE 


Hood Length 


FIG, 2.—HOOD INLET DIMENSIONS 
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which represents the discharge, dimensionless coefficient, de- 
notes the area the barrel corresponding the head over the invert the 
entrance, and the gravitational constant. can put into semi- 
dimensionless terms both sides are divided the product the pipe area 
D2/4 and the square root the pipe diameter After rearranging, this 
results 


p»/2 Ww 4 A D 0220606 
from which can seen that 
5/2 


where has the dimension per sec English units. 

the barrel slope flat, the head-discharge relationship will deter- 
mined exit conditions, such critical depth the tailwater level, instead 
the entrance acting weir. 

Pipe Flow.—Pipe flow governs the head-discharge relationship when the 
conduit full water. The head between the headpool surface and the point 
which the hydraulic grade line pierces the plane the exit freely dis- 
charging pipe consumed entrance, friction, exit, and other losses such 
bends, contractions, and expansions. All these losses are related the 
conduit velocity that 

which the entrance loss coefficient, the exit loss coefficient (as- 
sumed 1.0 for these studies), the Weisbach friction factor, and the 
conduit length. Letting 


where the elevation difference between the entrance and exit inverts, 
cos the vertical distance above the exit invert which the hy- 
draulic grade line pierces the plane the exit, and the sine the barrel 
slope. Solving for 


Correction for Viscous Resistance.—The Froude model law was used de- 
fine the model-prototype relationship since the flow through the spillway 
caused gravitational forces. However, the resistance full-pipe due 
viscous forces and these forces are sufficient importance that the Reynolds 
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model law must considered. Because these two laws are incompatible when 
the same fluid used model and prototype, corrections must made for 
the viscous effects the models are operated accordance with the Froude 
law. This correction accomplished through the relation 
Vp2 
The slope the hydraulic friction grade line for uniform established flow 
rel sufficient permit reliable evaluation the friction slope therefore was 


made part the experimental apparatus. Multiplying both sides Eq. 
the conduit length gives 


the friction loss term Eq. The expression for Eq. evaluated 
Eq. Since takento 1.0, Eq. can evaluated. Thus, although 
the model operated with the Froude correctionfor Rey- 
nolds law effects permits the evaluation the entrance loss coefficient and 
makes the entrance loss coefficient independent the viscous resistance. 

Similarity During Water-Air Flows.—Quantitative simi- 
larity water-air mixtures between model and prototype was not attempted. 
Fortunately, the presence air during slug mixture flows does not cause 
serious difficulty predicting the quantitative performance full-size struc- 
tures. This has already been demonstrated for the case the dropinlet closed 
conduit 

For the hood inlet, steep conduit slope, and increasing discharge, will 
shown subsequently that air flow begins when the inlet seals and slugs form 
near the inlet. These slugs suck air through the inlet they travel down the 
conduit. The amount air determined the water seal the inlet which 
the air has penetrate and the suction from the traveling slugs. 
the water flow increases, the headpool rises slightly and the seal the in- 
let becomes thicker. the same time, the slug length and the resulting suc- 
tion become greater. The net effect increase air flow. The air flow 
balances the slug air demand and the slug air demand governed the sum 
the slug lengths. The air flow the model will different for the proto- 
type than simple scaling would indicate, but the actual head over the inlet will 
little, any, different than that expected from Froude law relationships. 

Maximum air flow willoccur when the slugs are sofrequent that they com- 
pletely fill the conduit. This also the beginning water-air mixture flow. 
The head through the spillway for mixture flow determined density 
the water-air mixture and this effective head, less the friction head due vis- 
cous resistance flow, determines the amount air drawnthrough the water 
seal the inlet. Increasing water flow results decreasing air flow until the 
air flow stops completely and ordinary pipe flow controls the head-discharge 
relationship higher flows. the case slug flow, the air flow rate 
the form water-air mixture will, without doubt, different the proto- 
type than would predicted from model-prototype relationships. But the ac- 
tual air flow immaterial. The important consideration that the head- 
discharge relationship and the differences between the model and its prototype 
will, all likelihood, minor importance the Froude lawis used 
the model test results prototype quantities. 
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HOOD INLET 


The head and the discharge the beginning and end air flow certainly 
can predicted from the models with the usual reliability. For weir flow, 
air has effect the model-prototype relationship. For increasing water 
flow, the end weir control the beginning slug (and and the head 
and discharge this point can definitely reliably projected 
the prototype. Similarly, the beginning full pipe flowis the end air flow 
and the head and air flow can established and 
similarly projected the prototype. Therefore, only during air flow there 
any uncertainty all the model-prototype relationship. The prototype air 
flow can expected exceed the relative air flow the model because the 
greater actual prototype velocities. The relatively greater air the 
prototype will replace and thereby reduce the flow water or, considering 
the headpool level, greater head the inlet will required achieve the 
same water flow. The hood inlet model tests show that large change rate 
water flow causes only small change the headpool level during slug and 
mixture change relative discharge producing only unit rela- 
tive change head. Therefore, quantitative dissimilarity air flow between 
model and its prototype produces relatively much less dissimilarity the 
headpool levels. 

The fact: that exact quantitative similarity according the Froude law 
exists the beginning and the end air flow, that duringair flow the in- 
crease headpool level due dissimilarity air flow quantities over that 
predicted according the Froude relationship relatively small, and that 
previous studies different sized drop inlet closed conduit spillway models® 
have shown that the Froude law can used predict prototype water dis- 
charges during the flow water-air mixtures, sufficient justification for 


the earlier statement that air flow does not cause serious difficulty predict- 


ing the quantitative performance full size structures from tests made ontheir 
models. 


EXPERIMENTAL PROGRAM 


program experiments was planned develop the hood inlet could 
used under wide variety field conditions. The experiments were de- 
signed determine the optimum hood length for conduit slopes 36%, 
the effects conduit wall thickness and various approach conditions, the amount 
scour the vicinity the inlet, and the type andsize anti- 
vortex device required. All tests were carried out with laboratory apparatus 
especially designed for the closed conduit spillway study. 


SPILLWAY PERFORMANCE 


The test results will presented two sections. The first section will 
deal with the overall performance the spillway. The second section will 
devoted the spillway capacity. each section the various items investigated 
will discussed turn. 

Optimum Hood Length.—The optimum length the hood, measured parallel 
the conduit centerline, was determined both quantitatively and observation 
for six different conduit slopes. From five toeight different hood lengths were 
tested each evaluate the relationship between hood length and prim- 
ing head. All tests were made with thin-wall, re-entrant simulate 
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FIG, PRIMING HEAD FUNCTION HOOD LENGTH 
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HOOD INLET 


what was thought the most severe operating conditions. The results 
these tests are summarized Fig. 

Priming ordinarily took place sealing off the conduit near the inlet. This 
process illustrated for good hood inlet Fig. The headpool level 
low Fig. 4(a) the inlet acts weir. Close observation shows that the 
water surface inside the inlet very slightly above the headpool water 
level. The head Fig. 4(b) 0.91 and “ears” are developing just inside 
the inlet. These ears are local rise the water surface that triggers the 
priming action. The head has increased 0.98D Fig. the ears have 
grown until they almost meet. Twenty seconds later the head had increased 
0.99 the two ears had joined together, and the inlet had primed shown 
Fig. 4(d). another sec the head had increased 1.01 the inlet was 
completely full water, and the pipe was full fora lengthof about Das shown 
Fig. 4(e). the pipe creates additional head due the slope the 
pipe. This increases its capacity. Since there insufficient seal over the in- 
let, air sucked shown Fig. 4(f). Further increases flow and head 
result the formation slugs that travel through the pipe and suck air 
shown Fig. 4(g). These slugs are hydraulic jumps that form near the inlet 
and travel through the conduit. The slugs increase frequency until the con- 
duit full water-air mixture. the water flow increased, the air flow 
increases then decreases. Eventually the air flow stops, the 
conduit flows completely full water shown Fig. 4(h), and the laws 
pipe flow control the head-discharge relationship. 

The development the ears which determine the priming depends the 
hood length. This illustrated Fig. where the head for each photograph 
approximately equal the pipe diameter. The contracted jet expands and 
hits the inside wall below the horizontal centerline for the D/4 and the 
D/8 hood lengths shown Figs. 5(a) and 5(b). The point impingement 
higher for the D/2 hood length shown Fig. 5(c) but the air-filled entrance 
contraction readily apparent. There appears air the entrance 
contraction for the D/8 hood length shown Fig. 5(d). Also, ears are ap- 
parent. The ears have reached across the entrance crown and have sealed the 
inlet for the D/4 and hood lengths shown Figs. 5(e) and 5(f). 

The average priming head decreases the hood length increases until 
minimum reachedat hood length D/2 D/4. This shown Fig. 
There seems advantage using hood length longer than D/4, and 
this length recommended the minimum satisfactory length. This mini- 
mum hood length could also have been closely determined visual 
vations alone, illustrated Fig. 

Conduit major effect the slope the conduit the per- 
formance inthe additional head which becomes available the slope and the 
total dropis increased. The effect slope large for the range flows when 
air passes through the for full pipe flow but small for weir con- 
trol. 

The slope affects the spillway performance for low discharges the weir 
control range only for slopes flat that the weir the entrance not the 
control section. This shown Fig. where the only weir flow control data 
that deviate from the composite data that for the zero conduit slope. 

Full-flow discharge through the spillway increases with the conduit slope 
shown Fig. The conduit had the same length for all slopes, the 
total drop through the spillway increased with the slope. The increase total 
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(a) H/D 0.77. Weir flow. ears 0.91. Weir flow. Ears 
on water surface inside inlet. are just beginning to appear. 


(c) 0.98. Just before inlet seals. 0.99. Just inlet seals. 
Ears from opposite sides almost meet. Ears from opposite sides have joined. 


Just after inlet has (f) H/D Air bubbling in. 
sealed. There is no air entering. Pipe perimeter is wet about 2D. 


Hydraulic jump traveling (h) Pipe flow. 
and sucking air enters inlet. 


FIG, HOOD INLET HEAD- 
POOL LEVEL RISES AND INLET PRIMES 
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head due increasing the slope actually what causes the increased flow. 
There was change the entrance loss coefficient with conduit slope. 

Vortex Inhibitor.—A device control the formation vortices the en- 
trance closed conduit spillways necessity. The greatest need for the 
anti-vortex device low submergences the inlet where vortex may let 
air into the spillway; vortex inhibitor not needed for weir control the 
entrance and not essential when the inlet deeply submerged. 


(e) Hood length (f) Hood length 


FIG, HOOD LENGTH FLOW INLET FOR H/D 1.0 


The air entering through the vortex replaces the water and thereby reduces 
the discharge capacity the spillway. The magnitude the potential reduction 
capacity result vortices apparent Fig. where data are plotted 
for two series tests. The open circles and the curves represent data ob- 
tained withan anti-vortex wall. The anti-vortex wall was removed for the data 
identified dots and run numbers. The data obtained for each run are con- 
nected showthe order which they were obtained. Continuous curves could 
have been obtained enough points had been read from the recorder chart. 
However, the data plotted are sufficient illustrate the wide variation the 
discharge the vortex strength varies. 
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Series Symbol Slope 

| Full pipe flow 
| 
| 
flow and mixture flow Air sucked and carried 
through conduit. Eq. 
a/p¥2 
FIG, 6.—EFFECT CONDUIT SLOPE HEAD-DISCHARGE 
RELATIONSHIP. HOOD LENGTH 3D/4 
Series Symbol Anti-vortex Device 
run number None 
ANTI-VORTEX DEVICE CAPACITY 
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The discharge unpredictable either magnitude any given time 
when vortices are present. contrast, the completely reliable and predicta- 
ble head-discharge curve shown Fig. can obtainedif anti-vortex de- 
vice used. Since the capacity the spillway Fig. was reduced only 
one-third its potential capacity the elimination the anti-vortex device, 
the importance using some typeof anti-vortex device amply demonstrated. 

The anti-vortex device used during the series plotted Fig. splitter 
wall mounted the outside crown the inlet shown Fig. 8(a). This anti- 
vortex wall, the alternate shown Fig. 8(b), was used throughout most 
the test program with uniformly satisfactory results. Although vortices some- 
times formed, they admitted air the conduit only briefly and did not cause 
measurable reduction the spillway capacity. 

Objections the splitter type anti-vortex wall were voiced before the tests 
were completed and flat plate was suggested Paul Jacobson, ASCE, 
Soil Conservation Service State Engineer for Iowa. Observational tests 
that the minimum size anti-vortex plate that dimensioned Fig. 
8(f). 

Conduit Wall Thickness.—A wide range pipe wall thicknesses was tested 
with the D/4 hood inlet length and 20% conduit slope. The thickness the 
thickest wall tested was 0.224 This somewhat thicker than concrete pipe 
having the greatest relative wall specified bythe American Society 
for Testing Materials. The thinnest pipe wall tested was 0.015 At- 
tempts were made totest wall thicknesses 0.011 and 0.006 but the Lucite 
model was thin (0.024 in. and 0.013 in., respectively) that the inlet broke 
before data could obtained. Even these thin walls thicker than 
the walls metal pipe. Seventeen different wall thicknesses were tested. 

The wall thickness had apparent effect the priming action general 
inlet performance. Wall thickness did affect the loss the entrance. The 
magnitude this effect described subsequently. 

Approach Conditions.—A limited number different approach conditions 
were tested. Most tests were conducted with the inlet re-entrant, that pro- 
jecting into the headpool, since this was assumed the most severe test 
inlet performance. The inlet performance was checked with the inlet invert 
intersecting the sloping face dam and with berm the dam the inlet 
invert elevation. The face had slopeof 3.37 and the berm had 
width 4.22 Two conditions each dam shape were tested: For one 
condition erosion was prevented plastering the sand dam with neat cement 
the other condition scour hole was allowed develop under maxi- 
mum flow and the test was run with the scour hole. 

The various approach conditions had effect the general performance 
the spillway although the presence the dam did reduce the entrance loss 
somewhat. The range approach conditions tested limited and deserves 
further study. Nevertheless, unless the approach greatly restricted, ap- 
approach conditions will not have major effect general 
ing characteristics the type flow the spillway. 

Scour Hood Inlet.—The velocity the sharp-edged inlet theoretically 
infinite. only logical expect scour the vicinity the inlet the ma- 
terial there erodible. Fortunately, the velocity decreases rapidly with dis- 
tance away from the inlet and this limits the size the scour hole toler- 
able amount. series tests was performed evaluate the dimensions the 
hole that would scoured various sizes noncohesive materials. second 
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(c) Splitter Type Anti-vortex 
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(f) Anti-vortex Plate 
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objective these tests was determine the size riprap required pre- 
vent scour. Anda third objective was the determination the height the 
hood inlet above the erodible material which scour would occur. These 
objectives were achieved conducting tests four different discharges 
eight different sand sizes. 

The surface radius the scour hole plotted Fig. the case 
the berm, the radius taken one-third the sum the scour hole width and 
upstream length. When there berm, the radius taken one-half the 
width the scour hole. Variation the results expected. The re- 
markable thing about Fig. that the variationis small enough reasonable 
curves can drawn through the data. The curve for quite well 
defined. The same slope was used for the other values This slope 
fits the data well any other slope. The slope the curve 1/5. in- 
spection, the equation the curves Fig. 

1/5 
where the actual bed material size. 

The data can used determine the grain size, dj, the bed material 
which will just picked given discharge. These data are plotted 
Fig. for those flows when only few grains were picked for which 
movement was classified imminent. surprisingly well-defined curve 


obtained which, for the data available, applies grains greater 
diameter. This curve has the equation 


0.075 


The depths the scour holes are plotted Fig. 11. Although there 
considerable scatter the data, this not unexpected. The curves drawn 
Fig. were started the zero depth scour given The slope 
the curves was determined the better defined data obtained for the two 
smallest grain sizes. The same slope was used for all curves. The curves 
Fig. define the depth the scour hole and have the equation 


Eq. gives fair results for bed material sizes greater. The 
precision this equation low but may serve aid judgment. 

Eqs. 10, and represent the data. safety factor should 
used for design purposes. For example, stone 0.89 size was picked 
when slightly disturbed. passed through the pipe and smashed the inlet. The 
stone was heavy, having specific gravity 2.9. This the necessity 
using riprap adequate size, paving prevent scour, allowing the 
scour hole develop. The minimum riprap size can determined from 
Eq. 10, the minimum area paved can determined from Eq. and the 
approximate depth the scour hole can determined from Eq. 11. Appropri- 
ate safety factors should applied the computed values. 

Eq. can used determine the height the inlet invert above the bed 
material scour occur. This height computed for one bed material 
and flow. Tests made both with and without the berm confirmed the prediction 


that scour would take place the inlet located proper distance above 
the erodible bed. 
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0.00! 0002 0.003 0.005 002 003 005 Ol O3 05 1.0 
Bed Size Pipe Diameters 
FIG, 9.—SCOUR HOLE RADIUS 
Symbols 
With berm 
x Without berm ve 


Bed Size 
Pipe Diameters 


| | | 20 

ooo fo) fo} fo! 


FIG, 10.—GRAIN SIZE IMMINENT 


MOVEMENT 
FIG, 11.—SCOUR HOLE DEPTH 


‘al 

50 
Q/ 13) . 
2) 

02 

: 

4 


HOOD INLET 


The scour equations can used determine the minimum size riprap 
required the inlet crest located various distances above the surface 
the riprap. the distance below the inlet crest and the discharge are sub- 
stituted Eq. 11, the size riprap that will just picked can deter- 
mined. Also, the depth below the crest which the available riprap will just 

Finally, smaller sizes riprap can used with increasing distance away 
from the inlet. These distances for various riprap sizes can computed from 
Eq. 

These scour tests have defined reasonably well the relationships between 
the bed material size, the discharge, and the scour hole dimensions for vari- 
ous sizes conduit. Minimum dimensions are given Eqs. 10, and 11, 


and essential that adequate safety factor applied these equations 
are used for design purposes. 


SPILLWAY CAPACITY 


low discharges the control the inlet acting weir. intermediate 
discharges the control changes from weir pipe and air carried through 
the spillway with the water the hood length sufficient cause the inlet 
prime. high discharges the control for full pipe flow. The spillway 
rating curve may determined once the weir rating curve and the entrance 
loss coefficient are known. 

Capacity capacity the spillway with the entrance acting 
weir may affected the hood length, the conduit slope, the wall thick- 
ness, the approach conditions. The effect each these items will 
discussed turn. 

Hood Length.—The effect hood length the weir flow head-discharge 
relationship was determined for all conduit slopes, the data for all hood lengths 
being plotted single sheet for each conduit slope. 

There some scatter the data but trend that indicates any effect 
hood length the head-discharge relationship. This somewhat surprising 
view the wide range hood lengths, but other conclusion possible. 
These comments apply the data for all conduit slopes tested. 

Conduit Slope.—The conduit slope was found have small but significant 
effect the weir flow head-discharge relationship. The data for one hood 
lengthand six conduit slopes The other hood lengths 
show similar results. 

The control for the 0.025 and greater conduit the hood inlet crest 
and the family curves. However, the control for the zero con- 
duit slope was critical depth the conduit exit. The head-discharge relation- 
ship the inlet therefore dependent the conduit friction loss addition 
the entrance losses. will different for each conduit length and rough- 
ness. Although the data for the zero slope are shown Fig. 12, they should 
not used for design purposes. The data can used for design purposes 


when the conduit slope sufficiently steep the control always the 
inlet. 


Conduit Wall Thickness.—Head-discharge data covering the complete range 
wall thicknesses tested were plotted and compared. The data scatter some- 
what but the scatter not greater than one can reasonably expect. There 
systematic trend the scatter and single curve can drawn that ade- 


quately represents all wall thicknesses spite the wide range thicknesses 
tested. 
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Conduit Slope 
FIG, CONDUIT SLOPE WEIR FLOW HEAD- 
DISCHARGE RELATIONSHIP. HOOD LENGTH 3D/4 
FIG, 13.—COMPARISON HOOD INLET WEIR FLOW HEAD- 
DISCHARGE EQUATIONS WITH MAVIS CURVE 
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The conclusion that the wall thickness the square-edged re-entrant 
hood inlet has apparent effect the weir flow head-discharge relationship. 

Approach number different approach conditions tested 
too limited permit any general statements regarding their effect wier 
flow. The approach conditions tested include the sloping dam face and berm 
the dam face the crest elevation with both fixed and scoured surfaces. 

The data show effect preventing scour allowing the scour hole 
develop. However, the presence the dam appears increase the dis- 
charge approximately 10% over the condition for re-entrant inlets given 
Eq. but the present data are too few justify general recommendations 
increased flow for certain approach conditions. 

Discharge Coefficient.—The weir flow coefficient Eq.3 was computed 
for the data shown Fig. 12. constant, was found that linear 


5/2 
relationship existed between and different linear relationship 


was found for each conduit slope The relationship between and H/D 
for the square-edged hood inlet was found 


which, when substituted Eq. gives 


1/15 


The curve representing Eq. shown Fig. 12. The agreement with the 
data seen excellent about 0.9 H/D. However, the equation starts 
deviate from the data about H/D 0.8 and the curves have been dotted 
above this value. The dashed extensions the head-discharge curve above 
H/D 0.8 were drawn based the shape the curve given for the equation 
below H/D 0.8 and can seen that they represent the data within the ob- 
tained precision. Thus, Eq. may used obtain the weir flow head- 
discharge curve for square-edgedinlets and conduits steep heads 
between and H/D 0.8. The curve may extended obtain the head- 
discharge relationship for the higher heads. 

shown later that weir flowis not the control for heads those 
given Eq. 13. The curve Eq. shown Fig. define this upper 
limit weir control. 

Comparison with Mavis’ Curve.—The work Mavis® culvert hy- 
draulics has been widely used and often quoted. comparison Eq. with 
the weir flow head-discharge curve presented Mavis for square-edged en- 
trances therefore interest. The principal difference between the inlets 
that Mavis used plane headwall while here the inlets were re-entrant. 
Also, Mavis’ range slopes from 0.007 0.0657 and has indicat- 
effect slope. The comparison made Fig. 13. interesting 
note that the agreement generally good. 

Slug and Mixture Flow Capacity.—The gently sloping line Fig. gives 
the head-discharge relationship for slug flow and water -air mixture flow through 
the conduit. also limit curve giving heads and discharges below which 
the control will weir flow and above which the control will pipe flow. 


Hydraulics Culverts,” Mavis, Bulletin No. 56, The Pennsylvania 
State College, Engrg. Experiment Sta., State College, Pa., 1943, Fig. 23. 
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The equation the slug and mixture flow curve for the D/4 hood length 
may approximated 


Extrance Loss for Full Conduit Flow.—The loss coefficient for full 
flow affectedconsiderably the hood length, the conduit wall thickness, the 
shape the entering edge, and the approach conditions. The conduit slope has 
effect, and the effect the various anti-vortex walls not large. 

Conduit Slope.—The average entrance loss coefficient for each series 
shows trend when plotted against the conduit slope The inlets were re- 
entrant and the water had free access all sides the inlet for all conduit 
slopes. Therefore, the finding that conduit slope has effect the entrance 
loss entirely reasonable. 


Hood Length 


FIG, HOOD LENGTH ENTRANCE LOSS CO- 
EFFICIENT. SQUARE-EDGED INLET, 0.015 


Hood Length.—The effect the hood lengthon the entrance loss coefficient 
definite and considerable. The results the tests are shown Fig. 14. 
This significant effect for full flow contrasts with the lack any such effect 
for weir flow. Fig. demonstrates that there small but definite hydraulic 
advantage using the shortest hood length that will insure satisfactory hy- 
draulic performance. This hood length D/4. 

Conduit Wall the hoodinlet may classed thin 
thick. The entrance loss coefficient constant and the wall classi- 
fied thick when the ratio the pipe wall thickness the pipe diameter 
greater than 0.04. The loss coefficient variable and the wall classi- 


fied thin when less than 0.04. This shown Fig. 15. The data 
plotted include dataon walls thinas was possible test. Thin- 
ner walled inlets were built but they broke before data could obtained. The 
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three low were obtainedearly the test program. The reason why 
they are low cannot explained. The same inlets were tested again with the 
results shown points plotted directly above the low points. 

The theoretical and experimental curves obtained Harris? are 
shown Fig. because they have the same shape the curve obtained for 
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Entrance 


Harris [5] 
Theor etical 


Experimental 


Conduit wall 


FIG, CONDUIT WALL THICKNESS ENTRANCE LOSS 
SQUARE-EDGED INLET. HOOD LENGTH 
3D/4. CONDUIT SLOPE 0.20. 


the hood inlet. The wall thickness obtained Harris which the coefficient 


ceases vary with wall thickness seen the same that obtained for 
the hood inlet. 


Influence Pipe Thickness Re-entrant Intake Losses,” Harris, 


Bulletin No. 48, Univ. Washington, Engrg. Experiment Sta., Seattle, Wash., November 
1928. 
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The test results show, for the square-edged hood inlet having hood length 
D/4, that 


(thin-walled inlets), and that 


(thick-walled inlets). 
Anti-Vortex Device.—The various types anti-vortex devices that were 


tested are shown Fig. Their effect the entrance loss coefficient 
shown Table 


TABLE ANTI-VORTEX DEVICE ENTRANCE LOSS COEFFICIENT 


Series Anti-Vortex Device 
Fig. 0.96 
Fig. 0.95 
Fig. 0.98 
Fig. 1.01 
Fig. 0.85 

106 Fig. 0.83 


TABLE 2,—EFFECT APPROACH ENTRANCE LOSS COEFFICIENT 


Series Approach Condition 
dam face slope, fixed 0.80 
3,37 dam face slope, scour hole 0.93 
3.37 dam face slope, berm wide, fixed 0.84 
3.37 dam face slope, berm wide, scour hole 0.93 
re-entrant 0.95 


Hood Length 3D/4 Conduit Slope 0.20 t/D 0.0156 
Anti-Vortex Device Fig. 


There difference the behavior the two splitter type anti-vortex 
walls shown Figs. 8(a) and 8(b). They are alike that difference should 
expected. Extending the anti-vortex wall down front the inlet and bring- 
ing closer the inlet, Figs. 8(c) and 8(d), causes small increase 
the entrance loss coefficient. The headwall type anti-vortex wall shown 
Fig. 8(e) apparently suppresses the contraction the inlet and causes con- 
siderable decrease the entrance loss coefficient. However, the efficiency 
this wall vortex control somewhat less than the splitter walls. The 
coefficient for the anti-vortex plate shown Fig. 8(f) the low- 
est obtained. This, then, the best anti-vortex device, considering both its 
low loss coefficient and its vortex inhibiting ability. 
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Approach Conditions.—The use pavedapproach the hood inlet results 
significant reduction the entrance loss. This shown Table The 
entrance loss coefficient lowest berm used. 

The entrance loss almost equal that for the re-entrant inlet the ap- 
proach the hood inlet not paved and scour hole allowed develop. 

Such data are presently available different approach conditions show 
that significant reduction the entrance loss possible the approach 
paved prevent scour the dam the vicinity the inlet. 


CONCLUSIONS AND RECOMMENDATIONS 


Only the desirable weir and pipe controls exist for the hood inlet the hood 
length The undesirable orifice control may exist for short hood 
lengths. The shortest hood which will insure full conduit flow minimum 
headwater elevation D/4. 

There influence conduit slope the spillway performance the 
slope steep. 

vortex inhibitor some type necessity. The anti-vortex devices 
shown Fig. are recommended. However, the device shown Fig. 8(e) 
the least satisfactory. 

The capacity the spillway acting Eq. 12;the capacity 
for flow mixture water given Eq. 13;and the entrance loss 
coefficient for full pipe flow and square-edged entrances given 14. 

The composite head-discharge curve for particular spillway having good 
hood inlet can drawn plotting the weir flow curve, the slug and mixture 
flow curve, and the pipe flow curve. The extensions these curves beyond 
their intersections with each other are imaginary that the complete rating 
curve formed tracing, the headand discharge increase, along the weir 
flow curve its intersection with the slug and mixture flow curve, along this 
latter curve its intersection with the pipe curve, and along the pipe curve for 
all higher heads and discharges. 

The thickness the conduit wall did not affect the hood inlet performance 
its capacity aweir. The effect wall thickness the entrance loss 
coefficient given Eqs. 14. 

Although there was little variation approach conditions, the approach had 
influence the spillway performance. The presence the dam face does 
somewhat reduce the entrance loss coefficient. 
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APPENDIX NOTATION 


The following symbols, adopted for use the paper, conform essentially 


with “American Standard Letter Symbols for Hydraulics” (ASA 
prepared committee the American Standards Association with Society 
representation, and approved the Association 1942. 


orifice area taken the plane the entrance 

orifice discharge coefficient; 

weir discharge; 


diameter the pipe; 

Weisbach friction factor; 

acceleration due gravity; 

friction head loss; 

local pressure deviation from the friction grade line; 

head over the invert the entrance; 

total head water surface the point which the hydraulic 


entrance loss coefficient; 
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exit loss coefficient; 

length the conduit; 

discharge through the conduit; 
rate change headpool storage; 
sine the slope the barrel; 
pipe wall thickness; 

velocity the pipe Q/A; 


difference elevation between the conduit invert the 
the exit; 


proportion the pipe diameter above the invert the pipe exit 
which the hydraulic grade line pierces the plane the exit; and 


angle between the face the hood and perpendicular the conduit 
axis. 
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UNIFORM WATER CONVEYANCE CHANNELS ALLUVIAL 


SYNOPSIS 


Methods designing uniform alluvial channels are developed and 
ed. Design theories which are recommended include: (1) modification the 
regime theory India, and (2) modification the tractive-force theory. 
Special emphasis given the modified regime theory. The results this 
investigation are based field study stable alluvial irrigation channels, 
and other existing alluvial channel data which are applicable. 


INTRODUCTION 


The design silt stable regime channels has been the object consid- 
erable research during the past four decades. known that any design 
theory must recognize the effect the many variables involved. Originally 
only single equation the Chezy Manning type was utilized. recent 
years, however, various investigators have concluded that width, depth, and 
slope are all variable alluvial channels and that this implies mathematically 
the necessity least three design equations. The problem appears even 
more complex when one recognizes that sediment transport also affects chan- 
nel stability. 

Recently, two schools thought, one empirical, the other primarily theo- 
retical, have each evolved theories more capable adequately predicting chan- 


Note.—Discussion open until October 1960. extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE. This paper part 
the copyrighted Journal the Hydraulics Division, Proceedings the American So- 
ciety Civil Engineers, Vol. 86, No. May, 1960. 

Presented the 1958 ASCE Hydraulics Division Conference Atlanta, Ga. 

Proj. Chf., Branch, Geol. Survey, Fort Collins, Colo. 

Prof. Civ. Engrg. and Dir. CSURF, Colorado State Univ., Fort Collins, Colo. 
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nel behavior than any heretofore. yet, however, these theories are incom- 
plete and should subjected further study and investigation order 
broaden their scope and possibly reduce them into one general comprehensive 
and complete theory. 

Definition Stable Channel.—An excellent definition stable regime 
channels was presented Lane (14),3 1953, follows: 


stable channel unlined earth channel for carrying water, the 
banks and bed are not scoured the moving water, and which 
objectionable deposits sediment not occur. 


The foregoing definition does not exclude minor channel erosion accre- 
tion during the yearly cycle flow. does, however, require that these op- 
posing effects should cancel annual basis. 

Application Stable Channel Theory.—Stable channel theory widely ap- 
plicableto many phases work inthe fields civil and agricultural engineer- 
ing follows: 

Satisfactory design irrigation systems requires intimate knowledge 
channel design relationships. fact, considering the current quality these 
relations, considerable experience necessary bridge the existing gaps. 

Improper design may well introduce instability, artificial channels, 
such magnitude that not economically feasible operate them. 

Although canals constructed for the purpose conveying drainage and 
wastewater from given areas might considered under the foregoing 
ings, they have been subjected such gross neglect that they warrant special 
mention. Certainly lack suitable channel design theory, least lack 
application problems falling this group, has caused loss consider- 
able tillable acreage, many hydraulic structures, and the formation unsight- 
scars the earth’s surface. 

adequate knowledge the many variables influencing channel stability 
and their interrelationships would make possible more intelligent treatment 
of: 


Bank-scour problems caused increasing normal flow channels. 

Evaluation the effect slope changes, due the use cut-offs, 
channel stability. 

Erosion problems upstream from and downstream from bridges and 
other types hydraulic structures that constrict the channels. 

Stabilization non-regime channels armor coating banks with mate- 
rial more resistant scour than the existing natural material and/or in- 
troducing structures control slope. 


particular importance able design stable power canals 
alluvial materials, since the quantity sediment being transported the 
power plant penstocks must controlled—in fact, for the most part, 
nated, 

Research.— Based the inadequacy current design methods and lack 
understanding the regime theory the United States, field study stable 
canals was proposed 1953, with the purpose mind attemptingto clarify, 
expand, and perhaps combine the foregoing theories. Specifically, the major 
objectives this research are follows: 


Numerals parentheses, thus (1), refer corresponding items the Bibliography. 
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investigate the validity and applicability the regime theories, de- 
veloped India, canals the United States. 


investigate, expand, and possibly improve the tractive-force method 
stable channel design. 


relate the regime theories the tractive-force theory, insofar 
possible. 


this paper, the results the study the applicability the regime-type 
relationships channel design are emphasized. 


REVIEW CURRENT DESIGN METHODS 


Satisfactory design and construction artificial channels alluvial ma- 
terial affected many factors, some which are extremely complex and 
hence are only vaguely comprehended. brief history the development 
currently used empirical rules and theories presented herewith. 

Regime Theory.—The regime theory India was initiated Kennedy, 
1895, when produced his classic empirical equation 


whereC and werethought tobe constants and were originally assigned values 
0.84 and 0.64, respectively. Kennedy concluded that channels having veloci- 
ties based his formula would neither silt nor scour their beds. The equa- 
tion was empirically formulated based collected the Bari Doab ca- 
nal system, the Punjab. 

shallow one carry the same discharge. Actually, this far from the truth 
the matter. far design concerned the only time the Kennedy equa- 
tion can possibly yield correct results when shape properly selected, 
when dealing with very stable materials. 

Other investigators attempting prove disprove Eq. soon established 
that the constant and the exponent varied from system system. 
spite this, the Kennedy equation, its original form, can still found 
many recent engineering texts and was applied extensively design, or- 
iginally presented, until about 1930. 

1919, Lindley introduced the following regime equations: 


These equations were derived correlating data obtained from 786 observa- 
tions branch surveys the Lower Chenab Canal. This was the first time 
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that bed width and depth were introduced regime variables. his reply 
the discussion his paper stated: 


The existence these relations meant that the dimensions width, depth, 
and gradient channel carry given supply loaded with given silt 
discharge were all fixed nature, that is, uniquely determined. 


The variables bed width, depth, and slope were observed. Velocities, however, 


were not observed. They were computed means the Kutter and Chezy 
equations, the Kutter equation being: 


The magnitude Eq. was assumed constant 0.0225. The Lindley 
equations, although never popular the Urited States, were used extensively 
India, until about 1935. 

1927, Gerald Lacey was commissioned the Governments England 
and Indiato systematize all had been collected relativeto stable chan- 
nels. summary the results the study were published 1929 and 1933. 
Lacey’s original equations, presented the first paper, were: 


and 


1936, Bose (12), and the staff the Punjab Irrigation Research In- 
stitute, presented the following formulas: 


and 


Eqs. 10, 11, and represent the results several years painstaking col- 


lection and statistical analysis the data. Note the similarity Eqs. and 
10. 
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According Lacey, depends size sediment particles but not af- 
fected concentration the silt load being transported the 


applying the Lacey theory, known, estimated, and bed width, 


depth, and slope are calculated. 

1951 Thomas Blench presented his concept the regime theory. His 
equations are expressed terms and where the mean width and 
the depth average line through the channel bottom such that thearea 
the cross section WD. The design equations presented Blench are: 


least three American engineers, Pettis (15) and Leopold and 
Maddock, Jr. (16), have investigated the applicability the regime theory 
American rivers. The equations presented were the same form the 
regime equations but did not have the same magnitude constants and expon- 
ents. Hence, one might conclude that regime equations depend the condi- 
tions which they are based and are valid only within the range the ob- 
served data. 

Inadequacy Regime Methods.—The Lacey equations, and modifications 
them, have been very popular India, but have never been used extensively 
elsewhere the world. This group equations undoubtedly provides sound 
basis for design currently exists they are used under circumstances 
similar those from which they were obtained. The major disadvantages 
the method are: 


has not based onthe wide variety conditions encoun- 
tered practice; 


and 

involves factors that require knowledge the conditions upon which 
the formulas are based they are applied successfully. 


The regime equations presented Blench modify the Lacey equations 
such way that the effect the side and bed the channel can evaluated 
separately means the side factor and bed factor. This approach seems 
basically more sound than averaging the two effects the Lacey equations 
since effect bed and side conditions flow are vastly different. 

Tractive-Force Method.—Thetractive-force design theory formulated 
the basis that stability bank and bed material function the ability 
the bank and bed resist erosion resulting from the drag force exerted 
them the moving water. 
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This concept has been widely applied the theory sediment transport 
both the United States and other countries but only limited extent 
connection with design channels alluvial material. Use this method 
for design has been suggested Williams (14) and Armin Schoklitsch (20). 
The latter suggests that the relations given Table exist between type 


soils and permissible tractive force, and that these data are suitable for de- 
sign purposes. 


TABLE 1.—VARIATION PERMISSIBLE TRACTIVE FORCE WITH TYPE SOIL 


Permissible tractive force, 
pounds per foot 


Loam 0.062 
Sand 0.102 
Stony and Loamy soil 0.082 
Course gravel 0.205 
Very compact soil 0.256 


The foregoing tractive-force concept was developed further under the lead- 
ership while withthe United States Bureau 
The USBR procedure based the hypothesis that practical canal design 
tractive-force problem—beyond that, the approach new and theoretical. 

Three distinct classes instability have been defined Lane follows: 


Channels subjected scour that not silt; 
Channels which objectional deposition occurs but not scour; and 
Channels which objectional scour and silting both occur. 


Class instability the simplest the three proposed, and fortunately, 
also the one primary importance, since most the present and future 
canal problems are and will clear-water problems. 

The recommended design procedure was developed considering: 


Distribution tractive force over the channel periphery for different 
side slopes with special emphasis the magnitude shear exerted the 
sides compared the bed; 

Relative stability soil particles the bed and the sloping sides 
the channel; and 

Magnitude safe tractive force for different mean sizes and gradations 
non-cohesive materials. 


The shear distribution was worked out mathematically for rectangular chan- 
nels, and membrane analogy and the method finite differences for trape- 
zoidal sections. that maximum shear onthe bed was approximate- 
equal and the sides was about 0.76 this value. 

Various theories have been developed making possible estimate magni- 
tude tractive force resulting from the movement fluid with respect toa 
fixed boundary with which contact follows: 
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which the tractive force exerted the bed point where the water 


depth The derivation Eq. based equilibrium forces tangent 
the channel boundary; 


which the average tractive force exerted the channel periphery (Eq. 
derived using the same procedure that for Eq. 16); 


5.75 log 


which defined the tractive force the channel bed directly beneath 
the vertical which and are measured (The derivation Eq. 
involves use shear theory and the Karman logarithmic velocity distribu- 
section (isovels) establish lines zero momentum exchange. Then, com- 
puting the tractive force the bed between two zero momentum lines (lines 
drawn right angles the isovels) considering the weight water above 
the section bed question that confined between the two zero momentum 
lines. 

Inadequacy the Tractive-Force Method.—The tractive-force concept 
basically sound insofar has been developed. That is, for the design 
channels for the conveyance essentially clear water coarse non-cohesive 
materials. However, its application qualitative for the: 


Design channels fine non-cohesive soils (sand range and finer); 
Design channels cohesive materials; and 


Design channels that are required transport appreciable (in excess 
500 ppm) sediment load. 


FACTORS INFLUENCING STABILITY OPEN CHANNELS 


detailed study existing design methods immediately verifies the com- 
plexities stable channel theory. cover adequately each existing design 
case, any Suitable theory must include the effect all the pertinent vari- 
ables. date, theory has been conceived capable adequately consider- 
ing all them and their influence channel behavior. list the major 
variables involved follows: 


Discharge condition; 
Steady discharge; 
Variation discharge with crop needs; 
Periodic operation; 
Design discharge exceeded; 
Slope; 
Shape channel; 
Boundary material: 
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Cohesive materials; 
Non-cohesive materials, sand range; 
Non-cohesive materials, larger than sand; 


Water temperature; 


Magnitude and range size total sediment load; 
Berms; 


Deposition wash load; 
Shape channels possessing berms; 
The influence berm formation bank vegetation; 


Bank stabilization; 


Secondary circulation transverse flow; 
10. 
11. 


Seepage inflow outflow; 

Effect wind action stability channels: 

Estimate tractive shear exerted the bed channel due 
waves; 


Effect percolation caused waves the stability apermeable 
bed; 


Effect waves bank stability; and 
Wind effect during the non-operating season. 


FIELD DATA 


The data upon which this paper based were obtained from canal studies 
India Punjab (18) and Sind (11)] and the United States Imperial Valley (7), 
San Luis Valley (14), and canals Wyoming, Colorado, and Nebraska (2) and 


The study irrigation canals Wyoming, Colorado, and Nebraska was 


jointly sponsored the Corps Engineers, the United States Geological Sur- 
vey (USGS), the USBR, Colorado State University, and University Wyoming. 
Field data from selected canal reaches these states were taken 
Simons and Bender during the summers 1953 and 1954 (22), hereafter 
referred Simons and Bender data. 


The data taken the field from the selected straight stable canal reaches 
included: 


(1). 
(2). 
(3). 
(4). 
(5). 
(6). 
(7). 
(8). 
(9). 
(10). 
(11). 


Magnitude discharge; 

Velocity distribution; 

Slope water surface; 

Shape canal cross section; 

Suspended sediment distribution; 

Total sediment load, whenever possible; 
Samples bed material and side material; 
Armor coat samples; 

General condition the bed; 
Temperature the water; and 
Photographs. 


total reaches were investigated. Complete information each 
reach presented Tables and 

USBR Data.—Because the limited amount data tractive force for 
coarse non-cohesive material, the USBR (14) investigated canals constructed 
this type material the San Luis Valley Colorado. The canals were 
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located alluvial fan. The size the natural material varied consider- 
ably, decreasing size from the apex outward. The canals constructed the 
cone were stable, straight, and uniform cross section. Fifteen reaches 
canals were investigated. Values varied from cfs 1,500 cfs and 
slopes from 0.79 10-3 0.97 10-2. 

The primary purpose using these data increase the range condi- 
tions considered and establish more points for the correlation phases 
this study. 

India Canal Data.—A rather thorough study the available literature 
regime theory, its conception, and evolution was undertaken. During the course 
this investigation, considerable information was found the canal systems 
India which was used help establish and develop the regime theory. Two 
groups these data were sufficiently complete andpertinent warrant inclu- 
sion for use the theoretical analysis. 

The first group the data for forty-two stable Punjab canals (reference 
(18), pages through 64). Their capacity varies from cfs 9,000 cfs. 
Slopes are the order 0.12 10-3 0.34 10-3, andtheaverage diameter 
the bed material approximately 0.43 mm. 

insight the magnitude sediment load carried the Punjab canals 
was obtained (reference (17), page 87). The mean silt intensity for the chan- 
nels listed 0.238 per lor 238 ppm. This value within the order mag- 
nitude sediment concentration inthe canals measured Simons and Bender. 
Considering the small difference magnitude sediment load for these two 
groups canals, anticipated that they will behave single group ex- 
cept possibly where major differences bed and bank conditions exist. 

The second group canal data was found statement (reference (10), 
pages and 71) and statement (reference (11), pages and 75). These 
data were collected different reaches canals, 
and according the foregoing references are stable. Their capacity ranges 
from 311 cfs 9,057 cfs, slopes vary between 0.0592 10-3 and 0.0995 10-3 
and mean size bed material within the limits 0.0346 0.1642 mm. 

The Sind canals seem carry larger amount sediment than those 
the Punjab, least during part each year. This verified studying the 
data accordance with observations and taken 1934, the magni- 
tude the suspended sediment load ranges from 3.59 down 0.156 
per from 3,590 ppm down 156 ppm. The silt the Sind canals has 
smaller mean diameter than that found the Punjab canals. With the some- 
what larger sediment load, anticipated that these canals will behave dif- 
ferently from the Punjab canals canals studied Simons and Bender 
unless sufficient wash load occurs the Sind group automatically compen- 
sate for the difference conditions. 

Imperial Valley Canal Data.—Data from four irrigation canals the Im- 
perial Valley canal systems were obtained from atechnical bulletin Samuel 
unique that their sediment concentrations are relatively high, ranging from 
2,500 ppm 8,000 ppm (most this concentration must wash load), and 
their bed and bank conditions are similar those found the Punjab canals, 
the Sind canals, and the canals investigated the writer. 

Water Temperature Data.—An average effective water temperature for the 
twenty-four canal reaches investigated the writer has been worked out for 
each canal. These data are given Table 
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TABLE INFORMATION CANALS 


Maximum 

Sustained 

Location Cubic 

Feet per 

Second 

Bijou, West Morgan, Colo, 
T3N, R59W, S10 


Farmers miles east Mitchel, Neb. 950 
T23N, 

Ft. Laramie West Torrington, Wyo. 1125 
Mile 31.65 

Laramie South Lyman, 475 
Below mile 91,2 

Laramie South Lyman, 530 
mile 87.18 

Ft. Laramie Southwest Torrington, Wyo. 1030 
mile 

Ft. Morgan West Ft. Morgan, Colo, 137 
T3N, R58W, 

Ft. Morgan miles west Ft, Morgan, Colo, 190 
R58W, S18 and 

Ft. Morgan III mile west Ft. Morgan, 160 
T4N, R58W, S36, R55W, S30 

Ft. Morgan West Ft. Morgan 170 
T4N, R58W, sec 

Ft. Morgan West Ft. Morgan 200 
T4N, R59W, 

Garland miles west Powell, Wyo. 880 
14, T54N, R100W, S17 

Garland 4.5 miles Powell, Wyo. 750 
Parallel H.W. 14, 

Interstate North Morril, 1200 
mile 

Larimer Weld Immediately west intersection 600 


the Canal and H.W. west 
Ft. Collins, 


Lucerne miles west Torrington, Wyo, 
Barnes Siding, T26N, R63W, S32 

Lucerne mile west Barnes Siding 
Parallel 26, T26N, 
R63W, S30 

Parallel 26, T25N, 
R62W, S36 

Bijou, West Morgan, Colo, 190 
T3N, S10 

CNPP and 0.7 miles Funk, 240 
below mile 33.8 

Lat 0.5 miles Holdrege, Neb. 110 
mile 

Cozad South Gothenburg, 240 
R25W, S17 

Dawson South Cozad, Neb. 360 
R23W, 

Taylor Ord West Taylor, Neb. 180 


mile 3.4 


Discharge, 


Discharge, 
Cubic 
Feet per 

Second 


199 
236 
113 
227 


363 


181 Heavy 


Extent 
Bank 
Vegetation 


Moderate 


Light 


Light 
Light 
Light 
Light 
Moderate 
Light 
Light 
Light 
Light 
Light 


Light 


Light 


Light 


Moderate 


Moderate 
Heavy 


Moderate 
Heavy 
Heavy 


Moderate 


Light 


Stable, except where otherwise None all cases. 


| 
137 

| 


WATER CONVEYANCE CHANNELS 


INVESTIGATED SIMONS AND BENDER 


Bank 
Material 


T 


Maintenance 
Required> 


Bern 
Formation 


Method Operation 
During Season 


No, 


Moderately cohesive 
Moderately cohesive 
Cohesive 

Cohesive 

Cohesive 

Cohesive 

Cohesive 
Non-cohesive 
Non-cohesive 
Non-cohesive 
Non-cohesive 
Gravel 


Gravel 


Moderately cohesive 


Moderately cohesive 


Moderately cohesive 


Moderately cohesive 


Moderately cohesive 


Moderately cohesive 
Cohesive 

Cohesive 
Non-cohesive 
Non-cohesive 


Non-cohesive 


Marginal 


Marginal 


bank 
None 
Very little 
Very little 
Very little 
None 
Very little 
None 
Very 
None 
None 

None 


None 


None 


None 


Very little 


Very little 


Very little 


bank 
Very little 
Very little 
None 

None 


Very little 


Intermittent 
Continuous 
Continuous 
Continuous 
Continuous 
Continuous 

Slightly 
Slightly 
Slightly 
Slightly 
Slightly 
Continuous 


Continuous 


Continuous 


Intermittent 


Continuous 


Continuous 


Continuous 


Intermittent 
Continuous 
Continuous 
Continuous 
Continuous 


Continuous 
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TABLE 3.—SUMMARY 


Water 
Tem- 
perature 
F° 
2.42 2.92 79.0 
2.48 5.85 69.0 
8.29 72.0 
1,92 6.01 73.8 
2.10 73.0 
1.79 74.0 
3.51 77.0 
2.63 77.0 
1.39 2.93 77.0 
2.91 79.0 
79.0 
2.26 7.88 62.0 
2.57 5.73 61.5 
2.51 8.50 
2.55 4,90 63.0 
1,84 2.61 69.0 
3.01 71.0 
1,56 2.64 70.0 
2.42 3.31 73.0 
5.25 78.0 
2.01 4.33 79.0 
3.33 82.0 
1,90 2.95 83.0 
1,86 3.67 82.0 


Based the analysis one large sample. 


The temperature the San Luis Valley canals not known, but, based 
climatological conditions, effective temperature has been assumed 
for each the fifteen reaches reported. 

specific information temperature individual canals the Punjab 
the Sind was found. has been pointed out (reference (6) page 55), however, 
that temperature variation the Punjab canals ranges from and 
that good average for the entire year. addition, Blench (3) reports 
that the climate the Punjab similar that desert Arizona and that wa- 
ter temperatures vary from Based the foregoing information, 
base temperature was utilizedwhenever temperatures water were 
involved forthe Punjab andSind canals. This base temperature probably intro- 
duces some scatter those relationships which involved. 


RELATIONSHIP BETWEEN SEDIMENT LOAD, 
FROUDE NUMBER, AND CHANNEL STABILITY 


understand the importance sediment load and Froude number onchan- 
nel stability necessary consider the behavior alluvial channels 
complete detail. 

Regimes Flow and Forms Bed Roughness.—In alluvial channel the 
forms bed roughness are function the bed material, the sediment 


No. 
113.0 56.3 


WATER CONVEYANCE CHANNELS 
SIMONS AND BENDER DATA 


Estimated 
Dune 


Height 


0.0515 
gravel 
0.0806 
0.098 
0.143 
0.166 
0.109 
0.060 
0.149 
0.074 
0.182 
0.286 
0.036 
0.034 
0.177 
0.271 
0.067 


to & 
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transport, and the flow. Based the foregoing field study and laboratory 
study alluvial following regimes flow andforms bed rough- 
ness have been observed. 

Tranquil flow regime (Fy <1) 


Plane bed without movement bed material; 
Ripples; 

Ripples superposed dunes; 

Dunes; 

Transition from dunes rapid flow; 

Plane bed with movement bed material; and 


Rapid flow regime <1); 


Standing water waves and sand waves; and 
Antidunes. 


These forms bed roughness have been described Albertson, 
Simons, and Richardson (1). 

Sediment Load.—The magnitude the total silt-sand sediment load varies 
with form bed roughness, accordance with laboratory data (1), and ac- 


0.022 0,096 0.274 0.017 0.161 443.0 
0.161 0.465 0.945 0.479 40.30 28.8 
0.015 0.158 0.014 0.069 13.38 


cordance withdata from studies, the relationships between total sand- 
silt transport and form bed roughness is, for sand beds well graded ma- 
terial having median diameter 0.2 mm, roughly follows: 


Flow Regime and Form 
Bed Roughness Sediment Load, ppm 


Tranquil flow regime 


Ripples 
Dunes 1,000 
Transition dunes 

rapid flow 1,000 3,000 
Plane bed with movement 

bed material 2,000 5,000 


Rapid flow regime 


Standing waves 


should emphasized that these limits total load are qualitative and will 
influenced many factors—such asthe characteristics material, 
channel geometry, shear the bed, and characteristics the suspended sedi- 
ment. For example, considering Kalinske and Hsia data (1), the total load was 
the order 17,000 ppm with dunes. However, this case, the median di- 
ameter mm. Very fine suspended sediment which sometimes re- 
ferred wash load has been defined that material transported the 
flow which not found appreciable quantities inthe bed. result field 
that channel stability very intimately related wash load when wash load 
concentrations are large. The foregoing figures illustrate that with the ripple 
and dune forms bed roughness the presence very fine sediment reduces 
resistance flow and, also, drastically influences channel geometry. The ef- 
fect fine sediment fluvial mechanics sufficiently important warrant 
further study. 

Channel Stability Related Froude The field study stable 
channels conducted the writers and subsequent discussions channel sta- 
bility with others has emphasized that channel stable (no appreci- 
able bed and/or bank scour accretion occurring with time) the Froude num- 
ber defined must most cases less than 0.3 for alluvial mate- 
rial the sand-size range and finer. When 0.3, the banks scour the 
straight reaches channel and the bends must stabilizedto confinethem 
their rights-of-way. 

Summary.—Summing the preceding concepts they apply design 
stable channels, can concluded, for alluvial materials the sand range 
and finer, that excessive bank erosion avoided the Froude number 
must not exceed 0.30. This immediately confines the design problem the 
tranquil flow regime and excludes the transition and plane bed form bed 
roughness within this regime because the magnitude the Froude number 
associated with it. That is, the designer restricted such that the channel 
should designed with the ripple, dune form bed roughness, depending 
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channel slope limitations dictated the terrain, and the magnitude 
sand load which the channel must transport. Once recognized that de- 
signs must confined range small Froude numbers and the forms 
bed roughness cited, immediately apparent that the magnitude sand- 
silt sediment load which can transported is, likewise, severely limited— 
usually less than 500 ppm. 

This discussion brings out still another important point. Canals must uti- 
lize head works structures which not allow sediment concentrations, the 
sand-silt range, exceed 500 ppm for the dune form bed roughness, and 
the concentration sediment must even smaller for ripples. Otherwise 
the channel will unstable—aggrading with time. 

channels are designed carry larger sediment loads (concentrations 
500 ppm) the channel will unstable because banks and bend erosion 
0.3) unless some form bank stabilization utilized. cannot 
overemphasized that one must transport more sediment than can handled 
within the dune regime will necessary provide bank stabilization such 
gravel rip-rap, some other form armor plating, masonry con- 
crete. The tractive force method design, presented Lane (13), pro- 
vides satisfactory means designing stable protection, using coarse non- 
cohesive material for this condition. 


ANALYSIS DATA 


The tabulated basic data and the parameters computed therefrom can now 
utilized investigate the theory stable channels. The shape character- 


istics these canals will investigated first. 

Relationship Between and and and W.—The Lacey theory expressed 
terms wetted perimeter and hydraulic radius. The Blench theory 
terms average depth the channel bed and average width. 

Considering the Punjab canals the average depth the bed and average 
width were not given. The only measurements pertaining shape were wetted 
perimeter, hydraulics radius, and top width. overcome this deficiency 
data, and for design purposes, hydraulic radii were correlated with average 
bed depths and wetted perimeters were correlated with average width chan- 
nel, shown Figs. and Data from the canals investigated the 
writer, the Sind canals, and the Imperial Valley canals were utilized es- 
tablish these curves. 

The Imperial Valley canal data have been plotted Figs. and show 
that increased charge has little effect the relationships between and 
and and 

Based Figs. and now obvious that average width and bed depth 
can estimated rather accurately and are known. 

The actual method used compute and for the Punjab canals in- 
volved estimating knowing and then computing average bed depth means 

expected that either wetted perimeter some channel width dimension, such 
top width average width, should correlate with rate discharge. Based 
this type correlation, Lacey arrived Eq. 

Estimating P.—Using Lacey’s procedure, values and corresponding 
values were plotted Fig. first for the canals investigated Simons 
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and Bender (referred Simons and Bender data). The range materials 
forming the periphery these canals extends from fine cohesive material 
coarse non-cohesive material. The effect soil type clearly exhibited 
Fig. The sand channels all require relatively large for given 
while cohesive materials reach stability relatively small for given 

illustrate more fully how the Simons and Bender data plot relative 
the India data, values versus have also been plotted Fig. for the 
forty-two Punjab canals. Three curves have been fitted these data based 
type bed and bank material. The equation the arbitrary straight line 
representing canals with sand beds and cohesive banks 


Note the similarity Eqs. and 19. 

The effect bank and bed materials such relationship illustrated 
further considering the USBR data. With these data seen that values 
for given are even smaller than the case cohesive materials. 
This simply illustrates the ability the coarse material resist greater 
tractive force. 

The Imperial Valley canal data have also been plotted illustrate the quan- 
titative effect increased charge. From the viewpoint type bank and bed 
materials, these canals are similartothe other canals, excluding those formed 
coarse non-cohesive material. The trend line, however, falls approximate- 
the trend line representing the relationship between and for canals 
formed coarse non-cohesive material. This indicates that fine sedi- 
ment load, the type found the Imperial canals decreases, the stable wetted 
perimeter, and consequently stable width, decreases. 

Relationship Between Average Width and Top Width.—In most the rela- 
tionships involving width, the average value has been used. Under certain 
circumstances, may advantageous convert average width, top 
facilitate this conversion, see Fig. This relates with WT, 
based the data from the canals investigated Simons and Bender and the 
Punjab canal data. good straight-line correlation exists between these vari- 
ables top width approximately 300.0 ft. However, for the larger val- 


ues the relationship based onlimited data. The equation relating these 
variables 


Estimating R.—A relationship between and for stable irrigation canals 
was given Fig. now desirable relate other quantities pa- 
rameters. Lindley showed that both and were very closely related 
for all canals have been utilized, the foregoing relation, between and 

The effect natural bed and bank material apparent and consistent with 
the preceding relations. The canals constructed cohesive materials fall 
one group, those having sand bed and natural berm second group, those 
sand material third group, and those coarse non-cohesive material 
fall fourth group. 
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The basic relationship relating and for canals having sand bed and 
natural berm 


The same type relationship for the canals coarse non-cohesive material 


The influence sand banks has not been expressed equation form but 
their effects are clearly illustrated. 

The Imperial Valley canal data have been plotted illustrate the impor- 
tance charge stability channels (Fig. 5). Comparing the trend line for 
these limited data with similar lines representing the Simons and Bender, 
Punjab, and Sind data, apparent that the hydraulic radius, decreases 
with increased charge for given discharge, 

Relationship Between and interest note the close relation- 
ship existing between and for stable channels. Once again, however, 
expedient introduce type bank material third variable (Fig. 6). Val- 
ues and been plotted from all canals. Four curves drawn, 
indicated bank material and bed condition. The short, uppermost curve 
for canals having sand The intermediate curve for all other bank 
materials finer than sand. The lower curve for coarse non-cohesive mate- 
rials, represented the USBR canal data. 

Equations relating and for two the foregoing three conditions follow. 
For banks material finer than sand extending into the very cohesive range 


and for coarse non-cohesive materials 


The relationship between and for the Imperial Valley canal data shows 
that magnitude charge increases the required area, corresponding 
given discharge, decreases and the permissible average velocity, 
increased. This consistent with the trends indicated Figs. and 
also important note this case, the variation vertical displacement 
each point relative the line representing sand beds and cohesive banks. The 
effect variation charge shows very clearly. These expressions cover 
wide range discharge and boundary conditions and should very useful 
practical design work. 

Expressions Involving Velocity, Hydraulic Radius and Slope.—A variety 
regime-type equations have been, recommended determine design slope 
one individual another. One the most significant the Lacey formula 
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The general validity this expression shown Fig. extremely 
wide range discharge covered. Data used include the Punjab data and 
other miscellaneous India canal data which are not presented utilized else- 
where report. The difficulty with this expression that, although the 
trend very definite, some slopes computed this relationship vary 
siderably from the measured slopes used establish the correlation. 

Fig. for all the canals included this report. These data fall three sep- 
arate groups the plot, line has been drawn through each group. The 
upper line correlates with for the coarse non-cohesive materials rep- 
resented the USBR canal data. The intermediate line does the same for 
canals with sand beds and slightly cohesive cohesive banks. The third line 
represents sand bed conditions, least insofar asthe Simons and Bender data 
are concerned. The equations these lines are: 


respectively. Eq. identical with the preceding Eq. corresponding 
the plot Fig. 

The Imperial Valley canal data show that for given value the per- 
missible average velocity, increases charge increases and knowing and 
from Figs. and the slope consistent with stability can estimated. 


The Blench-King Regime Slope regime slope equation 
mended Blench for design 


5/6 
(29) 
2080 


This was derived plotting W/D against variety non-dimensional groups 
and would have been found once, according Blench, plotting V2/g 
against W/v. The basic regime slope formula 


1/4 
(30) 


where the average channel width. 

Fig. values V2/g have been plotted against W/v. value 
corresponding has been assumed for all India canals. The Punjab ca- 
nals yield points that plot quite close straight line log-log paper be- 
tween the limits Beyond this upper limit, the V2/g 
terms are nearly constant andthe slope theline flattens until lies approx- 
imately parallel the horizontal axis. 
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The twenty-eight points corresponding the sind data lie more less 
extension the straight-line portion the Punjab data. 

The Simons and Bender data give points that generally intermingle with the 
India canal data, except that seven the points corresponding canals with 
non-cohesive banks fall lower, one would anticipate. 

The significance the basic regime slope equation, presented King and 
Blench, quite apparent from study Fig. The India data group plot 
close straight line between the limits 105 Beyond 
more points fall below the straight line than above it, however, 
this may function the canals sampled. 

Tractive-Force Relationships.—The various procedures used estimate 
magnitude tractive force the bed and sides channels were discussed. 
Values shear the channel periphery were computed for the canals in- 
vestigated Simons and Bender, using all the methods described. sum- 
mary these data given Table 

Magnitude Tractive Force.—The magnitude boundary shear varies 
with method computation. Fig. 10, typical canal cross section given, 
including the shear distribution its boundary, indicated the various 
methods computation, size side and bed material, and standard deviation 
side and bed material. illustrated Fig. 10, shears computed the 
various methods for given channel are means close agreement. 
general, shears computed Eq. 16, and shears computed use isovels 
and the concept zero momentum transfer are both larger than shears indi- 
cated the velocity gradients measured normal the boundary across the 
channel. The lack agreement results, and the fact that shears based 
velocity gradients are smaller, leads one believe that something being 
neglected when shears are computed based the latter method. may that 
this results because the energy required totransport the sediment load and/or 
the energy involved secondary circulation are being neglected. 

Another interesting observation based shears computed from knowledge 
velocity distributions the way shears vary across the bed section each 
canal. Generally, one might expect more uniform shear distribution such 
may intimately related secondary circulation the has been 
proposed that secondary circulation studied the canals investigated 
Simons and Bender some future dateto provide knowledge its ef- 
fect shear and shear distribution, sediment transport, and channel stability 
general. 

Correlation Tractive Force With Mean Diameter Bed Material.— Lane 
(14), has developed relationship which relates size bed material and trac- 
tive force. This correlation provides very useful means establishing the 
design slope channels and canals coarse non-cohesive materials, pro- 
vided size bed and bank material canbe estimated with 

working interms tractive force based and/or the effect the 
Punjab and Sind canal data relations can shown. Fig. 11, 
have been plotted for all the canals involved inthis study, that is, the Simons 
and Bender data, the USBR data and the India data. 

Several facts interest are immediately apparent First gen- 
eral line extending through all the data can drawn. There is, however, 
considerable scatter about this line. Next, secondary lines crossing the major 
trend line have been drawn based intimate knowledge the Simons and 
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FIG, 11.—VARIATION TRACTIVE FORCE WITH BED MATERIAL 
TYPE CHANNEL AND DISCHARGE 
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Bender data and alimited knowledge the India data from study the liter- 
ature and existing data. Moving the upward direction, these lines 
associated with canals having cohesive beds andbanks, the second with canals 
having fine sand beds and probably berm banks natural banks cohesive 
nature, thethird with canals having coarser sand beds banks banks 
slightly cohesive natural materials, the fourth with canals having sand beds 
and banks, and the fifth with coarse non-cohesive beds and banks. Roughness 
bed seemsto increasetraveling from the bottom secondary line the fourth 
secondary line associated with sand beds and 

Next consider each the five secondary lines. Moving along these lines 
the direction increasing shear, that canal capacity increases. The 
points the extreme right end the secondary lines correspond large 
values, and the points the extreme left these same lines correspond 
small Q-values. The whole system lines shown Fig. are placed rather 
arbitrarily and would undoubtedly shift slightly additional new data were in- 
corporated into the plot. 

The Imperial Valley canal data have not been plotted Fig. because 
uncertainty regarding the mean size bed material. currently be- 
ing made secure these data, since will importance reflect the ef- 
fect increasing the magnitude charge permissible tractive force. 

Considering the lowest element the major curve, noted that 
were curved the right would fit the plotted points somewhat better. This 
indicates that allowable tractive force probably increases with decrease 
mean size sediment smaller than that size where the material starts be- 
come plastic. This aspect the problem was investigated flume study 
the University Wyoming (23). study verified that limiting tractive 
force increases with increase plastic index. The investigation was carried 
out under the supervision the writer using natural materials from canal beds 
investigated Simons and Bender. 

Channel Shape.—Considering the canals investigated Simons and Bender 
was observed that channel shapes varied widely. was also apparent that 
shape was appreciably affected thetype natural bank material, the amount 
and type bank vegetation, and type and concentration wash load. Most 
these canals, excluding those found sand, had sides that were very 
the water surface (this was possible because the reinforcement provided 
the plant roots) and that were asymptotic the channel bed. 

The shape these canals, with some exceptions, conform reasonably 
well the theoretical regime channel shapes described King (24) and 
Glover (8). 

Lane (13), and others have pointed out the possibility designing channels 
such shape and dimensions that the entire wetted perimeter ina state 
incipient motion. None the canals investigated have adjusted the forego- 
ing shape, implying that perhaps necessary construct this form in- 
itially desired. seems that this aspect design worthy more 
thorough investigation. 

Angle Repose for Non-cohesive conjunction with selecting 
suitable side slopes when designing channels ripraping canal banks, 
necessary able estimatethe angle repose the bank and/or rip-rap 
material. Lane (14), presented figure relating angle repose, median di- 
ameter the material and shape particles. Subsequent investigation this 
relationship Lane and Simons (25) verified that additional information 
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angle repose was needed. Asa result further study, Fig. was devel- 
oped. The angle repose, indicated this figure, for dry saturated 
material. the material saturated and subjected the action flowing 


water well, such canal banks are, the side slope should reduced below 
that indicated for static conditions 10°. 


DESIGN PROCEDURES 


The relationships the previous paragraphs, regardless their limita- 
tions, provides useful design procedure. The objective this section 
discuss the procedure and point out the respective limitations and advan- 
tages. 

most cases, because the general scatter data, the complexity the 
relations, and the time involved, curves have been fitted the data 
Where the scatter the points about the trend line not excessive, equations 
describing the relations have been determined. 

Selection P.—In Fig. and type and condition bed and bank ma- 


terial have been related. The four parallel curves presented are representa- 
tive of: 


Sand beds and banks; 


Sand beds and slightly cohesive cohesive banks; 
Cohesive beds and banks; and 
Coarse non-cohesive beds and banks. 


Although data for canals having sand banks are limited, Curve provides 
means estimating for this condition except possibly the very small, 
and, also, the very large, ranges Curve valid for the complete 
range covered the basic data and has been described mathematically 
Eq. 19. 

based very limited data and should employed with this fact 
mind. 

Curve4is recommended means estimating the magnitude the 
coarse non-cohesive range bed and bank materials. 

this point, knowing the wetted perimeter, the average stable channel width 
can determined from Fig. 

Estimating Knowing Vice Versa.—Wetter perimeter and top width 
are closely related, illustrated Fig. Based this relation, either 
can estimated for stable canals provided one the other known. 
The wetted perimeter, would undoubtedly correlate equally well with top 
The wetter perimeter estimated this manner should more repre- 
sentative true conditions than were computed based some initial 
trapezoidal shape. 

Selection R.—Referring the relationships between discharge and hy- 
from the appropriate curve can computed equation the sand 
bed-cohesive bank coarse non-cohesive range materials. These equations 
are Eq. for sand beds and cohesive banks, and Eq. for coarse non-cohe- 
sive materials. 


Equations were not developed forthe curves because these trends 
are based rather limited information. 
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Rate discharge could correlated with average depth bed depth, 
the need developed. However, the writers feel that hydraulic radius more 
meaningful measure for design purposes. 

Estimating When Known Vice Versa.—Corresponding values 
and are closely related one another stable canals, illustrated Fig. 
this case, the average bed depth. similar correlation could 
established relating average depth hydraulic radius average depth could 
correlated directly with bed depth. 

Selection Based and Soil Type.—The value Acan determined 
from the curves Fig. knowing and soil type, can solved directly, 
Eqs. 24(which depend soil type) apply. before, because limit- 
data, equation was writtenfor canals sand beds and sand banks. 

The area required transport given discharge maximum for sand 
banks and beds, somewhat less for slightly cohesive cohesive banks inclu- 
sive and minimum for coarse non-cohesive banks and beds. The explanation 
the preceding primarily difference stability bank ma- 
terials. 

Determination Design Slope.—As shown the preceding paragraph, the 
average velocity can estimated from soil type, and Knowing and 
possible evaluate slope referring Fig. which correlates and 
That is, for coarse non-cohesive materials Eq. applicable; for 
canals with sand beds and cohesive banks Eq. applies (Eq. the same 
Eq. suggested Lacey); and Eq. used for canals with sand beds 
and banks. 

Slope Determination the Blench-King Slope Formula.—The Blench-King 
slope formula, when applied dimensionless form the Simons and Bender 
and the India canal data, yielded Fig. The correlation quite significant 


for values <8,000,000, particularly for the India data. Beyond 


this point seems logical use some other method estimateS--since even 
the India data scatter badly when 107. 

evaluate slope, first estimate and Next compute the parame- 
ter W/vusing appropriate value for the viscosity. Then enter Fig. and 
obtain value for Knowing and the slope can determined. 

Determination the Tractive-Force Method.—The scope the origi- 
nal versus relation(14)has been broadened include conditions encoun- 
tered the canals investigated Simons and Bender and the India canals 
(Fig. 11). Using this family curves estimate can made first 
estimating mean size bed material, hydraulic radius, type bank conditions 
that will result, and knowing 

The major limitations this method are lack knowledge size bed 
material and the scatter occurring the plot—which is, course, related 
the accuracy with which slope can estimated. 

Summary.—Asummary the range the more important variables utilized 
establish the relationships which have been presented follows: 


varies from 9,000 cfs 
Slope varies from 0.000058 0.000388 
Average width varies from 264 
Depth varies from 2.8 10.5 
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Sediment concentration varies from ppm 500 ppm excluding four ca- 


nals which have concentrations ranging from 2,500 8,000 ppm the larger 
part which must wash load. 


conclusion, may stated that within the limits the canal data pre- 
sented: 


Both and can estimated with fair accuracy using Figs. and 
and/or the corresponding equations. 
and and/or the corresponding equations. 
reasonable estimate design slope can obtained using one 
more the following: 
Fig. representing the Blench-King regime equation 
form. 
Fig. relating tractive force, mean diameter bed material, type 
bank material, and discharge. 


The relative merits the respective methods available help estimate 
design slope can comprehended and appreciated more fully referring 
actual design problem. 


SUMMARY AND CONCLUSIONS 


Two basically theories are currently recognized because their 
superiority over other available existing methods used approximate the de- 
sign stable channels. These concepts are: 


The regime theory India developed Kennedy, Lindley, Lacey, 
Bose, Blench, and others, and 


The limiting-tractive-force theory proposed Lane and others, 


The applicability the regime theory has been emphasized this paper. 
investigation the regime theory verifies that the regime equations 


India are only valid for the limited range conditions upon which they are 
based follows: 


Channels having sand beds, and slightly cohesive cohesive banks, the 
banks which are usually formed the berming action the suspended sedi- 
ment. 

Channels that are not required carry heavy charge sediment for 
sustained periods time. That is, the canals yielding the data upon which the 
India regimetheory havetheir magnitude charge controlled sedi- 
ment exclusion and/or ejection structures that is, generally, less than 
500 ppm, which about the upper limit sediment load (sand and silt size) 
that can transported without appreciable bank erosion. 


The range conditions which this theory applies has been expanded 


result this study that canals each the following groups can de- 
signed this method. 
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(1). Canals formed coarse non-cohesive material the type studied 
the USBR (14) (charge <500 ppm). 

(2). Canals formed sandy material with sand beds and banks (charge 
500 ppm). 

(3). Canals possessing sand beds and slightly cohesive cohesive banks 
(Good results when charge <500 ppm, qualitative results when charge 500 
ppm). 

(4). Canals having cohesive beds and banks (charge <500 ppm). 


Within each the preceding four classifications possible evaluate 
area, average width, top width, bed depth, average depth, hydraulic radius, 
wetted perimeter, average velocity, and side slopes with ease and practical 
degree accuracy. This accomplished means the regime-type cor- 
relations involving these parameters, previously discussed (Figs. 11). 

With the exception Fig. none the slope equations relations inves- 
tigated reflect the effect sediment load slope. The major reason why the 
Simons and Bender data, well other data, have not cast more light this 
problem that the sand-silt sediment load measured canals proved 
fairly constant varying only within very narrow limits (concentration 500 
ppm)—thereby making virtually impossible determine any meaningful ef- 
fect. 

Utilizing all the basic data, attempt extend the range applicability 
the tractive force method determining slope include all canals 
was investigated. The results this study well summarized Fig. 
11. terms this immediately obvious that canals formed materials 
the sand range andfiner constitutes group within whjch sig- 
nificantly more complex than the coarse non-cohesive range. Note that five 
different curves have been drawn, each representative adifferent sub-group 
canals, and that increases moving from left right along any curve. Es- 
timating adesign slope inthe range sand andfiner means these curves 
sense precise approach, but the figure extremely useful that 
within limits slopes can estimated, and the presentation gives insight 
the complexity conditions withinthis range unrecog- 
nized the tractive-force 

Summarizing, the tractive-force method design seems valid for the 
coarse non-cohesive range. However, recommended that equations the 
regime type relating such terms and and and should used 
estimate width-to-depth ratio preference arbitrarily selecting and 
and then computing slope consistent with stability sides and bed, unless 
working with coarse material stabilization problem where coarse ma- 
terials will used armor plate the channel banks—the natural bank mate- 
rial being smaller 

Because the less significant correlation between tractive force, mean 
diameter bed material, etc., depicted Fig. 11, within the sand and co- 
hesive range particle sizes, other relations the regime type are perhaps 
equal superior the tractive force method for estimating design slope. 
any event, seems logical, based the validity regime type area, width, 
and depth relations, that they should used establish W/D regardless 
whether regime tractive force equations are used estimate design slope. 

The disadvantage imposing arbitrary W/D ratio design problem 
that when narrow widths are selected (that is, small W/D ratio) the mag- 
nitude slope must limited avoid bank erosion. This small 
average velocity and reduced ability transport sediment. 
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The foregoing brings out the advantages strong points the 
two theories, except possibly those cases where coarse non-cohesive mate- 
rials are involved. Then apply the tractive force theory directly. 


Dimensions 
L/T2 
L/T2 


F/L3 
degrees 
L2/T 
FT2/L4 
F/L2 


Definition Description 
Area water cross section 
Bed factor 

Bed width channel 

Chezy coefficient 

Median diameter 

Average bed depth 


Lacy silt factor 


Manning coefficient roughness 


Froude number 
Reynolds number 
Wetted perimeter 
Discharge 
Hydraulic radius 
Side factor 

Slope 

Average velocity 


Average width channel such 
that 


Top width channel the 
water surface 


Specific weight 
Angle repose 
Kinematic viscosity 
Mass density 


Tractive force shear 
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RESISTANCE FLOW ALLUVIAL 


Daryl Simons,! ASCE and Richardson,2 A.M. ASCE 


SYNOPSIS 


This paper presents the initial results flume study alluvial channels. 
detailed classification the regimes flow, the forms bed roughness, 
and the basic concepts pertaining resistance flow are discussed. 


INTRODUCTION 


The problem defining roughness alluvial channels dates back several 
centuries. The solution this problem, least the complete solution, has 
thus far eluded man. The principal reasons why only limited answers, which 
some cases are questionable value, have been developed that the scope 
the problem broad, and that there great number variables which 
influence resistance flow alluvial channels. Roughness alluvial chan- 
nels more complex problem than roughness rigid channels because the 
form the bed function the flow. That is, not only the alluvial chan- 
nel roughness elements resist the flow, but they, inturn, are shaped the 
flow. Thus far, mostscientistsworking this field have been limited time, 
facilities, instrumentation needs, and funds. result, only small parts 
the complete and complex problem have been thoroughly investigated. 


Note.—Discussion open until October 1960. extend the closing date one month, 
written request must filed with the Executive Secretary, ASCE. This paper part 
the copyrighted Journal the Hydraulics Division, Proceedings the American So- 
ciety Civil Engineers, Vol. 86, No. May, 1960. 

Presented the October 1958 ASCE Convention New York, 

Engr., Geological Survey, Fort Collins, Colo. 

Engr., Geological Survey, Fort Collins, Colo. 
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THEORY FLOW ALLUVIAL CHANNELS 


analysis flowin alluvial channels extremely complex because the 
many variables involved and the difficulty measuring them. Simons, 
Richardson and Albertson the pertinent variables and 
important dimensionless parameters related resistance flow alluvial 
channels, utilizing the principle dimensional analysis. 

Dimensional equation relating the selected dependent varia- 
ble and the independent variables states that 


Using and the repeating variables and applying the Pi-theorem 


Vdp 
The parameters involving are the Reynolds number and the 


Froude number respectively. Since for flow water 
nels, the Froude number can expressed the form 


The last term Eq. can modified drag coefficient for the parti- 
cles. That is, 


d Ays 
Eq. can rewritten 


This equation can simplified noting that: 


The width-depth ratio term which probably secondary impor- 
tance provided B/D When the effect side walls may appre- 
ciable. 

The relative standard deviation the size distribution the bed ma- 
terial undoubtedly influences flow under certain conditions which must ulti- 
mately determined. However, this stage the investigation only one 
sand has been used, eliminating variable. 

The shape factors for the channel and may eliminated because 
the bed completely alluvial, the cross section uniformly rectangular ex- 
cept for variation due bed roughness, and the channel straight. 


Utilizing these simplifications, Eq. reduces 


ardson and Albertson, Geol. Survey Water Supply Paper 1948A, 1959. 
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The parameters Eq. can modified the same manner 
that the last term Eq. was modified drag coefficient 


EXPERIMENTAL EQUIPMENT AND PROCEDURE 


This study (by the United States Geological Survey (USGS) the Colorado 
State University) roughness alluvial channels was conducted tilting 
recirculating flume 150 long, wide, and deep Any slope be- 
tween the limits 0.013 per could set and discharge could var- 
ied from cfs cfs. 

bed 0.7 deep was placed the flume. The sand was natural riv- 
sand obtained from commercial sand and gravel company located the 
Cache Poudre River. The median fall diameter the sand, determined 
using the visual accumulation size analysis was 0.45 mm. 

The size distribution plotted straight line log-normal probability 
paper (Hazen, Whipple and Fuller) with slope 1.60. 

Forty-five runs completed with this sand form the basis this report. 
These runs ranged from the condition sediment movement antidune 
flow. the writers’ knowledge, this the first set data, excluding that 
Gilbert,5 which covers such broad range flow conditions. 

General Procedure.—The general procedure followed for each run involved 
recirculating given discharge the water-sediment mixture the flume 
given slope until equilibrium conditions were established. 

Slope selection was accomplished general sense. any flow system 
where discharge, depth, and slope can varied, only two the three varia- 
bles can considered independent. natural stream the discharge and 
slope are normally independent with depth dependent. the flume the dis- 
charge was independent, slope was independent within limits, and depth was 
dependent. The slope was preset the beginning run adjusting the tail 
gate. Such adjustment indirectly influenced the depth dependent variable. 
Generally, and especially the flatter slopes (0.0014 0.006), the slope 
the water surface was adjusted parallel with the bed. With the development 
bed configuration, the slope and depth adjusted the new condition bed 
roughness. Thus, for these experiments the slope and depth are function 
and the roughness which develops regime flow. The non-uniformity 
flow caused change bed roughness was eliminated continuing the 
run until the bed slope and the slope the water surface again became parallel 
each other natural adjustment the sand bed. 

the tranquil flow regime the slope each run was adjusted the begin- 
ning first setting backwater curve, and then alternately adjusting 
the tailgate and measuring the water surface slope until reached the desired 
value. For the flat slopes (0.00014 0.0006), the bed was carefully screeded 
the desired slope before the run. With the steep slopes the bed was only 
raked remove the dunes from the previous run before starting the new run 


Accumulation Tube for Size Analysis Sand,” Colby and 
Christensen, Proceedings, ASCE, Vol. 82, No. June, 1956. 

Prof. Paper 86, 1914. 
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because the increased sand movement assisted obtaining equilibrium condi- 
tions quickly. the rapid-flow regime the slope was changed altering the 
bed control the lower end the flume and sediment movement quickly es- 
tablished equilibrium conditions. 

Equilibrium flow was considered established if: 


The same bed configuration existed over the full length the flume, ex- 
cluding the section influenced entrance conditions. 


The water surface slope remained essentially constant with respect 


The period time required establish equilibrium conditions varied with 
the slope and the discharge. Some runs with flat slopes required and days 
achieve equilibrium, whereas, for runs with steep slopes, equilibrium was 
established 3hr. Every run involved continuous flume operation until 
was completed. Whether not equilibrium conditions were established was 
based the measured data and the judgement the experimenter. assure 
the achievement equilibrium conditions, most the experiments were con- 
tinued longer than the required time indicated measurements. 

After equilibrium conditions were established, water-surface slope, dis- 
charge, water temperature, depth, velocity profiles, total sediment concentra- 
tion, and suspended sediment concentration were measured. Water-surface 
slope was measured with Lory point gage and precision level determining 
water-surface elevations every along the flume. Discharge was measured 
with calibrated orifice meters. Depth was determined subtracting mean bed 
elevation from mean water surface elevation. Velocity measured 
three verticals with calibrated Prandtl pitot tube. Total sediment concen- 
tration was measured traversing with width-depth integrating sampler 
through the nappe the downstream end the flume. Suspended sediment 
concentration was measured with DH-48 hand sampler. After the run was 
completed and the flume drained, profile the bed configuration was ob- 
tained and bed material was sampled. The basic data are presented Table 


OBSERVED FLOW PHENOMENA 


The form bed roughness alluvial channels function the sediment 
characteristics and the characteristics the flow. That is, the bed configu- 
ration can changed altering either discharge (which affects the depth), 
slope, temperature, the median diameter size distribution the bed ma- 
terial. the flume experiments the bed material was not changed and the 
slope, temperature, and depth were varied. Under these conditions the follow- 
ing forms bed roughness were observed. 


Tranquil flow regime, 
Plane bed without movement. 
Ripples. 
Dunes with ripples superposed. 
Dunes. 
Transition from dunes plane bed. 
Rapid flow regime, 
Plane bed and water surface. 
Standing waves. 
Antidunes. 


TABLE 1.—BASIC VARIABLES FOR 0.45 BED MATERIAL 


. 


No. ft/ft fps ppm Roughness 
3.94 0.81 Plane 
.016 6.22 0.80 Ripples 
.017 §.11 0.79 Ripples 
.019 1,84 0.65 Plane 
5.07 0.79 Ripples 
3.62 0.78 Ripples 

7.90 1.20 Ripples 

3.84 0.88 Ripples 

7.85 1.23 Ripples 

7.94 1.44 Dunes 

3.83 0.93 Ripples 

7.98 1.43 Dunes 
.088 1.95 0.75 Ripples 

3.90 1.07 Ripples 
4.24 1.30 Dunes 
12.12 1.58 Dunes 
13.54 Dunes 
.189 4.91 Dunes 
8.14 1.68 Dunes 
13.34 2.57 Dunes 
8.73 1.76 Dunes 
.301 21.41 3.32 Dunes 
.364 20.64 4.71 Transition 
11.19 3.52 Transition 
.369 4.54 1.89 Transition 
7.91 2.99 Transition 
3.15 2.04 Transition 
21.62 5.05 Transition 
5.33 Transition 
5.58 Transition 
8.44 Transition 
Transition 
21.38 Transition 
Transition 
656 14.96 Antidune 
862 5.58 Antidune 
898 10.83 Antidune 
3.43 Antidune 
1,01 21.42 Antidune 
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Typical ripple pattern 


Typical dune 


Typical standing waves 


Typical antidune pattern 


FIG, FORMS BED ROUGHNESS 
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The major forms bed roughness are sketched Fig. 

Slope and bed material were the dominant factors which determined the form 
bed Although, when slope and bed material remained constant, 
the form the bed could changed varying depth (changing changing 
temperature. However, changing depthor temperature did not change flow from 
tranquil rapid unless the slope for given bed material was close the 
critical slope. That is, for given bed material, slope and dune-bed 
form roughness, was impossible change standing waves changing 
depth With given depth, temperature, and bed material, all 
bed formscan observed changing slope. This last statement must 
qualified because with shallow depths three-dimensional flow and multiple bed 
roughness forms will occur. also possible have size bed material 
which will not produce some the bed forms regardless the depth, slope, 
water temperature. 

With given bed material, the change from one bed form another not 
necessarily abrupt. Neither does particular form occur the same slope 
(depth and temperature held constant) nor the same depth (slope and tempera- 
ture held constant) the process reversed. That is, the change from rip- 
ples dunes can occur different slope (depth and temperature held con- 
stant) than the change from dunes ripples. This gradual change and/or hys- 
teresis effect results transition from one bed form another. This tran- 
sition major importance when the bed form changes from dunes plane 
bed vice versa. Changes bed form occur many natural rivers re- 
sult achange discharge. This change will occur the slope the energy 
grade line close the critical slope that change inQcan cause change 
bed form. For example, the change from dunes low stage plane bed 
high stage. When this change bed form occurs there will break the 
discharge rating curve for that stream. Because there hysteresis effect 
the transition from dunesto plane bed plane bed dunes, the stage where 
the break the rating curve occurs depends whether the stage rising 
falling and also the rate change discharge. 

important note that the velocity the flow alluvial channel 
not zero the bed. The velocity the water through the bed, even though 
slow, important, especially when considering seepage force, which force 
form drag exerted the liquid the sand grains within the bed. The 
magnitude the seepage force per unit volume bed material equal 
the hydraulic gradient the flow and acts the direction flow. When the 
magnitude the seepage force equal magnitude and opposite direction 
the submerged weight unit volume bed material (inflow condition), 
the effective weight the bed material tends toward zero. Conversely, the 
seepage force and submerged weight bed material act the same direction 
(outflow condition), the effective weight the bed material increased. The 
direction flowin the bed canvary from vertically upward vertically down- 
ward depending bed configuration, direction flow the bed, water table 
condition and regime flow. The velocity flow the bed largest with 
the dune-bed form roughness because the porosity the material rela- 
tively large due the relatively loose packing the bed material resulting 
from the dune activity, and smallest the rapid flow regime where the 
bed material more closely packed. the result the floor, the 
velocity the flow through the bed was probably smaller than most natural 
streams with deep sand beds; consequently, the seepage force was smaller. 
natural stream the velocity has vertical components because inflow to, 
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outflow from the bed. The rate inflow outflow, which depends the po- 
sition the water table the surrounding area, can cause large seepage 
forces. These forces may have considerable local effect the form bed 
roughness and transport sediment. 

Tranquil Flow Regime. 

Plane Bed Without Movement.—With the plane bed, and movement bed 
material, the bed was soft and easily disturbed. The word plane used 
emphasize that the bed was not hydraulically smooth. For hydraulically 
smooth rigid bed (the ratio median diameter bed material the 
thickness the laminar sublayer) must less than 0.25. The minimum d/o' 
for the flume runs was 0.40. Obviously, possible obtain hydraulically 
smooth boundary without bed movement decreasing the slope the depth. 
smooth bed plane bed without sediment movement itself has 
little practical significance. The important factor consider the shear 
which movement bed material will begin the slope depth both are 
increased. was interesting note that the shear which movement began 
with this bed material was not the shear which the hydraulically smooth 
boundary became hydraulically rough. That is, the thickness the laminar 
sublayer decreased beyond the shear corresponding the transition before 
movement began. must remembered that the bed material was non- 
uniform and computations were based median diameter the material. 

Movement bed material began with slope 0.00015 when the depth was 
increased (by increasing from 0.61 0.71 the value the be- 
ginning movement the bed material was between 0.48 and 0.53. The 
Manning for bed material movement was approximately 0.015, and C/Vg 
varied between 14.0 and 15.0. 

When approached 0.53 movement the material began and ripples 
started form. The ripples first formed the upstream end the flume 
and any small depression ridge the bed, and then continued develop 
downstream direction from these points minor disturbance until the bed 
was covered. ripples formed the bed, the slope and depth the flume 
increased from 0.00015 and 0.71 0.00023 and 0.80 respectively. The Man- 
ning increased from 0.022 and decreased from 15.0 10.3. 

Ripples started form when movement began. Other investigators have 
reported movement with plane bed without the formation ripples. Whether 
this movement without ripples physical fact the result the experi- 
mental procedure and equipment unknown. Possibly, movement without 
ripples results sorted bed material used. The material used these 
experiments was unsorted. The width-depth ratio also may important. 

Ripples.—The form the ripples illustrated Fig. and photograph 
typical ripple formation shown Fig. For all ripple runs the height 
the ripples was less than 0.10 and their longitudinal spacing was less than 
2.0 ft. Their length-height ratio L/h varied from 20. The ratio depth 
flow ripple height ranged from 24, Manning ranged from 0.019 
0.027, and ranged from 7.8 12.4. The upper limit n-values, and the 
lower limit general, were associated with the shallow depths. 

There was little suspended sediment movement with the ripple bed 
form. The water was sufficiently clear that the visible all times. 
The total sediment load ranged from ppm 101 ppm and the sediment moved 
more less continuous contact with the bed rolling over the crest 
the ripples and coming rest their downstream face. This sediment did 
not move again until became exposed the upstream part the ripple 
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FIG, 3.—VIEW RIPPLE CONFIGURATION, LOOKING UPSTREAM. 


FIG, DUNE CONFIGURATION, LOOKING UPSTREAM. 
WATER-SURFACE SLOPE 0.00289. 
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the ripple migrated downstream. Thus, the movement the sediment parti- 
cles took place steps. The length the step and time interval between 
steps depended the length the ripple and ripple velocity. 

With ripples, the bed was soft but not soft dunes. The water sur- 
face was very smooth with little visible difference from the water surface for 
the plane bed without sediment movement. The separation zone downstream 
from the crest the ripple caused little jet impingement the downstream 
ripple. The reverse flow the trough moved only the finest particles, and the 
turbulence was small the interface between the main current and the sepa- 
ration zone. 

When the slope and/or increased beyond certain limit the rip- 
ples were modified dunes. Because the difference relative roughness 
and resistance flow necessary treat the two forms bed roughness 
separately. Compare the ripples shown Fig. with the dunes Fig. 

The bed changed from the ripple pattern dune pattern when the slope 
and/or depth were changed such that was approximately 1.0, and the 
Froude number was approximately 0.28. The change was abrupt, and the dunes 
established themselves over the full length the flume ina few hours. 
change from dunes ripples constant), was necessary decrease the 
slope and increase depth from the values which ripples changed dunes. 
With the reduced transport capacity caused the change slope and depth. 
considerable time was required for the flow convert completely dune con- 
figuration ripples. However, ripples superposed the backs 
the dunes short time. 

Dunes.—The form dunes illustrated Fig. and photograph typ- 
ical dunes shown Fig. Considering all dune runs, the dunes ranged 
average height from 0.15 0.52 and length from 4.0 7.5 ft. The 
length height ratio L/h varied from 28. The ratio depth flow 
dune height D/h ranged from Froude numbers varied from 0.38 0.60 
fixing the upper limit the unmodified dune bed form. Manning’s varied 
from 0.018 0.035, although observing the bed, much larger values Man- 
ning would have been expected (Fig. 4). The reason for small n-values 
that with large flume, with natural stream, the dunes form pattern 
that allows part the flow sideslip bypass the dune. Sideslippage 
meandering the flow around the dunes was very obvious when smaller dis- 
charge was run over dune bed formed larger discharge. 

Dunes move downstream result sediment toppling over the crest 
and accumulating the downstream face the dune. The larger the ampli- 
tude the dune, the more sediment stored the dune and the smaller the 
dune velocity for constant transport rate. Considering particular run, 
dunes with small amplitude and with their inherent high velocity overtake the 
larger dunes. This results much larger dune and decreases the dune ve- 
locity. Velocity the dunes varied from 0.02 fpm 0.70 fpm. The large 
dune velocities were associated with relatively steep slopes and shallow depths. 
Conversely, the small dune velocities were associated with relatively flat slopes 
and deep depths. knowing the average amplitude and velocity ripples and 
dunes, the bed load transport can 

The (L/h)-ratios for dunes and ripples were similar. The ratios varied 
from for ripples and for dunes. This indicates that there 
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correlation between dune height and dune length for given sediment. That 
is, the amplitude the dune ripple increases with change slope 
depth, the length also increases. The variation the (L/h)-ratios may the 
result the randomness the dune configuration—which makes difficult 
obtain true average length and amplitude. Some the variation L/h may 
result from including the separation zone measured lengths the dunes. 
Whereas the true length the dune may only extend from the end the sepa- 
ration zone the next crest downstream, the measured length, used, was 
from crest crest. The magnitude the separation zone increased slope 
increased until the transition regionwas reachedwhere the dunes began 
wash out. 

The ratio depth flow dune height (relative roughness) D/h was 
different for ripples than for dunes. Ripples have larger (D/h)-ratio than 
dunes. This parameter may help explain why the n-values are smaller for rip- 
ples than for dunes. Also with the ripple and dune-bed forms, the largest 
values occurred with the smallest (D/h)-values. 

appears from the flume data which for essentially two-dimensional 
flow that maximum dune height limited the depth flow. That is, dunes 
will only grow until their crests are within certain relative distance the 
water surface. phenomenon does three-dimensional 
flow. Also, with given slope and continuous increase depth the dune 
height may reach some maximum value for that size sediment and water 
temperature, and increasing the depth beyond this limiting depth will not in- 
crease the dune height. this true then the shear which dune height be- 
comes independent depth must determined predict accurately the re- 
sistance flow alluvial channel. This shear may function the 
width dune height becomes independent depth, the n-values 
may decrease with increasing stage. That is, h/D the relative roughness will 
decrease. 

With dunes, the bed was very soft and fluid and there was considerable 
segregation the bed material. The surface the upstream slopes and crests 
the dunes, which had negative slope approximately 4°, consisted 
the fine and coarse material moving downstream. The slopes the upstream 
faces were practically the same those the upstream faces the ripples. 
Ripples formed the upstream faces the dunes the smaller slopes. 
large slopes, the ripples the dunes were not apparent. The downstream face, 
which had positive slope between 31° and 39°, largely consisted the coarse 
fraction the bed material. Thecoarse fraction was the material which could 
not swept into suspension the crest and avalanched down the fore slope. 
was deposition avalanching which caused the dune bed soft and 
fluid. The material was deposited the force gravity and had relatively 
large voids. had been deposited the dynamic force the fluid, would 
have been packed into place, with relatively small voids ratio. Bagnold 
observed? and explained the same phenomena with wind-blown dunes. The 
troughs the duries contained horizontal layer 1/8 in. in. thick the 
fine fraction the material. Part the fines the trough came from the 
sediment which was swept over the crest and part came from the lower part 
the upstream face the adjacent downstream dune where the main current 
overrode the separation zone and impinged the bed (Fig. 2). the water- 
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surface slope was increased, with the dune form bed roughness, the toe 
each dune was blasted and eaten away the high velocity water that impinged 
it. The erosion provided source additional fine material which was de- 
posited the trough the reverse flow the separation zone. 

The slopes the downstream faces the dunes, which varied from 31° 
39° and averaged 34°, were clearly defined lines formed the dark-colored 
particles the dune moved downstream (Fig. 2). The angle repose air 
for dry non-cohesive sand (the bed material the flume was 95% finer 
than mm) varies from 29° 32° depending angularity the 
probable that the angle the forward face the dune was greater than 
the angle repose dry bed material air result pressure the 
fore-plane the dune the separation zone and the reverse flow the trough 
moving the face the dune. This was verified fact that the ve- 
locity reduced zero the flume, the downstream face the dune 
slumped smaller angle inclination. 

With the dune-bed form roughness, the water surface was uneven and tur- 
bulent and the water surface was lower over the troughs than over the crests 
the dunes. This was the result the acceleration and deceleration the 
flow contracted over the crest and expanded over the trough. This was 
illustrated for rigid boundaries Rouse.9 The roughness (turbu- 
lence) the water surface increased the bed roughness increased. The 
roughest water surface observed with dunes occured when there were large 
water surface boils. The water surface the boils was approximately 0.1 
higher than the surrounding water surface. However, the roughness the 
water surface was not great enough affect the accuracy the water surface 
slope measurements. The turbulence created between the main 
flow and the flow the separation zone dissipated considerable energy and 
normally the disturbance caused the interference was visible the water 
surface downstream. With the dune-bed form, the suspended sediment concen- 
tration, the intensity the turbulence, the relative roughness, and the velocity 
the reverse flow the trough increased the slope increased, Along with 
the dunes, potholes formed which had depth equal the height the dunes. 
These potholes, well the dunes, caused boils the water surface that 
were evident downstream from The potholes and boils moved down- 
stream front the dunes the velocity the dunes. Normally from ten 
twenty potholes and boils were evident the full length the flume with this 
type bed roughness. The pothole appeared develop result the in- 
creased strength the secondary circulation the separation zone. The 
combination the velocity the direction flow and the secondary circula- 
tion caused circulation similar whirlwind. This rotating motion 
scoured out additional material point the separation zone and produced 
pothole. The existence these potholes was indicated the formation 
very strong boils which visibly transported considerable suspended sediment 
upward the water surface. 

The segregation material the formation and movement dunes means 
thata large number bed material samples must obtained and analyzed be- 
fore the size distribution can established with confidence. With the flume 
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material, sixty samples were needed before another sample would not alter the 
median diameter significant amount. Segregation should have some value, 
though, because the field (Fr 0.6), might possible determine the 
characteristics the bed roughness core sampling the bed. 

Transition, Dunes Standing Waves.—The change from dunes standing 
waves was complex, and the form the bed roughness was erratic. The tran- 
sition occurred when the depth and/or slope were changed that and 
<1.0. 

With dunes, the maximum Froude number was 0.60. Runs with Froude 
number between 0.6 and 1.0 displayed multiple roughness, that bed form 
between dunes and standing waves which consisted washed-out dunes and sand 
bars. The washed-out dunes and sand length and had 
<h)-ratios between and 50. 

logical that the upper limit the Froude number, with the dune-bed 
form, considerably less than 1.0. The change roughness, and consequent- 
ly, the change resistance flow and the dissipation energy were large 
when the bed changed from the undulating form dunes, whereas the change 
energy from change the Froude number was small when =1. Thus, 
the velocity and depth were changed considerably when the bed form changed 
and small Froude numbers resulted. Conversely, the change from dunes 
rapid flow resulted smaller loss energy because the reduced rough- 
ness, which caused larger velocities, small depths, and larger Froude number. 

There was hysteresis effect the change bed roughness from dunes 
standing waves and back dunes. The value slope and/or depth for the 
change depended the bed configuration prior the change. the bed was 
covered with dunes, slope 0.0035 was perhaps required before the dune 
configuration was destroyed. the bed was standing waves, slope about 
0.0025 was required before dunes were formed. This hysteresis effect may 
have resulted because the energy changed the bed form changed from dunes 
standing waves standing waves dunes; that is, potential energy changed 
kinetic energy when the form changed from dunes standing waves, and 
kinetic energy changed potential energy when the form roughness changed 
from standing waves dunes. The changes energy resulted from the large 
changes roughness associated with the changes Dunes had 
values which ranged from 0.018 0.035, whereas the standing waves had 
values which ranged from 0.015. 

The fact typical dune pattern did not occur with this bed material, and 
Froude number equal greater than 0.60 important, especially the 
phenomenon for other bed materials. formfor 0.6 <1.0 
was washed out dunes with rather large (L/h)-values. The Manning n-values 
were 0.018 and 0.021 respectively; these n-values are smaller than n-values 
for dunes, but are larger than rapid flow n-values. With dunes, seems that 
the magnitude Manning related both spacing and amplitude the sand 
waves. 

Some runs had distinctly different types flow occuring side This 
was probably because discharges which were too small with respect slope 
were used which produced three-dimensional flow. This same phenomena was 
also observed the antidune range when the discharge was decreased below 
given value. With three-dimensional flow, the n-values ranged from 
0.023, whereas with standing waves, the n-values were approximately 0.012. 

Three-dimensional flow probably resulted because, the experiments, 
discharge and slope were controlled and with the steeper slopes the small 
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discharges set conditions favorable for three-dimensional flow. Other ex- 
perimenters may have experienced this flow condition because the limited 
discharge capacity their flume system which caused large width depth 
ratios develop when slope was increased. With the dune bed form, the 
smallest discharge used was sufficient insure two-dimensional flow. Pre- 
sumably, however, further decreasing the discharge, three-dimensional 
flow could develop. 

Rapid Flow Regime.—With rapid flow, 1.0, three forms sand bed 
and water surfaces were observed: plane, standing waves, and antidunes. 

Plane Bed and Plane Water Surface.—A completely plane bed and plane 
water surface, for the full length the flume was only produced once the 
total sequence runs. was anticipated from existing literature that this 
would the type bed configuration following dunes. However, for this sand, 
such was not the case. Based recent tests which fine sand was used, 
apparent that the development plane bed intimately related size 
bed material. With fine sand bed material plane bed runs are common 
phenomenon, and they develop the tranquil flow regime. 

Standing Waves.—The water surface consisted symmetrical standing waves 
small amplitude. The standing waves formed and gradually disappeared, and 
unlike antidunes they had tendency break migrate upstream. The 


forms bed roughness observed with standing waves, the order increas- 
ing slope, were: 


(a) diagonal dune pattern cross-laced like shoestring. 

(b) plane bed. 

(c) symetrical undulating sand wave similar form those observed 
the antidune regime. 


The standing waves formed when Froude number was approximately 1.0; 
the undulating sand wave bed developed when 1.2 The Manning 
for standing waves was relatively small, ranging from 0.012 The 
water surface waves were roughly 1.5 times high the corresponding sand 
waves. Concurrent profiles the sand bed waves and the water surface waves 
are illustrated Fig. 

the Froude number, which was computed the basis 
average velocity and average depth, was greater than 1.3,and was greater 
than 3.0, antidunes formed. 

Antidunes are defined train symmetrical sand waves which are 
phase with corresponding train symmetrical water surface waves. Both 
trains waves move upstream, grow height, then break (Figs. and 
causing cyclical fluctuation the water surface waves and the sand waves. 
The waves built from plane bed with plane water surface. They grew 
and moved upstream until one two the waves became unstable and broke. 
Normally, when one wave broke other waves broke immediately thereafter for 
distance waves upstream and waves downstream. That is, 
usually broke within short time interval. The chain reaction that followed 
the breaking the first wave was apparently triggered the action the 
first wave which broke. 

While making concurrent measurements the bed and water surface the 
troughs and crests antidunes, was observed that the antidunes became un- 
stable and broke whenever the water surface the trough the wave train 
was approximately the same elevation the crest the downstream bed 
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FIG, 5.—ANTIDUNE WAVE POINT BREAKING 
LOOKING UPSTREAM. 


FIG, 6.—ANTIDUNE WAVE AFTER BREAKING, 
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FIG, 7.—THE ACTUAL WATER SURFACE AND BED 
CONFIGURATION FOR RUNS AND 


| 
1.20 
| | | | 
ce 
| | 
| | | | | | 
| | | | 


wave (Fig. 7). The measurements also indicated that the water waves were 
1.7 times larger amplitude than the corresponding sand waves. The total 
height the water surface waves from trough crest was from 
times the average depth. When the waves started break their heights were 
about twice the average depth. The breaking the waves was very spectacular. 
There was considerable turbulence, dissipation energy, and mixing the 
flow. When waves broke they sounded like the surf the ocean and there 
probably some resemblance the breaking ocean wave. 

When the antidunes were building the bed was very firm and there was 
separation between the flow and bed. However, when the antidune broke, 
the crests the sand waves became soft and fluid and seemed almost ex- 
plode. When the waves broke, the flow was very turbulent, there was sepa- 
ration, and considerable sediment was thrown into suspension. Total loadcon- 
centrations were very large and ranged from 4,240 ppm 15,000 ppm. ap- 
peared that the upstream movement the antidunes resulted from scour 
the downstream side the sand wave and deposition the upstream side 
the wave—consequently, the wave moved upstream. 

Except for one run, which was considered three-dimensional, there was 
only one train waves the cross section. The train waves was located 
off center about the left 1/3-point the cross section. subsequent runs, 
however, two trains waves symmetrically located the flume were observed 
several times. The wave trains were not continuous throughout the entire 
length the flume. That is, train antidunes would either build up, break 
they would build and break two three discontinuous lengths the 
flume. The building and breaking the different reaches the flume could 

The period time took for one antidune wave build and break, and 
the number antidune waves which built one time, varied with the slope 
and the discharge. The antidune run with the flattest slope had only one train 
antidune waves which would build and break giventime. This 
train waves which was about long, would perhaps build the lower 
section, then the middle upper sections, and again some other position 
random Normally, about two complete cycles occurred every 
hour—a cycle includes starting with plane surface, the building waves, 
the breaking the waves, and back plane surface. For other runs with 
the steepest slopes there were three four reaches the flume which 
the train waves built and broke both phase and out phase, almost 
continuously. 

When the antidunes broke, considerable amount water was stored the 
flume. Based observations the glass-walled section the flume, when 
antidune broke the water the crest moved upstream, the water close 
the bed almost ceased move until the wave vanished, and shortly thereafter 
normal flow was restored. With the large slopes, when three four series 
antidune waves the train would build and break simultaneously, water was 
stored the flume such extent that the pumps would surge. the ex- 
treme case, antidune waves would build and break phase with the surging 
the pumps. This caused the discharge fluctuate and the antidune activity 
increase. Consequently, throughout the full length the flume much 
water was stored the flume that the tail-water level the tailbox dropped 
until the pumps lost their prime. the surging the pumps, the level 
the water the sump was raised and excess water was continuously added 
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the sump replace that stored the flume. When the antidunes were not 
breaking the surplus water was discharged through overflow the top 
the sump. 

The storage water caused the breaking the antidune waves probably 
accounts for the surging discharge that sometimes observed alluvial 
streams with steep slopes. That is, with antidunes water stored and re- 
leased the upper reaches the stream random haphazard manner, but 
the flow travels downstream the fluctuation discharge caused the anti- 
dunes break more systematic pattern until surges develop which cause 
the antidunes form and break regular time intervals. Some streams that 
surge are Muddy Creek Wyoming and Medano Creek the Great Sane Dune 
National Monument, Colo. 

Antidune flow using this bed material, hydraulically speaking, was very ef- 
ficient. The Manning ranged from 0.010 0.013. Even though considerable 
energy was dissipated when the antidunes broke, the period time and length 
flume over which the antidunes broke was small comparison with the to- 
tal time and total flume length. This explains the small n-values. was ob- 
served slope was increased antidune activity increased, the dissipation 
energy increased, the discharge coefficient C/Vg decreased, and Manning 
increased. 

Surface velocities, obtained timing floats along center line train 
antidune waves (when they were not breaking) and along nearby parallel line 
where the water surface was smooth, showed that the surface velocity was 
greater the train waves. The velocity flowin the troughof the antidunes 
was considerably greater thanthat inthecrest. possible thatthe antidunes 
break when the velocity decreases and depth increases over the crest because 
the wave height increases until Froude number less than 1.0 results. 


ANALYSIS DATA 


should mind that the relations presented are based only one 
group flume data, and that the size and gradation bed material remained 
constant except altered miscellaneous sorting which mayor may not have 
been significant. That is, were constant and varied only because 
temperature changes. 

Variation Velocity With Depth.—The variation velocity with depth and 
interrelationships these two variables with the six major bed rough- 
ness are given Fig. Lines equal Froude number have been superposed 
the figure illustrate the wide scope the data. Velocities range from 
those for movement bed material those far excess the critical 
velocity, from 0.65 fps 6.2 fps. Note that the transition roughness includes 
those runs which bed roughness was somewhere between dunes and the rough- 
nesses associated with rapid flow (washed out dunes Data rep- 
resenting runs this type plot just below the critical velocity curve. The 
lines dividing the plot into forms bed roughness are based the observed 
roughness and the Froude number. 

Variation Resistance Flow.—Considering the complete range condi- 
tions investigated the Manning ranged from 0.008 0.035. These n-values 
have been corrected for side wall effect accordance with the procedure pre- 
sented The magnitude increased the slope in- 


“River Channel Roughness,” byH.A. Einstein and Barbarossa, 
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creased and the depth remained constant, until the regime flow shifted from 
the tranquil regime the rapid regime; with rapid flow n-values decreased 
magnitude quite suddenly approximately 0.012. 

the tranquil flow regime the Manning was 0.015 for the plane bed and 
bed material movement, ranged from 0.019 0.027 for ripples, and from 
0.020 0.035 for dunes. the transition from tranquil flow rapid flow the 
and slope from large resistance, which was slightly smaller than the resist- 
ance for dunes small resistance, which was slightly larger than the re- 
sistance for rapid flow. With dunes, the magnitude the resistance flow 
dependent the spacing and amplitude the dunes. 

the rapid flow regime, Manning ranged from 0.0078 0.015. The fact 
that resistance flow much smaller than for dunes can attributed the 
change bed roughness. Dunes have large separation zone with large form 
drag, whereas the symmetrical sand waves had little separation zone, 
only the form drag the particles. Possibly some the decrease Manning 
can attributed the movement large quantities sediment bed load 
and not entirely the elimination the dunes. the rapid flow regime insuf- 
ficient turbulence was created hold large quantities sediment suspen- 
sion. Consequently, the sediment load which was carried suspension with 
dunes was carried near thebed. This large concentration near the bed, turn, 
markedly changed the properties the fluid-sediment mixture. The fluid was 
not homogeneous and sort stratified flow resulted which inhibited mixing 
and reduced the effective bed roughness extremely small values. The flow 
had some similarity plug flow pipes. During the single run for 
0.0078, the bed was plane and the resistance caused bythe standing waves was 
not effective. Under these conditions the plug-flow phenomenon was even more 

Variation Tractive Force and Total Sediment Load.— Fig. was obtained 
plotting the tractive which has been calculated 


against the total sediment load concentration. Although there appreciable 
scatter, definitely exists. The magnitude total load increases 
the tractive force increases. Fig. each run has been symbolized 
accordance with its form bed roughness. There discontinuity between 
the tranquil-flow data and the rapid-flow data because the relatively large 
tractive force associated with runs with large dunes. 

The deviation the runs with dunes from the general trend resulted from 
increase due the extreme size the dunes, withaut cor- 
responding increase inconcentration significant that sudden 
decrease occurred about 2,000 ppm the bed changed from 
plane bed standing waves. This was caused bya reduction rough- 
ness, which also decreased the ability the flow transport sediment. Con- 
sequently, was cut half from tun run 36, beyond which point 
the transition regime steadily increased. discontinuity occurred the 
change from ripples dunes from standing waves antidunes. There was, 
however, considerable scatter the rapid flow regime, which shows that other 
variables are needed define the phenomenon more completely. 

and the form bed roughness was investigated for both regimes flow 
(Fig. 10). importance note the precision with which the relation 
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describes and systematically groups the various forms bed roughness. The 
relation can modified include the effect size and gradation bed ma- 
terial. also significant that either both and can utilized 
the third variable Fig. 10. 

Additional study and adjustment parameters showed that for the bed ma- 
terial used, C/Vg can computed the product constant and the square 


the Froude number added the parameter This 
Fig. 11, from which 


The form bed roughness plots very well third variable except for 
minor intermingling ripple and dune runs the arbitrarily selected di- 
viding line. also possible utilize total sediment load third vari- 
able with useful accuracy. 

Based similar approach, but using the additive term in- 
stead figures similar Fig. which are only slightly less accurate 
can obtained. fact figure this type that involves could used 
estimate qualitatively total sediment load the rapid flow regime. 


CONCLUSIONS 


Resistance flow alluvial channels varies between wide limits and 
extremely comples becausethe form the bed roughness function fluid 
properties, flow and sediment characteristics, and channel geometry. The 


major bed roughness that were determined from experiments and that 
are discussed this report are: 


Tranquil flow regime, <1. 
Plane bed without bed material movement. 
Ripples. 
Dunes with ripples superposed. 
Dunes. 
Transition from dunes plane bed. 
Rapid flow regime, 
Plane bed. 
Symmetrical standing waves. 
Antidunes. 

the tranquil flow regime the resistance measured the Chezy 
coefficient discharge was 14.0 for the plane bed without movement, 
varied from 7.8 12.4 for ripples, and from 7.4 12.8 for dunes. The large 
values for the ripples dunes were associated with the deepest depths, 
and hence, the smallest relative roughness. 

the rapid flow regime varied from 13.9 largest 
value occurred with plane bed and plane water surface. Standing waves, 
general, had larger values than antidunes. This resulted from the dis- 
sipation energy the breaking waves. The value could reduced 
further for antidune flow but for the fact that the period time and length 
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flume occupied breaking waves was small comparison with the total time 
and length flume. 

was observed that with the bed material and the flume used the experi- 
ments, dunes started wash out and resistance flow decreased when 0.6. 
the transition toplane bed varied from 10.0 14.1. 
magnitude increased the percentage the bed which wascovered 
with dunes decreased. was ahysteresis when the dune bed changed 
standing waves which depended the form the bed prior the change. 
That is, with discharge constant, larger slope was required change from 
dune bed standing wave bed than was required change from standing 
wave bed dune bed. 

For certain size gradations bed material, the slope the energy grade 
line close critical slope, change stage causes dunes change 
plane bed standing waves vice versa. This phenomenon occurs many 
natural streams and produces discontinuity the stage discharge relation. 
However, because the hysteresis that associated with such change 
bed form, the stage which the discontinuity develops depends whether the 
stage rising falling and the rate change discharge with time. 


ADDITIONAL REFERENCES 


Simons and Richardson, Proceedings, ASCE, Vol. 84, No. HY1, Feb- 
ruary, 1958. 


Study Sediment Transport Alluvial Channels,” Barton and 


Lin, Dept. Civ. Engrg., Colo. State Univ., Fort Collins, Colo., No. 
1955, 43. 


“Mechanics Streams with Movable Beds Fine Sand,” Norman 
Brooks, Transactions, ASCE, Vol. 123, 1958, pp. 526-549. 


Kalinske and Hsia, Studies Engrg., Bulletin No. 29, State Univ. 
Towa, 


“Some Factors Stability Canals Coarse Gran- 
ular Materials,” Lane and Carlson, Internatl. Hydr. Conven- 
tion Proceedings, Minneapolis, Minn., 1953, pp. 37-48. 


“Mechanics Sediment-Ripple Formation,” Liu, Proceedings, 
ASCE, Vol. 83, No. Hy2, April, 1957. 


LIST SYMBOLS 


Symbols Description 


dimensionless form which 
equivalent 


Drag coefficient for the particle 


7 
: 
4 
4 
_ 
oo 
a 


ALLUVIAL CHANNELS 


Concentration total load 


Concentration fine material 
Median fall diameter bed material 
Average depth flow 

Froude number 

Average height bed roughness 
Any constant 

Average spacing bed roughness 


Mannings coefficient roughness 


Mannings coefficient roughness 
(Einstein) 


Discharge water-sediment mixture 
Reynolds number 

Shape factor the channel cross section 
Shape factor for the reach the stream 


Slope energy gradient equal water 


surface slope steady, uniform flow 
Time sec 


Average velocity based continuity 
principal 


Shear velocity which ft/sec 
Fall velocity sediment particles ft/sec 
Specific weight water 
Specific weight sediment 


Difference between specific weights 
air and water lbs /ft 


Difference between specific weights 
sediment and water lbs/ft 


Thickness laminar sublayer 


Kinematic viscosity /sec 


Dynamic viscosity 
Mass density water slug 
Mass density sediment 


Relative standard deviation the size 
distribution the sediment 


Tractive shear force developed 
the bed, 
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HYDRAULIC ANALYSIS SURGE TANKS DIGITAL 


Discussion Escande 
Closure Nicholas Barbarossa 


ESCANDE.!—Mr. Barbarossa’s paper constitutes important contri- 
the study results obtainable, electronic computation, the field 
certain phenomena involving oscillations water surge tanks. 

The examples studied relate system two surge tanks arranged 
series the same penstock. The treatment the problem takes into account 
large number factors which general are neglected classic theory, 
but does not introduce all such factors. Thus, the compressibility the water 
and the elasticity the walls the pressure conduit are neglected, and 
consequence, the corresponding waves overpressure. 

the case simple operation, such closing opening, the results 
obtained can considered very good. The same cannot said the study 
stability, for which the results are extremely sensitive the influence 
factors generally considered secondary, and which may lead results di- 
vergent that they cannot considered valid. 

Mr. Barbarossa concludes, from his studies, that programming for elec- 
tronic computations should carried greater degree precision than 
has yet been attained, and recognizes that the problem has not been solved 
this time. also seems advisable him very judicious accepting re- 
sults drawn from normal classic theory, because the simplifications neces- 
sarily introduced into their derivations. 

The writer has had occasion write paper2 analysis surge-tank 
phenomena the method finite differences, and quite recently have used 
electronic computers for the study certain problems relative the oscil- 
lations set successive periodic operations opening and closing. This 
study not yet completed, and thus cannot draw any formal conclusions. 
However, have noted the extreme sensitivity electronic computations and 
are not surprised the divergent results obtained the study stability. 
definition such phenomenon corresponds limiting case between two do- 
mains, and even minimum influence suffices orient toward one the 
other the two. Furthermore, hardly possible programming study 
stability take into account all the hydraulic phenomena, since the com- 
plexities would make the computation too involved for solution. 

The writer inagreement with Mr. Barbarossa onthe invaluable assistance 
which electronic computers lend the solution surge-tank problems. For 
sensitive problems such those stability, felt that Mr. Barbarossa’s 


April, 1959, Nicholas Barbarossa. 
prof., Faculte d’Electrotechnique d’Hydraulique. 


“Methodes nauvelles pour calcul des chambres d’equilibre,” Escande, 
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conclusions, concerning the use electronic computers and the possibilities 
practical accomplishment complete programming, are too speculative. 
the contrary, one should not too critical the classic methods for study 
stability. they are used with discernment, they may yield results suffi- 


ciently accurate serve practical needs, without introducing exaggerated 
margin safety. 


note the number and caliber the engineers who contributed discussions 
the paper. The author takes this opportunity express his deep appreciation 
all who have taken the time document their thoughts and suggestions 
the subject for the benefit designers and investigators. 

Since detailed physical knowledge the phenomena under discussion 
understandably inadequate, there can “exact solutions” computer. 
has already been pointed out, the computer merely servant with 
prodigious capacity for menial work comparatively small cost. enlight- 
ened exploitation this capacity, engineers may broaden their understanding 
problem phenomena from analytical standpoint but not, basically, from 
physical point view. Loosely stated, one may say that analytical research 
can done computing centers, but physical research has done the 
same old places—in the laboratories and field. seems rather 
common find that each type research stimulates the other. Thus, hap- 
pens, one example, that computer studies surge tank analysis described 
the paper, based essentially upon conventional methods long standing, have 
induced new lookat the whole physical matter power plant transients. The 
new look has prompted the same questions regarding physical relationships 
that must have been asked many times the past, when mechanical computing 
centers were made only slide rules and desk calculators. Some these 
questions have been raised discussions the paper. Before attempting 
deal with these questions and other points that have been raised, brief de- 
scription pertinent studies recently undertaken the United States Corps 
Engineers Dept.of the Army will given demonstrate what the new look 
that organization has produced that interest and promise engineers. 

The proposal for comprehensive study digital computer has been accept- 
basis development IBM 704 computational programs. These pro- 
grams are now being developed the Massachusetts Institute Technology 
contract with the Missouri River Divisionof the Corps Engi- 
neers. Computational programs are organized such manner 
preserve the general nature solutions the highest practicable degree and 
are compartmented way facilitate modifications that would 
necessary for application these programs power plant systems dif- 
ferent but typical geometrical arrangement. Machine programs are de- 
veloped for the and Oahe power plant systems the Missouri River 
and, possibly, for one other undesignated power plant system. two-year 
study, terminating sometime the summer 1961, has been set up. re- 
plannedthat the results this work will made avail- 
able the engineering profession through publication. 


Chief, Hydr. Sect., U.S. Army Engr. Div., Missouri River, Corps 
Neb. 


Located about river miles upstream Bismarck, 
New Approach Hydro Analysis,” Engineering News-Record, December 1959. 
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Since computational programs power plant transients should based 
actual operational behavior observed and measured the field, whenever 
conditions permit, detailed field tests the Garrison and Oahe power plant 
systems have beenscheduled the Corps Engineers order that computer 
results may evaluated sound basis practicable and order 
determine what extent operation corresponds with design. The planning 
these tests being done close coordination between and several of- 
fices the Corps Engineers, including the Waterways Experiment Station, 
Garrison and Omaha Districts, and the Missouri River Division. The first 
series field tests will made during the summer 
1960. For each test run, planned that both magnetic-tape andoscillograph 
equipment will utilized record electrical synchronization the instan- 
taneous values many pertinent variables (pressure and total head com- 
paratively large number points between intake and tailrace, wicket gate 
opening, turbine speed, power output, and forth). The response the many 
pressure transducers that will utilized and other instrumentation will 
have approximately 100 cycles per second order keep satisfac- 
torily with rapid pressure changes during the water hammer phase test run. 
Taken together, the analytical computer studies and directly related operation- 
tests should yield data and design information unprecedented nature. 
Mr. Jaeger, Messrs. Danel and Ransford, and Messrs. Herak and Rodriguez 
have all brought the question tee losses. The writer concurs that syste- 
matic experiments are urgently needed with regard tee losses under un- 
steady flow conditions and that computer studies could 
quite fruitful. Dynamical scale models which surges are actually pro- 
duced the same frequencies those expected the prototypes are likely 
expensive. However, the results thus derived may well the best that 
can obtained for particular installations stage. The assurances 
gained would probably permit closer design the prototypes, thus achieving 
economies which may more than offset the expense the models. 

While there likely serious disagreement the great value 
well-planned model studies applicable cases, questions arise sometimes 
the value prototypetests for plants which appear operating 
factorily. these plants are operating satisfactorily today, because 
the excellent work that has been done the past. The debt owed the past 
canonly discharged making contributions the future. Continuing proto- 
type tests, together with applicable model and computer studies, are urgently 
needed advance the engineers’ over-all knowledge power plant transient 
phenomena that tomorrow’s power plants may designed more economi- 
cally and may operate even better than they today. 

The original purpose the IBM 650 program described the paper was 
develop solutions load rejection, load acceptance, and surge-tank stability 
problems machine essentially the same way that these problems have 
been solved many times the past hand. demonstrated, the computer 
handles the first two types problems satisfactorily but not the third type. 
Several suggestions have been offered how the IBM 650 computer pro- 
gram might improved handle the stability cases. Mr. Jaeger proposes 
the inclusion the governor equation and Mr. Swiecicki the introduction 
speed drop based upon certain assumptions. Despite simplifying assump- 
tions, Mr. Swiecicki’s quantitative analysis provides revealing assessment 
the effects speed drop upon the stability problem. The writer’s proposal 
for the comprehensive IBM 704 program now under development included the 
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governor equation and speed relationships well many others. However, 
does not appear feasible include these relationships the writer’s 650 
program now stands, for several reasons. First all, the 650 program 
already very large, and major additions thereto might overtax the capabili- 
ties the IBM 650 (with minimum peripheral equipment) but, possibly, not 
those 650 computer with magnetic tape units, indexing registers, etc. 
Secondly, the IBM 704 program when developed will provide much better ans- 
wers because based upon minimum number simplifications. Thirdly, 
the relative facility solution characteristics the existing 650 program, 
depending, does, upon knowing the gate opening function time, 
would lost inclusion two more differential equations general form: 


these equations P(G,H) the power output the turbine, function and 
(and efficiency) or, alternatively, and L(t) denotes external electrical 
lar speedof turbine and C’s are the mathematical constants into 
which the physical constants the governor system are lumped;and other no- 
tation has been previously defined. Eq. characteristic the governor 
systems the Garrison and Oahe plants. applicable cases the third order 
and, possibly, the second order terms might dropped from 30. 

These remarks are not intended discourage any further development 
machine programs this field for applications IBM 650 comparable 
computers. Depending upon the objectives and individual judgment, simplifi- 
cations may justified, particularly during preliminary stages design for 
particular installation when broad outlines are being established and refine- 
ments are not necessary. such cases least, machine solutions the total 
problem surge tank stability might based upon making consistent use 
the following: the lumping process suggested Messrs. Evangelistiand Poggi 
and implied Messrs. Danel and Ransford; the analysis Mr. Swiecicki 
(based turbine efficiency and other quadratic law 
for head losses Messrs. Oldenburger, Evangelisti, and Poggi; 
the method upon which the results Fig. are based (discussed later); 
and other simplifications. 

states that would “of immense benefit the profession” 
prepare complete details the writer’s program including flow charts and 
machine codes. The suggested documentation has been accomplished. Mr. 
Sutherland advises caution regarding the use turbine performance curves 
IBM 704 709 applications. The writer well aware differences that may 
occur between expected and actual turbine performance. However, the engi- 
neer confronted with this possibility any case, whether machine appli- 
cations are involved not. Further, obtaining 704 solutions with two three 
different sets turbine performance curves may not expensive 


option was included the program described the paper for determining head 
losses quadratic relationship, but the results given the paper are based upon use 
the versus tables. not clear why Mr. Oldenburger assumed that the head 
losses the paper were taken proportional the “absquare” the discharge. 


| 
fly 
j 
4 
: 
4 
| 
3 
2 


DISCUSSION 


109 


imagined, particularly when closer design may result, thus permitting possi- 
ble construction savings that may more than offset computer costs. Another 
real prospect that 704 hourly rates will drastically reduced time be- 
cause better equipment being available. fact, this has already happened 
some computing centers. Expected turbine performance curves are normal- 
based upon steady state conditions the model. vexing but probably un- 
answerable question is, how well these curves apply for unsteady flow con- 
ditions 

The suggestion Messrs. Danel and Ransford test the “error control 
method used conducting least one blank computer run artificially- 
constructed fourth order equation which the solution known” worthy 
serious consideration all who might wish use the basic methodof solution 
described the paper. This has not been done the his colleagues 
for the reason given the paper. 

While true that independent check for determining computational 
accuracy exists for non-linear effects, Messrs. Danel and Ransford point 
out, the writer believes that one can usually circumvent this difficulty hav- 
ing the computer print particular case, all pertinent 
many time steps seem necessary, and then spot checks can made manu- 
ally. Table the paper (as example) was presented order permit 
such checks. 

Following the suggestion and the method proposed Messrs. Herak and 
Rodriguez, additional IBM 650 computer runs basedon minor program changes 
have been made Test No. Fig. 11, see whether not the path the 
transient would converge the new steady point (Mg). Part the results 
these runs are presented Table and Fig. 13. 1.6 sec, computer 
results are directly comparable with those Manual shown Fig. the 
discussion Messrs. Herak and Rodriguez. should noted that probably 
Manual Study does not include the effect the acceleration head the 
length whereas the computer study does every case. 

Fig. 13(a) only three the seventeen loops traced the computer are 
depicted for the sake clarity. Table provides some the values 
and that were used plot the three loops well some the values 
other interdependent variables. following the loops shown Fig. 13(a) and 
others not shown, the computer necessarily uses very small time increments. 
Consequently, takes about computer time obtain results for 
about sec problem time. Because the relatively large amount com- 
puter time required obtain results for several cycles surge-tank water 
surface oscillations (having quarter cycle about sec) and because 
other working demands the 650 computer available for this problem, has 
been possible within the time permitted for preparation this discussion 
obtain computer results for only about sec problem time. Hence, mea- 
sure the damping the hydraulic system over several cycles, desired 
Messrs. Herak and Rodriguez, not available this time. Nevertheless, the 
results shown Fig. 13(a) are interest indicating what can done 
artificial manipulation gate openings computer and showing some con- 
vergence toward the new steady state condition (point Mg). Under the artificial 
conditions assumed, the horsepower delivered the turbine oscillates between 
minimum about 74,000 and maximum about 83,000 (78,000 being the 


Thus, believed that the discrepancy between the Manual curve and Test No. 
curve primarily due the difference values tee loss coefficients. 
The use has causes Test No. curve follow wider path. Compare curves for 
Test No. and Test No. Figure for example. 


iz 4 
4) 
— 
J 
5 
We ig 


“AS 


110 May, 1960 


desired new value). the plant this wide variation would not permitted, 


course, well-designed governor system interconnected with basically 
stable hydraulic system. 


With reference Fig. 13(b), will noted that small-period oscillations 


TABLE 10.--OAHE STABILITY STUDIES TEST NO, (CONTINUED) 


First Loop 


0.00 1537.40 1537.36 
0.25 37.40 37.36 
0.45 37.39 37.34 
0.75 37.39 37.31 
1.00 37.38 37.27 
37.38 37.26 
37.37 37.23 
1.50 37.36 
1.80 37.34 
2.20 37.32 
2.22 37.32 37.17 


37.31 37.18 


1536.48 1536.14 


36.48 36.13 
36.45 36.08 
36.43 36.06 
36.42 36.07 
36.41 36.07 
36.41 36.08 
36.39 36.09 
36.38 36,06 
36.37 36.04 


36.37 36.02 


1535.62 1535.27 


117.37 35.61 35.26 
120.91 
123.13 35,57 
123.90 35.56 
117.69 35.55 
113.75 35.53 35.20 
113.62 35.19 
115.26 
117.36 35.16 


118.88 35,51 35.15 


Obtained from continuity equation 


are imposed upon the falling water surfaces, WS; and the surge tanks, 
due cycling the gate opening, which changed approximately every 1.2 
sec. and would probably reach minimum values about sec and 
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would then rise, overshoot the new steady levels, and oscillate until dampened 
quiescence. indicates that the turbine gates are opening and indicates 
that they are closing. Although time lapse indicated between and 


Fig. 13(b), governor dead time 0.05 sec was actually used the com- 
puter program, suggested Messrs. Herak and Rodriguez. For each 

and interval, gate movement was assumed linear; that is, the gate 
opening was assumed change the rate 20% per second (five-second 
governor). 


Mr. Swain’s general remarks the use computers are worthy note, 
particularly engineers who have not yet had opportunity directly 
associated with computer applications engineering problems. These remarks 
substantiate the writer’s convictions and experience this field. With 
ence the “table look up” routine, should noted once again that specific 
turbine performance values (Q, values) during computer runs were 
obtained the author’s program double interpolation the input table 
(Table 3). This point emphasized order focus attention the fact that 
not always necessary enter mathematical relationships into machine 
programs. 

Mr. Escande states, one shouldnot too critical the classic methods 
for study stability whereas Mr. Jaeger congratulates the author “for 
his criticism the assumptions which the present surge stability theory 
based.” criticism such the classical methods was intended. The au- 
thor simply desired known limitations classical methods 
might overcome great extent through new approaches, made practicable 
only recently the advent modern computers. The writer cannot agree 

with that “the use electronic computers and the possibilities 

practical accomplishment complete programming are too speculative.” 

Granting that some speculationis involved, the case inany new endeavor, 
should not the new potentialities for substantially better understanding power 
plant transient comparatively inexpensive computer analy- 
ses fully explored? not necessary, more than ever, continue both 
types research, physical and utilizing the best tools available? 

Mr. Jaeger states that “further research surge stability definitely 
wanted.” suggests that comprehensive paper summarizing all these new 
efforts the up-to-date the problem and show the way 
further research.” His recent advocating further research, 


excellent review onthat subject and recommended toall who desire pursue 
the subject further. 


Review Surge-Tank Stability Criteria,” Paper No. 59-A-270, ASME, presented 
ASME Convention, Nov. 29-Dec. 1959. 
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RESISTANCE PROPERTIES SEDIMENT-LADEN 
Closure Vito Vanoni and George Nomicos 


VITO VANONI,! ASCE, and GEORGE authors are 
indebted the writers who, contributing discussions which re-examine 
some the ideas the paper, have added greatly its value. 

Messrs. Laursen and Chou have advanced the idea that the resistance 
flow offered grains the bed different when they are moving than when 
they are fixed and suggested that this may account for the lowering the re- 
sistance observed for sediment laden flows. Although such aneffect may exist 
cannot account for all the observed change flow resistance. The fact 
that the Von Karman observed even for small amounts sus- 
pended that there substantial internal effect which modifies the 
turbulence. This effect always the direction reducing the friction fac- 
tor the flow resistance. 

Mr. Laursen’s suggestion that one can study, experimentally, the effect 
moving grains resistance represents interesting extension the present 
work. This raises the whole question the effect sediment flow during 
cases where the transport rate sediment low and most the material 
being carried bed load. investigation this problem further 
light the validity the idea that the principal action the sediment the 
turbulence occurs near the bed where the concentrations are high and where 
most the turbulence generated. 

Mr. Chou was correct inferring that the slope difficult measure. 
matter fact this one the most difficult measurements make both 
the field and laboratory and the one that probably has the greatest error 
it. making the experiment great care was taken avoid end effects which 
would involve either backwater drop down curves. When sample small 
one liter was taken from the flume determine sediment concentration 
change depth water could observed immediately. This was corrected 
replacing clear water the samples were being the flume. 

The authors are inagreement with general descriptionof turbu- 
lence and turbulence production. Turbulence producedat given rate the 
one hand and dissipated heat damping given rate the other hand. 
When the rates production and dissipation are equal, steady state has been 
established presumably occurred the flume during the experiments. The 
intensity and other characteristics the turbulence ina flow are principal 
factors determining the velocity distribution. When sediment added 
flow moved it, appears that the turbulence changed either 
modifying its production its dissipation rate. This also changes the velocity 


May, 1959, Vito Vanoni and George Nomicos. 
Prof. Hydr., California Inst. Technology, Pasadena, 
Specialist-Scientific Research Engr., Republic Aviation Corp., Farmingdale, 
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distribution has been observed. true Mr. Chou observed that 
flow “turbulence has been created the expense the energy stock the 
flow.” However, once stream energy, that is, potential energy, has been con- 
verted turbulence longer recoverable produce forward flow the 
stream. view this the authors not see how “dumping” its suspend- 
load stream can release the energy that necessary accelerate the 
flow suggested Mr. Chou. depositing its load modify its 
turbulence and hence velocity but cannot make available additional energy 
produce flow. Mr. Chou states, “Only fine particles less than 0.2 can 
thrown into suspension and move with turbulence generated purely hy- 
draulic resistance.” This not borne out observation much coarser 
material than this found suspension influmes wellas innatural streams. 

Mr. Chou states that velocities are higher rising than falling stages 
stream, because the slope higher during the rising stages and less 
loss This seems like reasonable statement, but 
the rates rise, and hence the previously mentioned effects, are extremely 
small and one must look for other explanations for the observed phenomena. 
example Carey and Keller report that the duration the rising 
part the hydrograph the lower Mississippi was about months. this 
case seems reasonable assume that the effect the rate rise very 
small and that one can consider the flow being essentially uniform. 

Mr. Chou’s suggestion that one way attack and clarify this problem 
through direct measurements turbulence certainly pertinent 
points the dire need for instruments for making such measurements which 
the present time are only the very initial phases development. 

and Richardson bring very interesting and important 
problem the effect the scale the system sediment transportation phe- 
nomena. This problem that workers sedimentation studies have hardly 
touched and one that must faced results laboratory work are 
applied field conditions. Certainly one the ways study the problem 
actually work larger and larger streams. appears the 
writers that considerable understanding transport problems can also 
gained continuing study the phenomena involved. Much this can done 
small flumes which are cheaper operate and more convenient work with 
than large ones. this way laboratories that not have large facilities can 
continue make significant contributions this field. 

The authors welcome the discussions bedforms. The terminology used 
the authors that outlined the subcommittee sediment 
the American Geophysical Union. This committee did not recommend term 
describe bed without dunes other forms. The authors used the term 
“flat” and Simons and Richardson use “plane” describe sucha bed. The 
authors prefer theterm flat over the absolute term plane because according 
definition allows for slight irregularities surface which, fact, occur 
and hence, more precise term. apparent that some additional forms 
have been found and Richardson and that the problem not sim- 
ple first appeared. the list bed forms should added the “sand 
wave” observed Brooks (10) occur between the condition “dunes” and 
“flat bed” shown respectively and Fig. 1(d) Simons’ and Richard- 


Numerals parenthesis refer the Bibliography references original paper. 
Report the Subcommittee Sediment Terminology, Transactions, Amer. Geo- 
physical Union, Vol. 28, No. December, 1946, 936. 
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DISCUSSION 


son’s discussion. The information reported the effect size and grading 
the bed material the form the bed roughness most interesting. 
experiments reportedin the paper the gradationof the two materials used was 
about the same; however there was change the mean size shown 
Table results reported previously (10,18) (for instance see 
Table ref. 10) indicate that for flows the same depth and velocity 
1/2-in. flume which form dunes, the grading sand does not greatly affect the 
observed friction factor. Detailed measurements the dunes were not made 
these experiments but photographs show only modest differences size and 
shape dunes. 

Simons and Richardson question the validity ascribing all the change 
friction factor between sediment laden flows and comparable clear flows 
the damping effect the sediment. They would ascribe some the effect 
the change the intensity turbulence the two cases. The authors have 
used the idea that the significant difference between the two cases, that is, the 
one with the loose bed and sediment motion and the other with the fixed bed 
and clear water, the sediment load and that therefore the cause the ob- 
served change roughness must some way due this sediment load. 
the turbulence intensity greater one the other then the argu- 
ment would that this caused the sediment the flow the lack it, 
whichever the case may be. 

These writers also refer their very interesting findings that extremely 
high concentrations fine material suchas silt andclay have very important 
effect the behavior The explanation these effects must in- 
volve one both the mechanisms for changing the roughness factor 
stream, that is, modification the turbulence and the bed forms. 

Simons and Richardson question the contention the authors that the small- 
friction factor the rising flood can explained difference bed 
configurations between the river rising and falling stages. They are more 
inclined explainthis difference time lag the change the bed forms. 
The authors feel that because during floods the stream actually has tremen- 
dous capacity for transporting sediment and hence changing the bed configu- 
ration there would not very large lag. This especially true the case 
the Mississippi River where the times involved are measured days and 
weeks and are not only matter few hours. further reason for not ac- 
cepting the idea “time lag” that the not willing discard the 
idea advanced Brooks that order get single-valued relation for stream 
characteristics one must consider the sediment load independent im- 
posedvariable. According tothis stream wouldadjust its bed forms, 
friction factor, mean velocity and stage enable transport the sediment 
load delivered it. The bed forms and other hydraulic quantities can 
thought resulting from the sediment load instead vice versa. This idea 
can now applied explain the rating curver with loop (15). One can 
first observe that the sediment load higher during the rising than during the 
falling stage. Then according the above concept the stream will adjust it- 
self carry the load. carry the heavier load delivered the rising 
stage the stream will reduce factor the bed modifying the bed 
form which will result reducing the stage observed. The validity this 
concept and the results deduced from well the idea “time lag” 
changing bed configurations need further evidence determine their validity. 
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PRESSURE CHANGES OPEN JUNCTIONS 


Closure William Sangster, Horace Wood, Ernest Smerdon 
and Herbert Bossy 


WILLIAM SANGSTER,! ASCE, HORACE ASCE, 
The contributions the discussers are gratefully acknowledged the authors. 
Additional data have thereby been brought light which tend essentially 
corroborate the authors’ findings. 

Mr. Geyer’s comments the desirability shaping pipe entrances and 
exits wellas junction and well-foundedon non- 
hydraulic bases. The authors were, course, basing their conclusions the 
hydraulic aspects the situation, particularly applied storm sewer sys- 
tems. Extension the conclusions include sanitary sewer systems was not 
intended and certainly not advocated. 

The verification the authors’ results for equal-size pipelines Messrs. 
Blaisdell and Manson gratifying. Comparison their with those 
the authors does indeed indicate remarkable agreement. 

The answers the questions raisedin the last paragraph the discussion 
Messrs. Blasidell and Manson the assumptions inherent 
Eq. stated, was assumed that the momentum the the later- 
ineffective maintaining flow the downstream pipe. This assumption 
obviously not valid the junction width transverse the main line only 
slightly greater than the diameter the downstream pipe. Figs. 
clearly show that the size the junction materially influences the magnitude 
the loss the junction. the smaller sizes the lateral flow apparently 
more effective maintaining the flow the downstream pipe with the result 
that lesser drop pressure across the junction required. 

Eq. should not expected close agreement withresults obtained 
with pipes joined directly.5 Rather Figs. and and the supporting Eqs. 
and are pertinent these cases. fundamental differences between the 
results presented those Manson and Blaisdell’s work5 are evident. 


June, 1959, William Sangster, Horace Wood, Ernest Smerdon and Her- 
bert Bossy. 

Assoc. Prof. Civ. Engrg., Univ. Missouri, Columbia, Mo. 

Prof. and Civ. Engrg., Univ. Missouri, Columbia, Mo. 

Assoc. Prof. Agric. Engrg., Texas College, College Station, Tex. 

Highway Research Engr., Bur. Pub. Rds., Washington, 

“Energy Losses Draintile Junctions,” Manson and Blaisdell, 
Agricultural Engineering, Vol. 37, April, 1956, pp. 249-252, 256. 
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DETERMINATION HYDROLOGIC FREQUENCY 


Corrections Discussion William Sammons 


CORRECTIONS DISCUSSION WILLIAM the Febru- 
ary Journal the Hydraulics Division page 111, for 6.80 under 
column, 0.58453 should 0.68453. page 112, for 11.90 under column, 
0.51510 should 0.52510. 


July, 1959, Ven Chow. 
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CONSISTENCY 


CLARK,! ASCE.—The author has presented unit hydrographs 
and data for one the unique streams the world the form plea for 
more consistency unit hydrographs. The paper illustrates the problems 
determining unit hydrograph the originally presented methods, ‘which re- 
quire floodrecord for unit duration storm, when there areno suchrecords, 
but only records much larger multiple duration events. 

The results such application produced noteworthy unit hydrograph, note- 
worthy because its two peaks and the conviction the author that these 
peaks are characteristic the watershed. 

number years ago, the writer introduced seeking con- 
sistency derivation, utilizing data not dependent having record from 
conventionally uniform storm unit duration. was also illustrated with hy- 
drographs exhibiting double peaks, both the records and the unitgraphs, 
likewise attributed and compared watershed characteristics, particularly 
watershed shape, length, siope, and factors contributing the storage and dis- 
charge capabilities the waterway. Perhaps the problems unusual streams 
and these basic concepts stream flow chronology could illustrated for 
other students comparison derivations and discussion based such 
comparisons. 

The author kindly furnished map, reproducedin Fig. 11, showing the domi- 
nant configurations the stream, Karnafuli River Rangamati, Pakistan. 
glance the pattern suggests several possible two peak unit hy- 
drograph for short duration rainfall. One such reason, watershed shape, 
seen the relatively small amount watershed contributing the main chan- 
nel the zone from Dimagiri Barkal and the relatively large contributing 
areas above and below this zone. This effect would even more the pro- 
portions indicated the author the hydrograph considered were 
Kaptai, for that case the areas above and below the constricted zone would 
about equal and the two distinct peaks would quite reasonable. This con- 
sideration may still applicable, however, because this stream unusu- 
ally small slope, about one foot four miles, that the stage Rangamati 
may indeed function flow farther down stream. fact, low stages, 
influenced, that case effects from tides the mouthof the river, 
many miles farther downstream. But this effect question itself, that 
could quickly answered, one had measured the flow the stream large 
number times and had determined that slope was always the same 


August, 1959, Bertram Barnes. 
Hydr. Engr., Tulsa, Okla. 


“Storage and the Unit Hydrograph,” Clark, Transactions, ASCE, Vol. 110, 
1945, 1419. 
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for the same stage, that was not. One can suspect the latter, but the 
stream truly unit hydrograph stream, stream constant chronological 
regimen, would have This, course, one the many 
problems. 

For this discussion, the premise will held that the stream susceptible 
unitgraph concepts analysis withinthe range accuracy desired, pre- 
dicated the author. Thus, possible either deny the downstream ef- 
fects lump them the many causes for troublesome erratic deviations. 

Shape may factor another significant way, too. There may sig- 
nificant difference the depth rainfall and hence runoff. This difference 
may either the studied thetypical patternof the basin, 
that may not watershed orientated the sea, and storms 
influenced such rugged ridges that parallel these many long tributaries. 


FIG. 11.—KARNAFULI WATERSHED 


the typical the specific distribution happened light the narrow zone 
the watershed and heavy wider portion, there could certainly strong 
influences the direction double peaks. 

The author also furnished some daily precipitation reports for July 
August 1947, showing in. in. the downstream portion the 
shed, but actually data all for the remote areas the watershed. The 
amounts are shown there were little rainfall the ex- 
treme limit the watershed, and these large quantities the half the 
shed nearest the gage, the writer has doubt that the basic hydrograph for 
short duration would exhibit two peaks the relative magnitude found the 
author. But would have raise the question whether this non-uniform 
distribution could produce multiple unit hydrograph such device 
defined. 

There another possibility which makes this stream even more worth know- 
ing about. could have double peak with the first peak being greater than 
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the second and for uniformly distributed runoff, the basically defined distri- 
bution for unitgraph, the slope the downstream waterway much great- 
than the upper one, cuts through gorge deep bank much narrow- 
than upstream. Improbable this might seem, the Karnafuli 
shed produced some kind uplifting and folding, and this action may have 
produced long flat ponding reaches above Barkal, series rapids and gorges 
below, both. Until something more known about the topography this 
watershed, these possibilities must kept mind. These are the things that 
are determined field reconnaissance and that often escape the desk-bound 
practitioner. One can hope that the future years can yield some records 
the things which now must surmised. 

Quite obviously, the Karnafuliis not asimple stream homogeneous nature 
sedimentary bed, both because the author has described otherwise, and 
because could not exhibit the kind stream pattern that does were. 
seen from the pattern Fig. 11, there are long narrow valleys the ui- 
rections adverse general stream direction, and abundance evidence 
stream piracy and other forms channels not equilibrium with the simpler 
forces stream building. were simple sedimentary stream, would 
not exhibit unitgraph with pronounced double peak. One may glimpse the 
factors which must defined before streamflow chronology can correlated 
with watershed characteristics and some the defining unitgraphs 
consistently, comparing the hydrograph for such simple watershed with 
that which reported for this stream. 

presentedin Fig. 12, the writer has deriveda unit hydrograph, storage- 
based instantaneous unit hydrograph (IUH), compatible with the furnished re- 
cords this stream, and the approximate shape this watershed defined 
this map, Fig. uniformly distributed unit quantity 
generated runoff. Engineers given the same watershed map, which di- 
mensionless, the same storage constants, and the same dimensional value 
the total watershed area, could get the same answer, achieving the independence 
from the personal equation sought the author the present paper. 

Comparing this uniform IUH functionfurther with the author’s derived unit- 
graph, the indicated peak near the magnitude found the author, and the es- 
sential duration nearly that found the author. Few currently active pro- 
cedures concern themselves with more, and work can done with 
this much accuracy and agreement. 

The agreement goes farther. Both falling legs have similar rates re- 
cession. The author has placed considerable emphasis the recession rate 
making his derivation, and the users the instantaneous unit hydrograph 
procedure recognize that one the two storage constants basic that method 
measure this phenomenon. With much agreement, almost axio- 
the gravity the hydrographs and their lags behind rain- 
fall are almost identical. 

There still remains the questions the exact hydrograph configuration and 
the sharp double peaks. The writer was genuinely disappointed that the basin 
shape, divided with real knowledge the location the watershed 
ridge lines, and many doubts about tributary connections, did not support the 
double peaks. 

Indeed, the enlarged lower basin comparison with the central zone near 
Barkal shows only shoulder the rising hydrograph, not peak. That 
which shows shoulder for uniformly distributed short duration runoff 
could peak for some non-uniform rainfall distribution, perhaps, particu- 
larly with little runoff from the remote portions the watershed. But 
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the rainfall distribution vaguely suggested the few numerical values rain- 
fall shown also Fig. 11, amounting in. somewhat above the mid-area, 
although only 3.1 in. the mid-zone, has tendency the opposite direction, 
shown the non-uniform IUH function also shown Fig. non- 
uniform IUH function computed much the same manner the uniform 


Derived Unitgraph 


Non-uniform 


Flow 1000 cfs per inch 


nN 


Days after beginning Runoff 


FIG, 12.—COMPARISON UNITHYDROGRAPH FUNCTIONS 


IUH function, but with the modification Victor account for 
varying depth runoff. 


The author’s search for consistency unitgraphs hinges two significant 
procedures derivation and the reduction the result dimensionless 


“Storage and the Unit Hydrograph,” Clark, Transactions, ASCE, Vol. 110, 
1945, Discussion Victor Cochrane, 1461. 
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basis. The procedure derivation involving determination volume moments 
the hydrograph, and the rainfall volume, and their difference mea- 
sure the phase difference, time lag, between runoff and stream flow 
sound mathematics and growingly popular procedure. has directly com- 
parable relation procedure storage-based IUH functions. 

The procedure developing the unit hydrographfrom multiple unit event 
unknown non-uniformity runoff distribution can necessary expedient, 
but mathematically indeterminate procedure principle, capable 
indicating obvious errors like the negative flow values treated the author, 
and producing large number sometimes conflicting sequences num- 
bers which evade the consistency sought. Only the author’s statement that the 
double peaks his unitgraph “were present all eight unitgraphs that were 
derived,” and the assumption that this refers unitgraphs from eight 
ent flood events, leads the writer hope for future evidence the watershed 
causes this indication. Some years ago the indeterminate aspect this 
procedure was hinges around the fact that Eq. when run- 
off sequence known, there are more simultaneous solution than 
there are unknowns, and one set unitgraph sequence values Ug, Ug... 
for each combination equations can found, all which compare with 
each other differently the forward and the reverse runs presented. The 
first and the last series the available equations make the “forward run” 
and the “reverse run” and are theseries most commonly selected for this kind 
analysis, but the intermediate series are available, and neither more nor 
less accurate than the initial and terminal series. appears probable that 
the solution values oscillate backand forthacross the correct solution some 
manner. The writer would unable estimate measure the amplitude 
this oscillating deviation, and has found large and occasionally em- 
barrassing. 

The author’s addition the procedure the recession value adds another 
equation those available, and another series solutions those considered 
above. While this addition increases the number solutions available, there 
little doubt that improves the quality all solutions substantially re- 
ducing the amplitude deviation about the correct solution the unitgraph. 

this process simultaneous equation solution extended very short 
units time, fractions the 24-hr unit illustrated, the number equations 
for solutionand the workof solution multiplies degree which precludes its 
performance until rapid electronic computers became available. The instan- 
taneous unit hydrograph functions, IUH functions, with their infinitesimally 
short unit periods and continuous curve solution, were natural, but 
evolution from these difficulties. 

The author recognizes lag, time-phase difference between runoff and 
the hydrograph, characteristic likely define hydrographand suscepti- 
ble successful correlation between watersheds and topographic measure- 
ments. Successful correlation with those watershed dimensions indicative 
storage and discharge capacity has been obtained, and correlation with slope 
also insofar this item can define storage, discharge, similarity 
sheds. The measures lag mentioned the author refer general toa 
mean-time displacement the centers centroids hydrographs, rainfall, 
and watersheds. Those contemplated IUH functions are measures the 
similar displacement the most remote elements rainfall and watershed. 
“Time concentration” the oldest presently used name representative 
such concept, and the one most likely lead consistency over the widest 
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range possible Presumably, for rectangular water 
shed, the former group displacement measures would about one-half 
the latter, and values thirdand two thirds the latter might apply 
the two kinds triangular watersheds whichdefine reasonable limits with- 
which orthodox watershed shapes lie. 

The author seems have selected 108 his measure lag. 
The illustrated IUH function for time concentration 144 hr. Hence the 
former two thirds the latter. The basin shape such that the centroid 
the watershed beyond the half-length the watershed. Greater consis- 
tency than this must await better means defining and measuring the items 
and concepts under discussion. 

The author defines the recession the Karnafuli unitgraph the dimen- 
sional time reciprocal, 0.69 per day. The writer defines the IUH recessional 
the storage-discharge parameter 2.75 days, for which the reces- 
sional decrease measured the ratio values apart 0.692. Ap- 
parently these are consistent, and are merely two ways writing the same 
thing. 

The author approaches consistency between watersheds via the dimension- 
less hydrograph. Merrill Bernard’s distribution graph had this merit and 
seems odd that did not survive place the unit hydrograph. Perhaps 
may revived for its real merit. The IUH function approaches the same end 
via the dimensionless map watershed shape diagram, commonly called 
distance-area curve. The shape diagram and the two storage parameters, 
days and 2.75 daysfor Karnafuli River above Rangamati, area about 3800 
miles, define the IUH function, either the practically dimensionless units 
percentage total volume per hour, English units cubic feet per 
second per watershed inch volume. For this particular watershed the re- 
sults have been shown similar many respects, highly consistent 
either, both, any criterion consistency. 

The author demonstrates application which also includes the much 
smaller Castle Creek having one thousandth the watershed. For IUH functions, 
storage values min and min have been found define instantaneous 
unitgraph for 1/4 acre near Temple, Texas, which dry the time. 
Values days and 4.5 years similarly define the wholly theoretical unit- 
graph the 300,000 miles involving the Great Lakes the St. Lawrence 
River, where backwater effect from downstream points extends almost two 
thirds the waterway length low water, and 30% the watershed inun- 
dated low water. The seeker consistency may find equally these ap- 
proaches and over wide ranges possible use. The seeker greater con- 
sistency must look also for more precise words, thoughts, and concepts, and 
first conceive and then measure the much smaller deviations which these 
two approaches are not entirely identical and not completely consistent. 

All not alike, and some cannot understoodin the exact terms 
unitgraph, but that item the best presently available tool for measuring 
concept. Even where unitgraph not therule the exact nature par- 
ticular stream, the unitgraph the best hope yardstick about which and 
from which some other behavior can defined consistent manner. 
understanding the factors presented the author against the background 
the reader’s observed personal experience will lead the consistency which 
has been attained only among civil engineers and within the Society. More 
records from this unique stream Pakistancan help close remain- 
ing gaps complete path consistency unitgraphs. 
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REVISED COMPUTATION VELOCITY HEAD WEIGHTED 


Corrections Discussion Manual Benson 


CORRECTIONS DISCUSSION MANUAL BENSON, ASCE.— 
the February Journal the Hydraulics Division page 128 the following 
corrections must made: the equation the top the page the term 
should changed Table Col. change the head- 
ing “Conveyance, “Conveyance, k.” 

Table the notation “Approximate Values of” that presently runs over 
columns through should only encompass columns through 
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WAVE INDUCED MOTION BOTTOM SEDIMENT 


Discussion Robert Miller 


ROBERT MILLER.!—The sequential papers Ippen and Eagleson2 
Eagleson, Dean and Peralta? and that under present discussion, have direct 
including the present writer. The interaction shoaling waves and sediment- 
fluid interface the near shore region the oceans presents very complex 
problem even under greatly simplified conditions. Specific questions this 
category include: equilibrium conditions for shorelines, particular, beach- 
es; the mechanism entrainment and subsequent deposition sediment 
particles; and the resulting statistical patterns average (or median) size 
and the variance (sorting) the sediment particles. 

The first two papers referred offer, the marine geologist, the 
tunity test and interpret the null point hypothesis directly nature. This 
was done analysis detailed field data collected specifically for this pur- 
pose. The analysis was made number field localities from 
the coast California, Lake Michigan, the New England coast. 
Expectationof the areal pattern for median size, and sorting (using the 
ficient variation) was developed from the empirical relation given Ippen 


and Eagleson2 as, 


assuming initial size distribution consistent with geological expectation. 
The final result was compared with field data subsequently smoothed least 
squares. Consistently good agreement with expectation was found, provided 


Assoc. Prof., Dept. Geology, Univ. Chicago, and Assoc, Marine Geology, 
Woods Hole Oceanographic Inst. 

and Eagleson, Hydrodynamics Lab. Tech. Report No. 18, 1955. 

Mechanics the Mition Discrete Spherical Bottom Sediment Particles due 
Shoaling Waves,” Eagleson, Dean and Peralta, Beach Erosion 
Board Tech. Memorandum 104, Washington, 1958. 

The average size, and the variance around the average sample sediment 
particles, constitute two very useful and widely used parameters Geology. 

Model Relating Dynamics and Sediment Pattern Equilibrium the Region 
Shoaling Waves, Breaker Zone and Foreshore,” Miller and Zeigler, 
Journal Geology, Vol. 66, No. July, 1958. 

Study the Relation Between Dynamics and Sediment Pattern,” Miller 
and Zeigler, Eclogae Geologicae Helvitiae, Basel, Switzerland (in press). 

“Comparison Theoretical Mass Transport Velocities with Observed Sediment 
Size and Sorting Patterns,” Miller and Zeigler, First Internatl. Oceano- 
graphic Congress Preprints, New York, Y., 1959. 
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certain simplifying assumptions were made and limiting conditions satisfied. 
The conclusion that the relationship expressed Eq. appears hold 
nature under the restrictions referred above. 

their present paper Eagleson and Dean present two rational equations for 
predicting median sediment diameter versus bottom slope, that are based, 
part, their experimental data, and part theoretical considerations in- 

view the encouraging results achieved applying the rational equation 
(Eq. natural conditions, planned repeat the procedure with the new 
prediction equations Eagleson and Dean. 

However, there remains one important aspect that felt has been neg- 
lected Eagleson and Dean. There definite patternto the sortingor vari- 
ation the sediment the near shore region, well the median size 
distribution, the natural cases which were studied. 


Eagleson and Deanassume that most nearly represents the 


cle geometry equilibrium natural beaches, which the sediment parti- 
cle diameter, and the mean diameter the roughness particle. This does 
not appear justified nature. The ratio may more realistically ex- 
pressed average natural beach, but will also vary function 
sorting. Since there definite gradient the sorting pattern, average 


will exhibit the greatest divergence from unity, the seaward direction. 


The movement particles down-slope direction seaward the vari- 
ous null points has been recorded under laboratory conditions. The existence 
this phenomenon, any appreciable degree nature, open question. 
The following argument given support this contention. 

Suppose start with initial uniform particle distribution. Then, null 


points for the various sizes are subsequently established under wave motion, 
two possibilities are considered: 


Particles smaller than given null point move shoreward. Particles 
larger than given null point move seaward. 

Particles smaller than given null point move shoreward. Particles 
larger than given null point have appreciable net movement. 


the first possibility, the implication that the null point itself the 
particle sizes remaining willall point diameter. evalu- 
ation variation dimensionless form, that is, standard deviation mean 
deviation, should result approximately equal values over traverse from 
beach out sea. 

the second possibility, the that the ratio standard deviation 
mean should show definite gradient, from high values seaward low values 
shoreward. 

Observations nature and those others such indicate defi- 
nite gradient from the dimensionless sorting coefficient seaward 
low values shoreward. Thus appears correct, and seaward movement 
particles down slope from null points appears negligible. 

“Mass Transport Water Waves,” Longnet-Higgens, Philosophical Transac- 
tions, London, Ser. No. 903, Vol. 245, March, 1953, pp. 535-581. 

“Areal and Seasonal Variations Beach and Nearshore Sediments Jolla, Cali- 
fornia,” Inman, Beach Erosion Board Tech. Memorandum 34, 1953. 


Be 
| 
3 
> 
| 
4 
| 
| 
| 
- 
ES: 


DISCUSSION 131 


the implications the gradient sorting are considered, then, the 
average D/k varies from the value with the greatest divergence sea- 
ward direction. such statistical chance that particle move down- 
slope under the influence gravity minimizedin nature, due the low like- 
lihood particle removal from “pockets” the interface, and the high like- 
lihood immediate retrapping particles “pockets” when they are moved 
downslope direction. D/k approaches unity the similarity between the 
natural condition and the experiments performed the Massachusetts Institute 
Technology will become stronger. 

summary, the two new prediction equations given Eagleson and Dean 
appear considerably refined over the earlier form given 
Ippen and Eagleson. particular importance for field investigation the 
present inclusion beach slope variable, and the implicit inclusion 
mass transport velocities. 

The writer’s collected field data will reanalyzed terms these new 
prediction equations with the expectation considerable improvement our 
results. felt that Eagleson and Dean have made distinct contribution 


toward solution the complex problem sediment movements versus wave 
motion the near shore region. 
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HYDRAULIC CHARACTERISTICS GATE SLOTS# 
Discussions Roy Tinney, Robertson and Bennett 


ROY TINNEY,! ASCE.—The analysis flow the vicinity 
gate slot has considerable practical significance even though the slot itself 
detail large structure. Excessive vibration the gate leaf and 
cavitation damage the downstream conduit walls can become serious 
ational hazards. 

The problem essentially one determining slot shape that will ac- 
commodate wheel, roller, guide without affecting the flow. Although this 
might seem straightforward problem, reality becomes somewhat 
complex. The fundamental problems involved are: (1) the spreading jets 
and, (2) flow separation discontinuities the boundaries. The offset dis- 
tance depends the rate spreadof the jet into the slot and the shape the 
downstream corner must designed that the impinging high velocity jet 
does not tend separate from the downstream wall. 

One experimental approach this problem study the flow patterns 
the neighborhood the slot with particular attention the spreading the jet 
and separation the flow. The author however, has chosen the alternate ap- 
proachof analyzing pressure patterns for several the most practical designs 
with systematic variation the dimensions and configurations. The suc- 
cess this latter approach has been well demonstrated. evident from 
the data presented that gate slot with small offset and very gradual tran- 
sition from the downstream corner the line the downstream wall 
cavitation-free even extremely high heads. The type transition does not 
12-in. radius arc the downstream endappears give results comparable 
with those for large radius, circular arc convergence (Fig. 14). Presuma- 
bly, ellipse would equally satisfactory since gives acceptable 
solution even when projects into the flow (Fig. 16). These findings are 
accord with designs for other cavitation-free bodies subject impingement 
high velocities, for example, baffles spillway aprons, nose-bodies for 
torpedoes, struts, etc. 

The writer questions the interpretation placed the data onincipient cavi- 
tation. particular, the writer disagrees with the statement that “the pres- 
sure the surface just downstream the downstream slot corner neednot 
vapor pressure magnitude for cavitation occur the slot,” and with the 
idea that “scale effect” exists. Two other interpretations are suggested that 
account for the apparent anomalies. First, the average pressure indicated 
piezometer-manometer tube apparatus may considerably higher than 


Head, The Albrook Hydr. Lab., Div. Industrial Research, Washington State 
Univ. 
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vapor pressure and yet the pressure may fluctuate strongly that will 
vapor pressure for significant percentage the time. Indeed, has been 
that model data may indicate average high -14.5 
water and yet serious cavitation damage occurs the prototype. The use 
pressure transducers and suitable recording instruments measure pressure 
fluctuations show that this same point the instantaneous pressure was less 
than -30 water for 26% the time. 

second explanation for the apparent discrepancy between relatively high 
pressures measured the model and the occurrence cavitation the 
extreme difficulty finding the minimum pressure when the pressure gradients 
are large. Fig. 11, for example, not all clear that the minimum 
-0.30. this latter value used with the data Fig. for the square cor- 
ner (hy 146, 14), the minimum value -30 ftof water 
well with the plotted data. does not appear, therefore, that “scale effect” 
exists the author suggests. 

The alternate method proposed the author using cavitation index 
determined special model tests course presupposes the absence “scale 
effect.” This technique does not require detailed pressure measurement but 
does require careful visual observation. compilation such indices for the 
would greatly simplify the designof cavitation-free 
gate slots. the meantime, the application the excellent set guide lines 
presented this paper most helpful and the author commended for 


the diligent effort that has made over long period time gather this 
wealth material. 


J.M. ROBERTSON, ASCE and author presents 
thorough discussion the hydraulic problems and their possible solutions, 
resulting from the necessity having gate slots large conduits. His test 
results and observations represent valuable extension our knowledge 
the flow past slots and generally agree with other studies. This dis-- 
cussion based primarilyon the writers’ experience with slot problem 
large hydraulic structure where the experimental results supplement the au- 
thor’s. The necessity for information the boundary-layer thickness the 
flow approaching the slot emphasized. 

The paper discusses two basic hydraulic gate-slot problems—‘turbulence” 
within the slopes and cavitation just downstream the slots. Although 
these may result entirely different damage effects, they are not unrelated. 
Both are intimitely associated with the flow conditions and near the slot. 
The pattern (that the author the slot func- 
tion the slot proportions—primarily the width depth ratio—as suggested 
the author’s observation that great amount circulation occurs the 
slots and that this varies with slot width. view other observations the 
flow patternin slots, the comment that the circulation seems greater with wide 
slots than with narrower ones must considered tobe limited applicability. 
Certainly the circulation vorticity small for very wide slots, since 


Model Studies Noxon Rapids Stilling Basin,” Howard Copp and 
Roy Tinney. Unpublished report The Albrook Hydr. Lab., Div. Industrial 
Research, Washington State Univ., February, 1960. 

Prof. Theoretical and Applied Mech., Univ. Urbana, 

Supt., Riverside Paper Corp., Appleton, Wis. 
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such case the flow appears expandinto the slot with definite vortices only 
near the ends. 

The nature vortex occurrences slots such that one must careful 
interpreting the results experiments conducted few width-depth 
ratios. The vortex structure and some the resulting effects may vary ina 
nearly discontinuous fashion with the slot proportions. For certain slot pro- 
portions (W/D ratio about one, one-half, three) rather stable patterns ap- 
pear likely. With proportions not greatly different, the vortex pattern not 
definite and considerable flow possible. unstable pattern ap- 
pears likely for depth width ratio about two. For rather deep slots the 
vorticity small and the flow disturbance negligible. This the rational 
explanation for the commonrule that piezometer hole for pressure measure- 


STABLE TRAPPED VORTEX 


PATTERNS 
W/D= 172 


W/D 


FIG, 1.—INTERNAL FLOW PATTERNS SLOTS VARIOUS PROPORTIONS 


ment should least twodiameters deep before changes size. Fig. 1(a) 
indicates possible stable vortex patterns for this and shallower (that is, wider) 
slots. Very few the vortexpatterns inslots are avail- 
able, although numerous authors indicate the stable trapped vortex found with 
ratio about unity. Thom and Appelt have recently the 
laminar -viscous flow pattern lowslot Reynolds numbers indicating the man- 
ner which the dividing stream line enters into the openings. connection 
with the flow past rough surfaces, when the protuberances re- 


c 


“The Pressure Two-Dimensional Static Pressure Hole Low Reylonds Num- 
bers,” Thom and Appelt, British Aero. Res. Council, ARC18798-FM2473, 
February 20, 1957. 
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sult large W/D values, few experimentors have pattern 
one, possibly two, vortices behind protruding obstacle with following 
region disturbed flow into which the outer flow intrudes. This general 
agreement with the writers’ observations® indicated Fig. 1(b). Folsom 
reported motion pictures the unsteady pattern showing successive for- 
mation and sweeping out vortices. his studies, the maximum energy loss 
occurred when the width was about twice the depth. One can rationalize that 
vortex pair pattern, which might expectedfor W/D about two the scheme 
stable trapped vortices indicated not possible due the man- 
ner which the flow mustoccur the upstream and downstream corners. 
least one other flow situation, which the writers are cognizant, has exhibit- 
almost discontinuous behavior (as evidenced loss coefficient) 
certain W/D ratio. The careful studies? Roshko show some discontinuities 
slot flow W/D ratios the vicinity 

Returning the observations the pressure distributions immediately 
downstream slots, shownin the author’s one might question whether 
the curves indicated might not actually contain discontinuities. Since the maxi- 
mum disturbance (from energy loss measurements) appearsat width depth 
ratio about two, some discontinuities the pressure curves might 
pated near this ratio value. That such are not evidence, either the curves 
shown the author presented later this discussion, may only due 
incompleteness the experimental programs from which the data was ob- 
tained. Certainly anyfuture experimental studies this problem should verify 
the possibility discontinuous behavior. 

the author indicates, the effect offsetting the corners the slots may 
either beneficial adverse depending upon the relative offset direction. 
Evidence offset effects obtained the Iowa Institute Hydraulic Research 
has that the upstream corner the pressure dipis quick- 
reduced while raising the downstream corner increases the dip. This 
qualitative agreement with the author’s Some interesting “field” 
prototype experience slot offset effects was obtained with the Garfield Thomas 
Water Tunnel the Pennsylvania State University. this research facility, 
water speeds tosome fps flows through 4-ft cylindrical section (about 
long) which hydroballistic shapes are tested for their cavitation and 
removable section cover the 4-ft conduit. The necessary 
ance slots the upstream and downstream edges this cover, although only 
about 1/16 in. wide, were great concern the design the tunnel due 
their being possible sources cavitation. Since offsetting was expected, 


“Frictional Resistance Artificially Roughened Pipes,” discussion Fol- 
som, Transactions, ASCE, Vol. 101, 1936, 707. 

“Relation Statistical Theory Turbulence Hydraulics,” discussion Ne- 
menyi, Transactions, ASCE, Vol. 105, 1940, 1582. 

Past Slot Surface,” Bennett, The Pennsylvania State Univ., 
Master’s thesis, 1949. 

3488, August, 1955. 

“Use the Low Velocity Air Tunnel Hydraulic Research,” Rouse, Pro- 
ceedings, 3rd Hydr. Conference, Studies Engrg., Bulletin 31, State Univ. Iowa, 1947, 
132. 


«Cavitation Sluice Gate Slots,” Fetouh, Master’s thesis, State Univ. 
Iowa, 1947. 
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design studies indicated that any problem could easily round- 
ing the downstream corners. 
Studies the cavitation characteristics the completed tunnel indicated 

the occurrence premature cavitation. This was finally traced offset 

the seating the cover, resulting from improper gasket construction. 

The offset was measured about 0.050 in., which rather small amount 
the 4-ft diameter. was the boundary-layer dis- 
placement thickness, that was about 0.11 the front the cover where the 
offset was out the flow and 0.24 in. the end where the offset was into the 
flow. After reworking the gasket eliminate this offset, the test section 
cavitation index number was found!2 reduced from 1.1 value 
about 0.3. probably coincidental that this final value the cavitation in- 
dex close the author’s value 0.28 for slot. One very 
important aspect the water tunnel slot cavitation occurrence was that was 
only evidenced, and presumably only occurred, speeds above The 
implications such discontinuous speed effect are serious when one con- 
siders extrapolation low-speed cavitation tests structures other 
scales.13 

The evaluation cavitation potentialities structures from experiments 
involving wall-pressure measurements are indicated the author his dis- 
cussion the results found with the standard slide gate. The minimum wall- 
pressure coefficient was found not the same the cavitation index, pre- 
sumably found visual observationof Such discrepancy 
may several factors, particularly those involving the effects physi- 
cal properties and occurrences, such are often lumped under the cognomen 
“scale One factor, which not really scale effect, that cavi- 
tation probably did not occur along the wall where the pressures were mea- 
sured, but inthe separated flow region away from the wall. This 
may not occur for well rounded contours which inhibit separation. Consider- 
ation the studies other hydraulic structures susceptible cavitation sug- 
gests that the best solution not always found reducing the magnitude the 
wall pressure dip, since this may lead cavity collapse the surface rather 
than off the surface (that is, the water). Cavitation collapse the separated 
flow off the surface may result less cavitation erosion, even though more 
severe cavitation occurs. Thus have the problem correlating our con- 
clusions based wall-pressure measurements with cavitation inception and 
erosion occurrences. Until these several problems are sorted out, must 
temper our use wall pressure measurements with qualitative estimates 
the nature the flow behavior and how may affect cavitation. 

Most experimenters hydraulic structures, including the author, appear 
neglect one factor whichcan have significant effect ontheir results and the 
proper interpretation thereof. This the nature the boundary layer flow 
near the contour they are studying. present case, prototype conduit 
structures not have uniform velocity profiles nor the flow always fully 
developed. Inthe laboratory studies consideration appears have been 
given simulating, even evaluating the effect of, the nature (primarily the 
relative boundary-layer thickness) the approaching flow. will shown, 


“Performance Characteristics inch Water Tunnel,” Robertson and 
Ross, Engineering, Vol. 184, July 19, 1957, 76. 

“Similitude Cavitation-on Several Laws Cavitation Scaling,” Robert- 
son, McGinley and Holl, Houille Blanche, September, 1957, 540. 
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the boundary-layer thickness slot can have significant effect the flow 
and pressure distribution. Relatively thick boundary layers tend smooth 
over sharp discontinuities contour. Thin layers not and hence result 
large effects approaching those predictable from potential flow analyses. 
the affects boundary-layer thickness the wall pressure 
distributions near wall slots was obtained® some time agoas part 
studies for the large high-speed water tunnel already referred to. The 
concern the prototype was that around the access cover the test section. 
The experiments were performed slot models much larger than the proto- 
type order that the flow could studied detail. small single-pass air 
tunnel was used for the study with velocities outside the boundary layer about 
The slot widths were in. in. leading testing slot Reynolds 
numbers about the same values expected prototype. The tunnel 


0.2 0.4 0.6 


FIG, 2.—BOUNDARY LAYER VELOCITY PROFILES 


the 3-in. wide floor. Since major concern was the boundary-layer 
thickness, three different boundary layers were natural and 
two artificially thickened. The natural layer, here termed the thin layer, was 
about 1.2 in. thick with displacement thickness about 0.14 in. The other 
two boundary layers had velocity profiles that agreed reasonably well with those 
for natural profiles; the ratio displacement momentum thickness, con- 
ventional indicator the nature the profile, was 1.14 and 1.29, not far from 
the usual value about 1.3. The medium boundary layer had thickness 
about 5.5 displacement thickness 0.64 in., while for the thick layer 
these thicknesses were about in. and 1.28 in. The three velocity profiles are 
shown Fig. three boundary layers, conjunction with the use 
slot widths in. in., permitted study the problem over wide range 
boundary layer thickness. Besides the flow pattern observations al- 


ready noted, the wall pressure distribution and near the slot was also ob- 
served. 
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Wall pressure distribution results similar those shown bythe author were 
found for various combinations slot depth, width, and boundary-layer thick- 
ness. Tests were made with square-edged and rounded downstream corners. 
all cases the upstream corner was sharp and without offset. The effect 
slot width, boundary-layer thickness, and corner rounding were clearly evi- 
denced. Over range inwidth depth ratios 1/4 the slot depth appeared 
have little shown Fig. The pressure coefficient used here 
similar the coefficient plotted the author except that the measured 
pressures have been corrected for the hydraulic gradient the conduit (that 
is, the test the dynamic pressure which the pressures are di- 
vided based the velocity outside the boundary layer. The proper sign 


THIN AYER, W=2 
0.07 


FIG, 3.—EFFECT RELATIVE SLOT WIDTH 
WALL PRESSURE FOLLOWING SLOT 


convention (in contrast the author’s inverted one) has been used with this 
coefficient. Thus, pressures below the hydraulic grade line appear nega- 
tive coefficients. Although the smooth curves drawn through the points 
Fig. appear contravene the earlier discussion slot proportions’ effects 
and discontinuous occurrences, they appear reasonable. Actually, more 
continuous variation should have been studied. The lack significant 
slot depth effect, noted, differs from the author’s observations shown 
Fig. 6F. will discussed below, pressures’ measurements downstream 
sharp corners probably only indicate the minimum they are within dis- 
tance about 0.01 slot widths the corner. This implies that the author’s 


piezometer number only indicated the minimum pressure when was in. 
and larger. 
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Slot widthand boundary-layer thickness had marked effect the pressure 
distributions. This clearly shown for one the rounded corners Fig. 
For the wider the two slots the pressure patterns have maximum located 
about 45° around the curve while for the narrower slot the maximum occurs 
further around the curve (in the flow direction), that is, about 70°. The ef- 


V 
90° 
BOUNDARY LAYER 
THIN 
MEDIUM 
THICK 
o° 
FIG, 4.—PRESSURE DISTRIBUTION AROUND RADIUSED 


CORNER END SLOT 


fect boundary-layer thickness seen modulating the amplitude 
the pressures. The higher pressures occur with the thinner boundary layers. 

Proper determination the minimum wall-pressure coefficient requires 
that the first piezometer following the corner lies the region minimum 
pressure. For the sharp cornered slots the extent this region must de- 
fined the separation the flow from the corner. This turn dependent 
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the location the point streamline separating the outer flow from 
that trapped the slot intersects the downstream wall the slot. This point 
stagnation point whose location indicated the location the maximum 
pressure plots such that Fig. Although much more difficult 
vation, similar positive pressure peak appeared for square cornered slots. 
The distance this stagnation point below the upper surface the slot 
proportional the width the slot, noted connection with Fig. ap- 
and rounded corners. From study various other flows with separation ap- 
separation should extend along the wall the flow direction 
distance about 1.5 times This suggests that the separation (and hence 
minimum pressure) zone extends x/W about 0.01. For the tests re- 
ported this discussion the minimum pressures were thus only indicated 
the first (at 0.013 in.) when was in. general, 
great care must exercised interpreting pressure measurements made 
fixed distance from the corner the slot the widthof the slot changed. 
thus appears that some the author’s results may misleading due the 
above consideration. 

Taking data from the first piezometer following the in. and in. wide 
slots representing the minimum pressure, possible indicate the ef- 
fect relative boundary-layer thickness the pressure dip for the sharp 
corner. This shown the upper plot Fig. together with some the 
results for the rounded corners. The importance the boundary-layer thick- 
ness determining the minimum pressure following slot, particularly 
sharp cornered one, clearly apparent. The thinner the boundary layer rela- 
tive the slot width, the larger the magnitude the pressure dip. 

The suitability providing some downstream corner radii order re- 
duce the magnitude the pressure dip indicated the lower plot Fig. 
This shows clearly the importance small roundingof the corner allevi- 
ating the pressure dip. should noted that only simple radiusing the 
corners was the tests reported here. The author shows clearly the 
effect curving this corner with transitions more gentle thanthe simple circu- 
lar arc. Asubstantial decrease inthe magnitude the pressure dip obtained, 
but secondary dip occurs curve did not joint the wall the con- 
duit properly. The water tunnel design studies previously referred included 
study transition curves!4 which was indicated that contours must 
utilized which have continuous first and second derivatives. Thus cubic tran- 
sition curve desirable. The first derivative should continuous avoid 
separation while the second should avoidpressure discontinuities pre- 
dicted potential flow theory. Electrical analogy studies conduit 
Rouse and Hassan also indicate the desirability cubic transition curves. 
course just how sophisticated curve necessary depends upon the rela- 
tive boundary-layer thickness. The cubic may only needed for very thin 
boundary layers such will encountered near conduit entrances. 

this discussion does nothing else, hoped shows im- 
portance the effect boundary-layer thickness the slot problem. Thus 
many the author’s results (those Fig. 20, for example) are uncertain 
absolute validity. The boundary-layer thickness for given test situation 

“Water Tunnel Transition-Section Flow Studies,” Ross, Ordnance Research 
Lab. Report NOrd 7958-285, The Pennsylvania State Univ., 1954. 

“Cavitation-Free Inlets and Contractions,” Rouse and Hassan, Mechan- 
ical Engineering, March, 1953, 213. 
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FIG, 5.—EFFECT BOUNDARY LAYER THICKNESS AND 
CORNER RADIUS MINIMUM WALL PRESSURE 
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much part specifying the test conditions the width the 
hoped that the author will ableto indicate some estimates the relevant 
thicknesses. Thus, the test data leading the curves Fig. 
was presumably obtainedat one station conduitor tunnel known geome- 
try. this had well-rounded entry and known length straight section pre- 
ceeding the offsets, should possible estimate the approximate 
layer displacement thickness. would seem from the information presented 
this discussion that any such estimate would make Fig. considerable 
absolute value whereas now little such value. Actually, course, the 


effect variations boundary-layer thickness this particular occurrence 
should studied. 


a 
. 
= 
: 
BE. 
| 
3 
area 
all 
x 


48 
4 
a 
& 
By 
° 
: 
1 
4 


CAVITATION DAMAGE ROUGHENED CONCRETE 


Discussion Rylands Thomas 


RYLANDS THOMAS, relationship between shear velocity 
and pressure head (Fig. provides useful indication the conditions under 
which cavitation likely occur. The margin safety against cavitation 
provided static pressure clearly illustrated. 

There is, course, clear-cut “threshold” cavitation. Turbulence 
phenomenon and therefore the onset turbulent stream 
discontinuous process which could defined frequency occurrence. 
may that the build-up cavitation considerable proportions rapid 
for practical purposes diagram such all that required. 

The tacit assumption the applicationof results obtained conduit in. 
deep larger conduits (see Fig. 10) that the instantaneous velocity the 
boundary the onset dependent only the boundary roughness 
and the shear velocity that point, and not flow parameters relating 
the whole cross section, such for Reynolds Number mean flow. could 
pointed out that the conditions slightly away from the boundary iden- 
tical because the rate decrease shear velocity with distance from the 
boundary greater small conduit than large one. Nevertheless, the 
assumption appears justified observational data the intensities 
velocities turbulence near the boundaries rough and smooth conduits, the 
root-mean-squares which bear constant ratio the shear velocity for 
any given boundary surface. 

The curves shown Figs. and relate only two specimens eroded 
concrete surface tested the author. Useful these are, the designer 
high-velocity tunnels and open conduits requires information relating the 
onset cavitation before erosion the surface has occurred. particular, 
requires information smooth concrete surfaces with protuberances 
steps due joints formwork and slight waviness due lack alignment 
forms. hoped that the author’s work will extended cover these practi- 
cal cases. 


November, 1959, Donald Colgate. 
Cons., London, England. 
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MOUNTAIN CHANNEL TREATMENT LOS ANGELES 


Discussion George Newhall and Frank Henry 


GEORGE NEWHALL,! ASCE, and FRANK Ferrell’s 
paper gives very clear background and summary the objectives, benefits, 
costs, and methods the valuable channel-treatment program underway 
Los Angeles River Watershed, California. 

This writer, having viewed the Brand Canyon, Cooks and Monroe 
Canyon works, and part the other works and reference reports, offers the 
following comments and questions with view advancing few small facets 
the great and largely inexact science watershed management engineering. 

Increased deep percolation ground water could be, Mr. Ferrell point- 
out, very important resulting benefit such channel treatment projects 
areas fractured porous rock formations. conceivable that, the 
locations where maximum significant deep percolationoccurs can found, 
the amount deep percolation could further beneficially increased, rela- 
tively small cost either increasing the height dam at, just down- 
stream from, that location, drilling tunnels wells from the channel 
just upstream from dam, preferably before the reach channel fills with 
debris. Such tunnels wells could horizontal, vertical, inclined fol- 
low the strike and dip rock strata and fault zones, and could extend one 
more directions from the channel, depending the geology, hydrology, and 
economics the project. Have such features been considered, investigated, 
planned with this program mountain channels Los Angeles 
County? And have estimates been made probable increased amounts deep 
percolation? 

The significant amount “inflow water” within many water supply tunnels, 
such Tecolote, Mission Doulton (Juncal) Tunnels the Santa 
Barbara area evidence that valuable amounts usable water can and 
come out deep percolation below some mountain watersheds. Amounts from 


mgd mgd flow fairly steadily throughout the year, year after year, with 


more inseasons following heavy precipitation seasons, and less drought pe- 
riods. seems logical that this “deep percolation groundwater” could 
artificially increased through tunnels wells barriers, just can 
extracted from tunnels wells springs. 

For such “recharge” purposes, floodwaters would need de-silted 
relatively free silt, has been discovered experience Los Angeles 
County Flood Control District and various irrigation districts with recharge 


November, 1959, Ferrell. 
Civ. Engr., Forest Service, Santa Barbara, Calif. 
Civ. Engr., Forest Service, San Francisco, Calif. 
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well operation. Such de-silting debris load reduction flood waters ap- 
parently accomplished some extent automatically channel stabilization 
projects suchas Brand Canyon and Cooks Canyon decreased ve- 
locities flood flows. This clearer water may inherent benefit such 
projects, due partly resulting natural revegetation channels, and leads 
the conclusion that any proposed “recharge tunnels” well intakes should 
several feet lower than dam spillway level, get clear “recharge 
water.” 

This benefit unpredictable without thorough study the geological 
mations and prediction quantity underground streamflows and their 
probable pattern. This can well sizable benefit and should given more 
consideration during planning. The feasibility and forth, 
for collecting and distributing underground waters would considered the 
cost studies. primary consideration would whether not the additional 
waters would seasonal (intermittent) continuous. Perhaps the studies 
Monroe Canyon (and others) may serve guidelines. 

Mr. Ferrell points out that anticipated, after channel stabilization, “that 
any debris that does enter the reservoir will fine material which can 
more readily sluiced.” Whether this material sluiced through debris dam 
excavated from debris basin, would seem that this finer material, hav- 
ing comefrom mountain-watershed topsoil and duff, would often contain 
able plant nutrients, and might economically salvaged and utilized nur- 
series, ranchers, gardeners, park departments topsoil can ef- 
fectively advertised, stockpiled, and distributed. This could some cases 
result income greater than the costs advertising, stockpiling, and distri- 
bution, which “net income” profit could applied against channel stabili- 
zation project costs. 

This may well local condition, that is, watersheds with ideal cover, 
humus, and forth, that will furnish high percentage desirable silt 
related the total eroded material. The cost excavating and transporting 
selected basin deposits may too expensive inasmuch removal costs may 
approach $3.00 per cubic yard. This not planning item but one 
considered the time basin sluicing debris removal. 
Mr. Ferrell quotes one policy item adopted about 1956 his District 
“the spillway height should not less than above the stream bed nor 
greater than ft.” not logical that some dams should higher than 
ft? These could the dams where groundwater conservation benefits would 
great enough influence the design and economics. perhaps they would 
dams located the better damsites, where favorable topography might indi- 
cate that dam higher than would more economical. recognized 
that other factors control the maximum desirable height these de- 
bris barriers, such safety, standard crib member dimensions, other de- 
sign model studies and stream profile studies. However, 
further study the benefits increasing deep percolation for water 


vation may show that dams barriers higher than may optimum 
certain locations. 
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ELECTRONIC COMPUTERS USED FOR HYDROLOGIC 


Discussion Alexandre Preissmann 


ALEXANDRE authors congratulated for hav- 
ing shown how important the use electronic calculating equipment can for 
very hydrologic problems. Calculations for which desk com- 
puters appeared hopelessly inadequate require comparatively little time 
and expenditure solved means electronic equipment. 

felt that the authors could have gone even further, for the use elec- 
tronic calculators can quite appreciably change the ways and methods ap- 
proaching hydrologic problems rivers, well setting new problems. 
hoped that the following examples may help present this idea more ex- 
plicitly. 

Before predicting the overallsummer flowfrom meteorological and physio- 
graphic data, multiple correlations have established. only desk com- 
puters are available, very small amount data, likely determine the 
overall summer flow (such estimate the snow cover), are introduced 
the regressions. However, electronic equipment available, one may 
tempted introduce large amounts data the forecast equations, and then 
pick out from all the possible regressions the one corresponding the high- 
est total correlation factor. Although such method would acceptable the 
available records covered extremely large number years, this require- 
ment unfortunately never met real life. The available records only ex- 
tend over very few years, and computed correlation factors are affected 
random errors that are quite appreciable. 

The highest computed correlation factor out hundred does correspond 
the best forecast equation ever likely correspond the highest true total 
correlation coefficient. 

Unless one chooses poor prediction method—or, extreme cases, one 
that actually absurd—one must possession criteria enabling rea- 
sonable choice made. For instance, one could ask oneself whether the 
formulas obtained for can reasonably applied other similar 
basins (that implies increased use calculating equipment) which case 
the considered formulas would become more reliable. alternative method 
would discard any formulas which cannot interpreted physically. 
Whatever the case under consideration, not true say that the use 
high-speed calculating machines supersedes the hydrologist’s experience. 
the contrary, gives its full meaning. 

Flood problems are approached analyzing the records and attempting 
use the observations find empirical relationships between the lag time 


November, 1959, Francis Swain and Herbert Riesbol. 


Engr. (Societe Grenobloise d’Etudes d’Applications Hydra- 
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between two points, the spreading out the flood, and certain data, such 
the seasonor the magnitude the flood. Where electronic equipment avail- 
able, this may not necessarily the most rational approach, for such equip- 
ment can used construct mathematical “model” unsteady flow the 
river, based the principles conservation mass and momentum. ad- 
justing certain data (coefficients frictional resistance), one can arrive 
rational understanding the observed effects. running such 
which would quite out the question with ordinary desk ef- 
fect man-made changes the course the river (such storage up- 
stream reservoirs, dykes, cutoffs, and forth) upon the flood regime can 
estimated within good degree approximation. 

would appear, therefore, that the available electronic calculating machines 
now may expected not only modify appreciably the details computation 
procedure, but also—and this much more give rise prob- 
lems new type requiring new methods solve them. 

particular, the use electronic equipment for systematic studies 
large number basins may well make possible, the long run, extend 
empirical rules originally established for individual basins any basins 
certain type, and thus gain clearer insight into the runoff phenomenon. 
this connection, felt that users electronic calculating equipment for 


hydrologic problems might benefit discussing difficulties and results 
periodic intervals. 
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HYDRAULIC DOWNPULL FORCES HIGH HEAD 


Discussions Robert Cox and Ellis Pickett, and Simmons, Jr. 


ROBERT COX,! ASCE and ELLIS PICKETT,2 ASCE.—The 
author has presented valuable andinteresting paper the vertical hydraulic 
forces face and cylinder gates. There definite need for more 
mationon this subject, especially the data can dimensionless 
parameters suitable for design purposes. 

This limited whichare normally locatedsome 
distance downstream from the sluice tunnel entrance. The control gate re- 
cess either open gate well bonnet-type recess. Fig. shows the 
features these two types recesses. 

Factors considered this discussion are: The effects gate clearances 
and venting the hydraulic force the top the gate;2) The effects gate 
bottom geometry and venting the hydraulic force the bottom the gate; 
The development dimensionless parameters for The effects 
gate bottom geometry and venting the stability the gate. 

The equation for the hoist load resulting from the hydraulic and gravity 
forces acting the gate shown Fig. 


which denotes the hydraulic and gravity tons, the dry weight 
gate tons, signifies cross-sectional area gate square feet, de- 
notes average unit-of-area downthrust top gate feet water, equals 
average unit-of-area upthrust bottom gate feet water, the 
specific weight water, 0.0312 tons per ft3. The friction load additional 
should considered actual design. The friction loadis dependent 
the pressure drop from the upstream the downstream side the gate, 
the face area the gate, and the rolling sliding coefficient the materials 
contact. 

The weight the gate not included the data presented here but re- 
view buoyancy and “wet weight” the gate interest. Fig. simplifies 
definition the involved. For the purpose discussing buoy 
ancy, gate may assumed rectangular parallelepiped with the verti- 
cal axis coincident with the direction gravity. the body completely in- 

November, 1959, Donald Cogate. 

Chf. Analysis Sect., Hydr. Analysis Branch, Army Engr. Waterways Experi- 
ment Sta., Vicksburg, Miss. 

Chf., Prototype Sect., Hydr. Analysis Branch, U.S. Army Engr. Waterways Experi- 
ment Sta., Vicksburg, Miss. 

Works and Spillway for Garrison Dam, Missouri River, North Dakota, Corps 


Engrs., Army, Waterways Experiment Sta., Tech. Memorandum 2-431, March, 
1956. 
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closed, the buoyant force still water equal the difference between the 
total pressure top (downthrust) and the total pressure the bottom (up- 
thrust). For such inclosed body, water pressure the upstream face has 
vertical component pressure. 

Some engineers use the expression, the “wet weight” gate. This 
simply the dry weight air minus the buoyant force. the body cellular 
lacks upstream skin plate, the wet weight that the complete- 
inclosed body shown Fig. The gate shown further considered have 
horizontal projections such gate seals. 

Tests and Measurements.—The Army Corps Engineer Waterways 
Experiment Station, has made number tests, both model and prototype, 
evaluate the vertical hydraulic forces acting control gates. Measurements 
gate oscillations resulting from fluctuations these forces have also been 
made. All test results presented this discussionare for fixed gate openings 
and all dimensions and forces are prototype values dimensionless. 

The total hydraulic force acting gate motion slightly greater 
less than that acting fixed gate depending upon whether the gate being 
raised lowered. moving gate accompanied lag changes the 
water depth pressure top the gate and the discharge under the gate. 
However, the rate travel gate normally less than two per minute 
and the effects gate motion the hydraulic load are believed small. 

The pressure top the gate, downthrust, can determined mea- 
surement the water-surface elevation the gate well measurement 
the pressure the gate bonnet. The pressure, upthrust, the bottom 
the gate can measured means piezometer taps electrical trans- 
ducers located the gate bottom. Electrical transducers also permit evalu- 
ation the pressure fluctuations which contribute gate oscillations. Gate 
oscillations can also measured accelerometers, strain gages and dyna- 
mometers installed the gate system. 

The average hydraulic load can measured models spring scales. 
Dynamometers and are suitable both model and prototype test- 
ing and permit evaluation fluctuations the hydraulic loadas well aver- 
age values. Measurements with all these instruments normally include the 
weight the gate and the effects friction damping. Model gates can de- 
that the friction load negligible. the force usually 
computed. However, can indirectly evaluated the hoist load, down- 
thrust, upthrust and weight the gate are known. models, the weight the 
gate can compensated for calibration the measuring equipment but 
would generally computed prototype tests. The prototype hoist load can 
computed from the cylinder pressures hydraulically operated gates. 

Downthrust.—The unit pressure ontop gate, downthrust, dependent 
the geometry and size the upstream and downstream clearances between 
the gate and the gate recess and the relative pressure differences across the 
gaps. The flow across the top the gate has some hydrodynamic effect, but 
this considered negligible. series tests were the Garrison 
Dam model determine the effect gate clearance onthe downthrust. these 
tests the downstream clearance area remained constant 6.3 while the 
shows the geometry and dimensions the clearances tested. Fig. shows the 
decrease indownthrust with eachdecrease upstream clearance. The curves 
figure reflect the effect changes gap geometry wellas gapareas. 
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Air vents are provided control relieve low pres- 
sures when back pressure from full conduit flow does not exist. air vent 
not provided, the reduced pressure behind the gate may cause lowering 
the water surface the gate well, reduction the pressure the bottom 
the gate, and decrease the stability the gate. 

The effect air vent the pressure differential across the gaps in- 
the curves Fig. The curves are the Ft. Randall Dam 
model tests which the downthrust was measured with the air vent down- 
stream the gate open and then closed. Pressures the conduit side the 
downstream gap were essentially atmospheric with the vent open. The pressure 
differential gap was greatly increased closing the air vent. The 
water-surface elevation the gate well lowered until the flows through the 


upstream and downstream gaps were again equalized. appreciable decrease 
the downthrust resulted. 


UPSTREAM 
CLEARANCE 
AREA 


0.0 FLOW TUNNEL 


ROOF 


SECTION THROUGH GATE WELL 


FIG, CLEARANCES TESTED 


The energy friction losses the flow through the gaps and, consequently, 
the downthrust dependent upon the gap surface roughness and flow Reynolds 
numbers well gap areas and geometry. Model results should trans- 
ferable directly prototype values the model upstream and downstream gaps 
have the same proportion relative eachother those the prototype. How- 
ever, the model and prototype flow quantities through the gaps not neces- 
sarily have the correct scale ratio. the designer has reliable data gap 
discharge coefficients, the gate well water-surface elevation and, consequently, 
the downthrust could computed for any gap geometry and pressure differen- 
tial. Gap areas and shapes vary with each structure designed and frequently 
with each gate opening. Therefore, series laboratory tests determine 
discharge loss number gap sizes and shapes would help- 


and Outlet Works, Fort Randall Dam, Missouri River, South Dakota,” 
Corps Engrs., U.S. Army, Waterways Experiment Sta., Tech. Report No. 2-528, Oct- 
ober, 1959. 
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ful the designer. The necessary tests could made large-scale, two- 
dimensional gap models. 

The gate well can sucked completely dry water with certain combi- 
nations upstream and downstream gap areas between the gate and the roof 
the conduit. then, the upthrust exceeds the weight the gate, the entire 
body the gate will vertically upward. Such phenomenon occurred 
one Corps Engineers project. Therefore, the gate recess should de- 
signed that the downthrust combined with the weight the gate will exceed 
the upthrust the bottom the gate. 

Another important design consideration the high velocity jet flowing 
through the downstream gap and the accompanying local low pressures which 
could cavitation damage. The first fewfeet the recessat the down- 


NOTE: DOWNSTREAM GAP 
AREA CONSTANT 


GATE OPENING FEET 


FIG, 4.—EFFECT UPSTREAM GAP AREA DOWNTHRUST 


stream gap sometimes steel-lined excessively low local pressures are 
anticipated. 

Upthrust.—The hydrodynamic effect water flowing past the bottom the 
gate substantial. reduction the bottom from the theoretical 
static head generally called “downpull” (Fig. Downpull may viewed 
ent upon the geometry the gate bottom. 

The Waterways Experiment Station made series tests number 
gate bottom designs part the preliminary design studies for Garrison 
and Ft. Randall Dams. Complete descriptions these tests are else- 
where. Seven the gate designs tested are shown Fig. Special care was 
taken reduce the gate friction load minimum the use steel roller 
bearings that contacted the gate slots all directions and the omission 
gate seals. The total vertical hydraulic load was measured independently 
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means spring scale and strain gages. The submerged weight the gate 
was compensated for during equipment calibration prior beginning each test. 
The downthrust was computed from the observed water-surface elevations 
the gate well. The upthrust was computed substracting the downthrust from 
the total hydraulic load. Fig. also shows the gate bottom geometry 
upthrust with the air vent open. Negative upthrust indicates less than at- 
mospheric pressures the bottom the gate. For type gate this pressure 
would amount approximately -28 water the prototype and probably 
cavitation would occur the gate bottom. Five the designs tested involved 
extension the downstream gate lip from ft. Each increase gate lip 
extension was accompanied increase total upthrust. Type the 
standard design adopted the Corps Engineers. structurally more 
desirable than most the other types and results favorable upthrust. Type 
hydraulically adequate but the parabolic gate bottom curve presents fabri- 
cation difficulties. This design was installed Denison Dam prior design 


DOWNTHRUST TONS 


OPENING FEET 


FIG, 


DOWNSTREAM VENTING DOWNTHRUST 


type 13. The curves shown Fig. supplement those presented the au- 
thor’s Figs. and 

presents upthrust curves for four the model gate designs (see Fig. 
with the downstream vent closed. Only type had positive upthrust all gate 
openings. should noted that upthrust values lower than about -40 tons for 
these model gates indicate prototype pressures lower than absolute zero. There- 
fore, the actual minimum upthrust prototypes these gates would about 
-40 tons. 

discussion upthrust would complete without mentioning determi- 
nation gate bottom pressures the flow net method. The technique for 
drawing flow nets have been discussed numerous publications. typical 
flow net for evaluation pressures the bottom gate Rouse. 


“Engineering Hydraulics,” Rouse, John Wiley and Sons, Inc., New York and 
Chapman and Hall, Ltd., London, 1950, 538. 
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FIG, 6.—EFFECT GATE BOTTOM GEOMETRY 
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Dimensionless Presentation.—The data far presented are dimensional 
and are not easily applicable design problems. The Waterways Experiment 
Station has prepared dimensionless design downthrust and upthrust 
for use the designoffices the Corps Engineers. Simplifications these 
charts are included Figs. and 


UPTHRUST TONS 


-80 


FIG, 7.—EFFECT GATE BOTTOM GEOMETRY 
UPTHRUST, AIR VENT CLOSED 


“Hydraulic Design Criteria,” Corps Engrs., Army, Waterways Experiment 
Sta. 
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FIG, 8.—DIMENSIONLESS PLOT EFFECT 
RECESS GEOMETRY DOWNTHRUST 
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FIG, 9.—DIMENSIONLESS PLOT EFFECT GATE 
BOTTOM GEOMETRY UPTHRUST 
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Fig. shows top the gate. The curves are 
basedon field gate-well water-surface elevations Ft. 
Randall Dam, South Dakota, and bonnet pressures the Pine Flat, dam 
California, and Norfork dam Arkansas. The observed pressures feet 
water were divided the velocity head the jet the vena contracta 
correlate different reservoir elevations. This velocity 
head was used because excludes all energy losses upstream from the gate. 
Details the clearances between the gate and gate recesses are also shown 


0.4 


OPENING FEET 


AIR VENT CLOSED 


OPENING FEET 


FIG, GATE BOTTOM GEOMETRY 
AND VENTING GATE OSCILLATIONS 


«Vibration and Pressure-Cell Tests, Flood-Control Intake Gates, Fort Randall Dam, 
Missouri River, South Dakota,” Corps Engrs., Army, Waterways Experiment 
Sta., Tech. Report No. 2-435, June, 1956. 

Pressure and Air Demand Tests Sluice, Pine Flat Dam, 
Kings River, California,” Corps Engrs., U.S. Army, Waterways Experiment Sta., 
Paper No. 2-75, February, 

“Slide Gate Tests, Norfork Dam, North Fork River, Arkansas,” Corps Engrs., 
Army, Waterways Experiment Sta., Tech. Memorandum No. 2-389, July, 1954. 
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this figure. The area used computing the downthrust should include 
the area the gate within the gate slots, the area between the conduit walls 
and the area the gate top seal. The large downthrust values for the Ft. Ran- 
dall gate result from the fact that the upstream clearance considerably great- 
than the downstream clearance. The smaller ratios the gap clearances 
for the Pine Flat and Norfork gates resulted appreciable reduction 
downthrust. 

Fig. shows similar dimensionless curves for upthrust the sloping bot- 
tom four 45-degree gate-bottom designs. The Ft. Randall gate has down- 
stream skin plate the gate proper, upstream skin plate the sloping 
gate bottom and downstream seals. The Pine Flat and Norfork gates have up- 
stream skin plates and downstream seals. 

The upthrust forces were computed from observed pressure data the 
sloping gate bottom. Pressure contours were drawn, integrated, and divided 
the projected gate area between the conduit walls determine the upthrust 
per unit area cross section. The velocity head the jet the vena 
contracta was used correlate data for different reservoir elevations. The 
curves indicate that the ratio conduit-width average gate-thickness 
factor the magnitude upthrust per unit area. The average gate thickness 
includes the gate bottom seal. 

Stability.—Vertical oscillations the gate caused pressure fluctuation 
the gate bottom are not usually consideredin hoist equipment design. These 
oscillations alternately decrease and increase the hoist-load. The friction- 
load should normally damp out minor gate oscillations. However, desir- 
able design gate bottom geometry for minimum vertical movement the 
gate. Fig. shows the maximum observed gate movement for type and 
gates with the air vent open and closed. The curves indicate the relative 
stability that will obtained with the 45° sloping gate bottom (type 13). 

agreement with the conclusions the au- 
thor. Vertical hydraulic forces gates based model tests are needed for 
economical design hoist capacity. Confirmation model test results 
prototype measurements always desirable. Considerable information 
available gate bottom pressures, but relatively little known about gate 
clearance criteria and the resulting downthrust. Measurement downthrust 
prototype structures relatively simple and hoped that such measure- 
ments will made whenever possible. For design purposes, necessary 
compute the height the water-surface the gate wellor the pressure head 
the bonnet. believed that discharge coefficients could determined 
from model investigation gap areas and geometry, and that large-scale 
tests should minimize any scale effects. Measurements the flow through the 
gaps together with pressure differentials would permit evaluation the re- 
quired discharge loss coefficients. Gaps for existing prototype structures 
should included any general model investigation. The results sucha 
series tests would provide valuable data for designers. 


SIMMONS, ASCE.—The author clearly points out number 
difficulties encountered with the direct weighing methodof determining down- 
pull tests high head gates that are fully submerged. The usual use 
pressure tank for these high heads place cumbersome head box 


Hydr. Engrg. Research, Bur. Reclamation, Denver, Colo. 
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introduces friction the point where the gate lifting stem emerges from the 
tank. this respect, the head tests differ from the simpler ones often 
conducted for spillway gates, land forth. the latter case, downpull mea- 
surements direct reading scalesor similar possible with- 
out any such frictional interference the stem itself. 
Strain gages similar the ones mentioned the prototype studies for the 
Shasta dam may placed the lifting stems models within the pressure 
tank. They, thus, offer partial solution the problem that they measure 
the the eliminate any stem frictionfrom measure- 
ments. They not eliminate another major problem posed the sliding 
rolling friction between the leaf itself and its supporting tracks guides. 
Several excellent methods for reducing eliminating this leaf frame 
friction have been However, all these methods seem have 
the effect releasing the seals from the seating surface, thereby permitting 
leakage. Leakage past the top seal can appreciably alter the measurement 
heads acting the top the leaf high head gates because will change the 
depth water within the bonnet shaft and, hence, over the leaf top. Leakage 
past either the top the side seals the gate can influence the flow pattern 
filled tunnels downstream from the leaf and the leaf bottom. Unfortunate- 
ly, the effects leakage are relatively indeterminant factors, but they are, 
nevertheless, real and must considered. 
Downpull determinations the pressure-area methodalso have limitations. 
Adequate coverage representative areas the leaf bottom many piezome- 
ters mandatory but not always easily obtained. Frequently, sizeable per- 
centage the gate leaf situated within the gate slots. This particularly 
true fixed-wheel and roller-train gates that are built from standard 
structural rolled beams and plates and require very wide slots. The pres- 
sures acting the portions the leaf within the slots are markedly different 
from those the main flow passage. Furthermore, the pressure gradient from 
point just inside the slot point deeply inside may steep and unpre- 
dictable pattern. Only extensive piezometric coverage will suffice produce 
reasonably accurate downpull results. further complicate the situation, 
appurtenances like box enclosures around trains flow 
below the fixed wheels must consideredand included. the direct measure- 
ment method, the effects these variables not require special treatment, 
other than being certain that the gate truly representedin the model, because 
their effects are automatically included the load measurement. 
Difficulty may also encountered deciding what the pressure acting 
downward the top the gate should and where will act. The problem 
simple the case fully enclosed, downstream seal leaf because the 
pressures the bonnet can measured and referred the elevation the 
gate top. But the case where the leaf built from structural members 
and has open upstream face, may impossible say where the down- 
ward force may considered act and, hence, what the effective force will 
be. addition there may appreciable downward drag force due the action 
high velocity water the exposed horizontal beams. Thus, the net value 


«Vertical Hydraulic Forces Water Controlling Gates,” Bryce and 
Foulds, Proceedings Eighth Congress, IAHR, August, 1959. 

“Investigations the Total Pulsating Hydrodynamic Load Acting Bottom Outlet 
Sliding Gates and Its Scale Modelling,” Abelen, Proceedings Eighth Congress, 
IAHR, August, 1959. 
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found subtracting the upward acting leaf bottom pressures from the down- 
ward acting leaf “top” loads and drag forces may not readily obtained. 
the direct measurement method these problems need not considered because 
their effects are included the stem load measurement. 

It, thus, appears that there nosingle best method for determining, through 
model studies, the hydraulic downpull forces acting all high head control 
gates. Each determination will require careful evaluation the principal fac- 
tors contributing downpull and the test arrangements most suitable for 
that structure. some cases, may desirable use both testing methods. 
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REFERENCE PRESSURE-LEAF BOTTOM PRESSURE 
VELOCITY HEAD BENEATH LEAF 


FIG, 


There are other examples!3 downpull studies based the pressure-area 
method described the author. These studies were made using air the 
flowing fluidinstead water. This practice which canresult 
able economies time, money, and testing equipment. One the studies con- 
cerned high-pressure, downstream seal, slide gate. The leaf bottom hada 
45° sloping surface that was followed short, horizontal section that seated 
the floor (Fig. 1). The gate was the type used the United States Bureau 


“Air Model Studies Hydraulic Downpull Large Gates,” Simmons, 
Proceedings, ASCE, Vol. 85, No. January, 1959. 
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Reclamation, Dept. Interior, Palisades Dam, Idaho for free-discharge 
regulation under heads much 237 ft. 

The data obtainedin the air tests for the condition sub- 
merged discharge because the airflow was discharging into limitless reser 
voir like fluid, the atmosphere. However, was stated that the pressure 
distribution upon the sloped portion the leaf bottom would the same for 
either free discharge submerged operation because the flow pattern this 
area would remain the same. The pressures themselves would affected 
any submergence the structure, but the difference pressure from one 
station another this portion the leaf would for either 
type operation given rate flow. 

the interim the writer’s and the publi- 
cation Mr. Colgate’s paper, the writer has studies with 
components the air model. these later tests, water was used the flow- 
ing fluid that both free discharge and submerged conditions could repre- 
sented. Considerable strengthening and waterproofing were required the 
model enable withstand the greater loads imposed the water and 
avoid excessive warpage. further limit distortion the model, the tests 
were conducted short possible. new leaf was used for the 
water tests place the original one, which had become damaged and un- 
serviceable. 

The test results were plotted two different forms. the first, the mea- 
sured pressures acting the leaf bottom, divided the velocity head the 
vena contracta the flowing water, were plotted against gate opening (Fig. 2). 
readily apparent that the free discharge and submerged discharge values 
differ greatly, and that the submerged values for these tests are lower than 
the free discharge ones. The previous air model data arealso plottedand show 
good agreement with the submerged tests using water. 

parameter that better describes the pressure distribution within flow 
the pressure any point the system;for instance, the leaf bottom slope, 
and the velocity head the flow passing beneath the gate leaf. The test 
data, when expressed this form, take into account the effect tailwater be- 
cause the upstream reference head, which also directly affected the tail- 
water, used the datum. The free discharge and submerged discharge data 
obtained with water tests plot together (Fig. 3). The submerged discharge data 
obtained with air tests show similar pattern, but with slightly lower values. 
Thus, the pressure distributions the sloped portion the leaf are shown 
the same for free discharge and submerged operation, and the reasonable 
agreement the air and water evidence that air models can used 
effectively for certain downpull studies. 
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DISCHARGE FORMULA FOR STRAIGHT ALLUVIAL 


Corrections Discussion Blench 
Discussions Lean, Lucien Brush, Jr., Don Culbertson 
and Paul Jordan, and Bruce Colby 


1960 Journal the Hydraulics Division: Page 131. Complete the footnotes 
bottom page as: 


“Regime Behaviour Canals and Rivers,” Blench, Butterworths Sci- 
entific Publications, London and Toronto, 1957. 


“The Fourth Root n-f Diagram,” Blench, Proceedings, ASCE, Vol. 86, 
No. HY1, January, 1960. 


proposed formula rests data drawn mainly from 
laboratory flume work, natural streams which are small enough rank 
flumes. exception the data Lane and for depths ftand 
pebbles with mean size between 0.8 in. and3.2 in. With this exception, the 
depths are seldom greater than ft. Thus, for rippled beds the relative rough- 
ness will large, and the effect small relative roughness which might 
quite important for deep natural streams not taken into account. So, also, 
the large-scale sand waves which can occur deep and which add ap- 
preciably the roughness are ignored. The general applicability the index 
and the coefficient is, therefore, doubt. 

The formula for plane beds composed sand and silt (size 0.06-0.6 mm) 
little direct interest because with long continued action such materials 
move ripples except certain which are seldom found under 
natural conditions. The evidence that for rippled bed the velocity for given 
depth not proportional not presented. For the index beless than 


(Fig. 10), the ripple roughness must increase (increase height) the ve- 


locity increases. Flume studies the United Kingdom’s Hydraulic Research 
Station 0.17 sandat well below the critical shear which ripples 
not support this. 

should, perhaps, stated that the agreement the calculated and ob- 
served velocities given Fig. does not prove the formula but verifies the 
fit the constants the formula the data. 


Prin. Scientific Officer, Hydr. Research Sta., Wallingford, Berkshire, 
Lane and Carlson, Transactions, U., Vol. 35, 1954, 453. 
“Systematic Changes the Beds Alluvial Rivers,” Carey and 
Keller, Proceedings, ASCE, Vol. 83, No. August, 1957. 
4R. Bagnold, Philosophical Transactions, Royal Soc., London, 1956, 249, 235. 
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The writer the opinion that field study the ripple and bed configu- 
ration from which extract equivalent Nikuradse roughness more re- 
liable than the use power formula the type advocated. this connection 


table ripple sizes versus equivalent uniform roughness would invalu- 
able. 


LUCIEN BRUSH, attempt establish discharge formula for 
alluvial channels commended because the urgent need equation 
for practical design. Quite justifiably, the authors have used empirical 
analysis due the difficulties they encountered attempting analytical 
approach. However, rational empiricism, means dimensional analyses, 
has sound analytical basis leading the development dimensionally ho- 
mogeneous expressions. often tempting sacrifice dimensional homogene- 
ity, the authors have done, for the purpose simplifying some unwieldy 
expressions, but doing, the resulting equation usually not generally 
applicable. Thus, the proposed equation, Eq. 31, may some use for la- 
boratory flumes, but cannot extended prototype alluvial channels, even 
“first approximation,” suggested the authors, without extensive field 
verification. 

The solution Eq. for the mean velocity requires that the particle size, 
slope, temperature, and the hydraulic radius, with respect tothe bed, known. 
This somewhat misleading because the hydraulic radius with respect 
the bed, not length parameter which can measured physically, but 
actually conceptual length, the determination whichrequires the knowledge 
the slope, width, depth, velocity (as the authors note), shape the channel, 
and hydraulic roughness the side walls. Thus, truly known, there 
only one possible solution for the mean velocity and this solutionis completely 
independent the empirical curves presented the authors. other words, 
obtained applying wall the hydraulic radius. the 
standard wall correction valid, only one solution exists for the velocity. 
plot computed versus measured velocity would have noscatter compared 
Figs. and shown the authors. show that Rp, width, depth, 
slope, temperature, and the shape the channel are known for particular 
flume, only one solution for the velocity exists, necessary only apply 
the wall correction outlined Vanoni and Brooks. 


coefficients for the bed and entire channel, respectively. 
For smooth walls and rectangular cross section, 


Research Engr., Iowa Inst. Hydr. Research, and Asst. Prof., Dept. Mechanics 
and Hydr., State Univ. City, Iowa. 


“Formula for theTransportation Bed-Load Movement,” Einstein, Trans- 
actions, ASCE, Vol. 107, 1942. 
“The Importance Considering Side-Wall Friction Bed-Load Investigation,” 
Johnson, Civil Engineering, Vol. 12, No. 1942, pp. 329-331. 
“Laboratory Studies Roughness and Suspended Load Alluvial Streams,” 
Vanoni, Report No. E-68, California Inst. Tech., 1957. 
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which pp, and are the wetted perimeter the bed, wall, and entire 
channel, respectively, and the resistance coefficient for the walls. 
Compute equals inwhich Rand are the Reynolds numbers 

the channel and walls, respectively. 

From smooth-pipe diagram find fy. 

Substitute values for fp, and the equation step 

Repeat this process until step balances. The resulting velocity then 
the unique solution for the given conditions. 


For other cross sectional shapes and different types wall roughness, the 
procedure more complicated but will also yield only one solutionfor the ve- 
locity. requiring that and known, the authors demand, the 
solution for the velocity predetermined and cannot vary unless the wall- 
correction procedure assumed incorrect. 

This result might have been eliminated had the authors chosen either the 
depth the hydraulic radius the final equation (Eq. 31) instead Rp. These 
parameters can measured assumed for different boundary shapes with 
knowledge the mean velocity. course, leaving out the wall correc- 
tion equally difficult problem introduced data from different flumes 
are compared. This would unavoidable for the particular approach 
used. 

Dimensional Analysis.—Had the authors considered all the assumptions 
which they eventually make pertaining grain shape, sediment uniformity, 
density variation and supercritical flow during the dimensional analysis, ap- 
pears that the computational procedure could have been greatly simplified. 
More than being useful tool for efficiently organizing series variables, 
the procedure used performing dimensional analysis forces the investi- 
gator examine closely the pertinent variables involved the problem. Only 
one the variables may dependent, while the others must independent. 
more than one dependent variable included, superfluous terms will 
generated which, used graphing the results, may yield correlations which 
have bearing the particular problem being considered. Obviously, the 
exclusion important independent variable will have more deleterious 
the analysis and the resulting graphical representations. The authors 
have chosen include terms which are superfluous and have omitted some 
which are importance. 


The variables which the authors feel are important are arranged the 
functional relationship, 


D= 9 (a, S, g, d, 0, n) 


The inclusion particularly requires explanation. Sigmao, the 
standard deviation the particle size distribution may effect the dune height 
and doing affect the depth flow. However, there seems less jus- 
tification including the particle-shape factor independent variable 
determining the depth flow. Compared the other variables, the shape fac- 
tor must be, best, second third-order effect. The authors admit this, 
course, because they state that: order improve the accuracy the 
method necessary determine the effect particle shape the mean 
velocity.” However, the inclusion allows the authors substitute fall ve- 
locity their list variables which already contains particle size. Without 
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only one the other the variables, particle size and fall velocity, may 
included because their interdependence (see page 74). Ultimately, both and 
are considered relatively unimportant the authors and are excluded 
from the analysis the determination the discharge coefficient C,. With 
the end result being the exclusion these variables, difficult under- 
stand why either the particle Reynolds number the relative-roughness pa- 
rameter should not excluded from the graphical representations (despite 
the authors’ comments the contrary), and why the original functional re- 
lationship should not 


For the assumption uniform steady flow, which was made the authors, 
the repeating variables, the functional relationship for the dimensionless 
parameters would 


the authors have considered only subcritical flow, the third term—the 
Froude discarded. However, any important inertial effects 
are retained the fifth term. Eliminating the third term, the functional re- 
lationship then reduces 


Although complex relationship still found exist, the resulting equations 
4(a) much easier work with thanthe extremely complicated fac- 
tor used the authors. However, for the assumptions ultimately made the 
authors, Eqs. 4(a) (b)are equal inefficiency the equations which they use. 
Eq. 4(a) suggests that the approach shown Fig. might yield fruitful results 
the fourth term included. Eq. 4(b), the other hand, suggests that Fig. 
might the starting point for equally efficient grouping the basic pa- 
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rameters provided the additional terms given Eq. 4(b) are used. The au- 
thors, course, have used Fig. for the purpose classifying the bed re- 
gimes. 

the authors incorporate bed regime their analysis, appropriate 
consider briefly another aspect the problem. functional relationship 


for the dune height, may obtained using approach similar that 
for arriving Eq. 4(a) 


Eq. 5(a) combined with Eq. 4(a), after multiplying through d/D, the 
following result may obtained 


Eq. contains two dependent variables each which function the 
same three independent variables. Although Eq. can have and effect 
significance, may advantageous group the data the form Eq. for 
plotting purposes, done plotting typical lift versus drag-coefficient 


diagram. For Reynolds numbers and rather conservative range d/D, 
Laursen relationship 


which both consistent with the rough-pipe analogy considered and Eq. 
The similarity Eq. Fig. with anadded d/h term worthy note. 

Finally, appears that the general approach used attempt- 
ing obtain equation for predicting the velocity flow alluvial channel 
should pursued farther. Conceptual lengths such Rp, course, must 
eliminated. Bed classification schemes must eventually replaced quanti- 
tative descriptions the bed geometry. Dimensional homogeneity should 
retained throughout the analysis and resulting equations. However, when these 
details are worked out, conceivable that useful equationfor the prediction 


velocity alluvial channel can obtained using the general outline 
given the authors. 


DON ASCE and PAUL pro- 
posed discharge although based largely laboratory flume data repre- 
senting shallow depths and relatively steep slopes, combines directly in- 
directly many the variables influencing resistance flow alluvial chan- 
nels. The authors are congratulated for developing discharge formula 


“The Total Sediment Load Streams,” Laursen, Proceedings, ASCE, Vol. 
84, No. 1958. 
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based empirical correlations the important flow, and sediment parame- 
ters and the bed configuration. The formula should interest laboratory 
workers and may, with additional research, have practical application chan- 
nel design and some special problems encountered natural streams. 

The authors point out that the bed-configuration classification shown their 
Fig. not representative for large natural rivers. They also point out the 
probable necessity modifying the discharge formula suit field conditions. 
Therefore, this limited comparing the bed configuration se- 
lected rivers with that shown the authors’ Fig. and comparing the field 
measurements velocity with those computed from the discharge formula. 

Representative examples the available river data are shown Fig. 
which V/V, the third variable. Although these data represent wide range 
depth, slope, mean velocity, water temperature, and sediment discharge, 
they represent narrow range median particle size bed material. All 
the rivers considered are reasonably straight the measuring sections; bank 
friction probably insignificant because the depth-width ratios are low. Be- 
cause the bed configuration usually difficult observe flow more than 
about deep and because accurate measurements bed configuration for the 
Loup River only are available the writers, the ratios V/V, are 
used indicate the bed configuration. Low ratios V/V, (about for 
the Cedar, Middle Loup, Niobrara, and Little Blue Rivers and about 
for the Mississippi River St. Louis) indicate dune bed. High ratios V/V, 
(about for the Middle Loup River and about for the Mississippi 
River) indicate transition flat bed. Intermediate ratios V/V, indicate 
dunes and transition different parts the cross section and are not plotted 
Fig. apparent from Fig. that the parameters and are 
inadequate estimate the bed configuration streams such the Middle 
Loup River, however, they may adequate for laboratory flumes and for 
some rivers. For median bed-material sizes 0.2 0.5 the river 
data indicate that the average curve the authors’ Fig. for transition bed 
too low. 

According the authors’ Fig. given shear velocity will produce the 
same bed configuration for fairly wide range median particle sizes, ex- 
cept the region near the beginning sediment motionand ripple formation. 
For example, constant temperature 65° and median particle size 
shear velocity about 0.18 fps required maintain transition 
bed. the particle size 0.80 mm, the required shear velocity still about 
0.18 fps. 

Data from measurements confirm the authors’ Fig. Al- 
though these data generally are outside the range values the laboratory 
data, they are good agreement with the related lines constant This 
agreement was further apparent the authors’ Fig. Agreement was good 
even when the bed did not have the appearance dune bed and the roughness, 
indicated Manning’s Chezy’s was much less than common for 
dune beds. The excellent agreement may due partly thesimiliarity the 
parameters V,d/v and The small deviation both laboratory and field 
data from the curve Fig. indicate that Eq. should give very good re- 
sults when obtained from Eq. 32. Eq. contains 
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which for dune bed equivalent A-5-9. Note that error can re- 
sult about 35% error Thus, the computation from Eq. not 
only tedious but can result large errors. 

Velocities computed with the coefficients and exponents for dune bed are 
generally within about 20% the measured velocities, whereas velocities com- 
puted with the coefficients and exponents for flat bed are nearly always much 
higher than the measured velocities. Table shows the pertinent data for two 
rivers widely different measured velocities. The bed configuration for the 
Middle Loup River near Dunning, Nebr., was determined with ultrasonic 
sounder, and the bed configuration for the Mississippi River St. Louis, Mo., 
was assumed the basis 

The data Table are typical many field measurements; computed ve- 
locities for flow over flat bed are unreasonably high. Also, velocities com- 
puted with Eq. may unrealistic for alluvial rivers having small slopes, 


TABLE 1.—DATA FROM FIELD MEASUREMENTS 


Middle Loup River Mississippi River 
Dunning, Nebr. St. Louis, Mo. 
June 26, 1959 Dec.6, 1959 Feb. 21, 1956 1952 
Bed configuration Flat Dune Flat 
Median particle size 
bed material. mm. 0.31 0.29 0.31 0.48 
Hydraulic radius foot 2.15 1.52 16.5 44.7 
Shear velocity feet per second 
8.2 16.4 10.5 19.3 
Measured velocity feet per 
second 2.20 3.95 2.20 7.69 
Computed velocity: 
Flat bed feet per second 6.05 17.2 
Transition (ripples bars) 
feet per second 2.90 9.54 


Dunes feet per second 1.86 2.13 


Average amplitude 0.72 ft, average length 12.5 ft. 


shallow depths, and bed material coarse sand. For example, the computed 
velocity flow over dune bed higher than the computed velocity over flat 
bed whenthe slope, depth, and median particle size bed material the Little 
Blue River are 0.0005, 3.0 ft, and 1.00 mm, respectively. 

summary, comparison the observed bed configuration selected 
rivers with the configuration alluvial beds, according 
the authors’ criteria, suggests that the authors’ Fig. inaccurate the 
ranges dune bed and transition. Velocities computed with the authors’ dis- 
charge formula were much higher than the measured velocity for flow over 
flat bed but agreed reasonably well for flow over dune bed. 


BRUCE proposed discharge formula, which here as- 
sumed mean the authors’ Eq. and Figs. 10, and 11, shows clearly 
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the wide range resistance flow for different bed configurations. im- 
plies, the omission sediment concentration from the formula, that any 
major effect sediment concentration mean velocity correlates with bed 
configuration, fall velocity the median particle, The authors are 
commended for bringing out these points. 

Presumably the discharge formula intended for use natural streams 
without extensive modification, although was derived from flume experiments 
and accuracy was not claimed for the formula when applied outside the labo- 
ratory. The purpose this discussion explain some the difficulties 
using the formula for stream outside the laboratory. 

One disadvantage the formulais that was developed for shallower depths 
and higher slopes than those most natural streams. Hence, the extension 
the formulato deeper flows and flatter slopes may lead toinconsistent results. 
For example, the computed mean velocity slightly higher over dune bed 
than over flat plane bed for median diameter mm, hydraulic radius 
ft, and slope 0.0006. (Fig. will indicated later, not accurate 
enough show which bed configuration would actually exist for these con- 
ditions.) 

Other disadvantages the proposed discharge formula become apparent 
when used compute the depth-discharge relationship for Pigeon Roost 
Creek near Byhalia, Miss. The median diameter the bed material about 
0.4 mm, and the surface slope steady flow about 0.0011. Near the gaging 
station, the flow dredged channel about wide andstraight for many 
hundred feet both upstream and downstream. water temperature 16°C was 
assumed for the computations. Discharges per footof stream based 
directly Figs. 10, and and are plotted Fig. this discussion 
solid line that has horizontal jogs changes from one bed configuration 
another. 

Streamflow measurements Pigeon Roost Creek, near Byhalia, define two 
curves stage-discharge relationship. One the curves for flow over 
dune bed and one for flow over the field equivalent flat bed. The curves 
are expressed Fig. depth-discharge relationships per foot stream 
width and are applicable for steady flow. Manning’s equation with n’s 0.032 
canbe used reproduce rather closely the curves for the low-velocity 
curve and the high-velocity curve, respectively. The curve for high-velocity 
flow has break even depths when antidunes sometimes occur along the 
middle the channel. Some streamflow measurements, made when the dis- 
charge was rapidly changing and depths were within the approximate range 
ft, plotted between the two curves. 

Several conclusions can drawn comparing discharges that are com- 
puted from the proposed discharge formula with discharges that are based 
streamflow measurements. First, the agreement between computed and mea- 
sured flows over dune bed very good. However, Fig. incorrectly indi- 
cates that dunes should not exist depths greater than 0.9 for the median 
particle size, slope, and water temperature that were used, whereas dunes 
nearly always exist depths and may exist depths greater than 
ft. 

Second, over flat bed the discharges computed from the formula are ap- 
preciably higher than the measured discharges. Perhaps flat bed Pigeon 
Roost Creek less regular than flat bed laboratory flume. Also, bank 
friction may have appreciable effect depths greater than ft, but its ef- 
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conclusion that the formula, generally, may give somewhat too high dis- 
charges for flow over flat bed natural stream. 

Third, Fig. gives too low shear velocity for the change from dune bed 
bed. The other changes bed configuration are flows low that 
neither the changes nor the velocities have muchpractical significance many 
natural streams. 

Incidentally, Fig. does not indicate the limiting shear velocity for the 
change from dunes bars direct change from dunes flat bed has 
been assumed. Bars and other bed configurations intermediate between dunes 
and flat bedare relatively insignificant for Pigeon Roost Creeknear Byhalia, 
partly for reasons that are associated with length upstream channel. Con- 
sider reasonably uniform alluvial channel several miles long whose upper 
end the flow being increased slowly. When the flow becomes slightly too 
great transported over dune bed, the velocity increases considerably 
even though the discharge increases only little. Hence, the water the high- 
velocity tends overtake the slower water ahead it. the 
rising side the hydrograph flow becomes increasingly steep the flow 
moves downstream. 

After the water has traveled few miles, the flow the changefrom dune 
bed flatter bed longer increases gradually but may double within few 
minutes. The rapid increase relatively high discharge causes abrupt 
change from dune configuration flat bed. the recession side the hydro- 
graph, the flow some point upstream may decreasing gradually when the 
flat bed begins change bars dunes. The reduction velocity causes 
the channel drain less rapidly. the flow progresses downstream, the re- 
cession side the hydrograph flow becomes increasingly steep the time 
flat bed rougher bed. This faster decrease flow makes 
the change from the flat bed dune bed more rapid than might have been 
upstream where the decrease more gradual. other words, gradu- 
ally changing rates discharge may unlikely, not impossible, the change 
from dune bed flat bedor from flat bed dune bed longand rea- 
sonably uniform alluvial channel exists upstream. 

Fourth, comparison the computed depth-discharge relationships with 
those that are defined streamflow measurements indicates (Fig. that the 
shape the transitions from one bed configuration another are unrealistic. 
the transitions from plane bed rippled bed and from rippled bed dune 
bed, the discharge formula and the authors’ Fig. indicate that slight in- 
crease depth will reduce the mean velocity. (The computed reduction ve- 
locity the change from plane bedto rippled bedis thought too small.) 
such reductionin velocity actually was appreciable amount, 
the bed configuration probably would return the form that has lower re- 
sistance flow. 

The shape the transition computed discharges from dune bed 
flat bed is, also, incorrectly shown single horizontal line. According 
the computed discharges Fig. steady flows 1.4 3.4 cfs per cannot 
occur. That is, both the velocity and the computed discharge are discontinu- 
ous the change from dune bed flat bed. Measured velocities are, in- 
deed, discontinuous for steady flows the change from dune bed flat 
bed, but measured discharges for steady flows have appreciable discontinuity 
because the depth flow can readily adjust the bed configuration and the 
discharge. Thus, the depth-discharge curves defined streamflow measure- 
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ments overlap depth flow. (See Fig. 2). Within the range overlap they 
supply field verification the Brooks that the velocity cannot 
“expressedas single-valued functionof the bed shear stress, any com- 
bination depth and slope, bed hydraulic radius and slope.” this range 
about depth Pigeon Roost Creek near Byhalia, neither depth nor 
This fundamental weakness any relationship suchas that the authors’ 
Fig. that makes bed configuration function flow parameter except 
shear velocity. 

Overlapping the depth-discharge curves for and for flat bed 
somewhat unusual natural streams. requires that the bed almost 
entirely covered dunes some times and flat bed the same depth 
flow other times. Such consistently uniform bed configuration one type 
theother inthe depthfor whicheither dune bed flat bed might 
exist seems possible only the lateral distributions depth, velocity, 
and bed-material size are fairly uniform. 

Most streams have irregular lateral distributions depth, velocity, and 
bed-material size and hence part the bed one type configuration 
while another part has another type. Also, part the bed may have some 
intermediate type The alluvial streams that have 
cut discontinuity the depth-discharge relationship generally have range 
depth which the resistance flow varies widely and somewhat uncertainly. 
For those streams having reasonably average transition bed roughness 
from dune configuration flat bed, the method Einsteinfor comput- 
ing mean velocity probably gives reasonably accurate discharges. 

Even for Pigeon Roose Creek near Byhalia, the discharges computed the 
Einstein method (Fig. are somewhat more accurate than those computed 
the authors’ discharge formula. The depth which the transition between the 
two types bed configuration occurs also better determined from the Ein- 
stein method. However, for depths 1.5 3.0 the discharges computed 
with the Einstein method are considerably different than those that are deter- 
mined streamflow measurements. 

compared the flow Pigeon Roost Creek near Byhalia, the authors’ 
discharge formula gives computed discharges that agree well with measured 
flows over dune bed but are too high for flow over flat bed. Also, the au- 
thors’ Fig. inadequate measure the range depths which tran- 
sitions occur between dune beds and flat beds vice versa. 
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«Mechanics Streams with Movable Beds Fine Sand,” Norman Brooks, 
Transactions, ASCE, Vol. 123, 1958, 547. 

“The Bed-Load Function for Sediment Transportation Open Channel Flows,” 
Einstein, Tech. Bulletin No. 1026, Dept. Agric., Washington, C., Sep- 
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EFFECT AQUIFER TURBULENCE WELL 


wish comment only one aspect this paper, namely 
the view implicit that “in laminar flow, the energy losses are directly 
proportional tothe velocity,” and that turbulence term embracing all flows 
“under conditions that relate the energy loss velocity exponent that great- 
than one.” 

fluid mechanics, the terms “laminar” and “turbulent” have quite definite 
meanings, feel that must not ignore these meanings and set new 
and quite different ones. Although some profound questions about turbulence 
remain very much unsettled2 can quite definitely assert that flows may 
laminar, the sense being ordered flows? that obey the Navier-Stokes 
equation, without there being linear relation between flow rate and energy 
loss. the writer shall outline here, there very strong evidence that flows 
media inthe Reynolds number range 100 over (these num- 
bers how Reynolds number defined) belong this lami- 
nar, but non-linear, regime. The confusion this point seems have arisen 
from false analogies between flow tubes and flow porous media. Unfortu- 
nately, these misleading ideas “turbulence” porous materials appear 
the text books flow porous media. 


The steady state form the Navier-Stokes equation for incompressible 
fluid 


the (microscopic) vector velocity, and the potential defined 


which the pressure, the fluid density and the potential the 
nal forces. (Note that P/dand Qare, ina gravitational field, entirely 


analogous the velocity, pressure and gravitational components “head” 
hydraulic engineering.) 

Textbooks the derivation Eq. from simple considerations about 
the mechanics fluid motion. The forces and accelerations small 
element the fluidare relatedon the supposition (rigorous for Newtonian fluids, 
such water)that the shearing stress sustained the fluid proportional 


November, 1959, Joe Mogg. 

Research Officer, Div. Plant Industry, Canberra, Australia. 

For example, with stream surfaces constant time. “International Symposium 
Atmospheric Turbulence the Boundary Layer,” Geophys. Res. Papers, 19, pp. 479-491, 
1952. 
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the velocity gradient normal tothe shear. The (constant) ratio shear 
ing stress tovelocity gradient the dynamic viscosity the fluid, which 
equal da. 

Navier obtained early 1822, basis being improved 
Saint-Venant 1843. Curiously enough, this general analysis the motion 
viscous fluid preceded the much less general, but better known, studies 
Poiseuille and Darcy. does not seem been recognized the time that 
their results represented simple cases 

The terms the left hand side Eq. are known the inertia terms. 
These terms are identically zero for certain special motions including axial 
flow straight tubes. this case Eq. reduces 


follows that, for suchcases rate proportional tothe ener- 
“laminar” give way disorderly unsteady “turbulent” flow. well 
this occurs for straight tubes circular cross section Reynolds 
numbers the order 2,000. 

The difference character between flow tubes and porous media is, 
mathematical terms, that porous media the inertia terms are not identi- 
cally zero, and that the linear regime porous media depends the approxi- 
mation that the inertia terms may neglected. fact, can that 
this approximation can expected fail when the Reynolds number unity, 
perhaps little more (depending again the definition used). This failure 
the linear (Darcy) regime porous media flows Reynolds number 
range, from just over unity times this, well known experiment- 
The fact that the non-linearity arises where predicted the 
Navier-Stokes equation very cogent evidence that has its origins the 
non-vanishing the inertia terms and not the onset turbulence. (Of course, 
the flow rate were continuously increased, point wouldultimately reached 
where genuine turbulent regime had been established.) 

Physically, this difference character between laminar flow tubes and 
parcel fluid moving through each system. the tube the parcel possesses 
constant linear velocity throughout, while the medium subject 
which continuously change both the magnitude and direction its velocity. 

The Swedish engineer Lindquist8 recognized the central role that the in- 
ertia terms play producing the non-linear laminar regime 1932. Muskat 
this interpretation because diffusion marker solution increases 


“Modern Developments Fluid Mechanics,” byS. Goldstein, Clarendon Press, Ox- 
ford, Vol. 1938, 319. 

“Hydrodynamics,” Lamb, 6th Ed. Cambridge University Press, 599, 1952. 

“The Flow Homogeneous Fluids through Porous Media,” Muskat, McGraw 
Hill Book Co., Inc., New York, 1937, 56. 

Physics Flow through Porous Media,” Scheidegger, Univ. Toronto 
Press, 1957, 

the Flow Water through Porous Soil,” Lindquist, ler Congres des 
Grands Barrages, Stockholm, 1933, pp. 91-101. 

Flow Homogeneous Fluid through Porous Media,” Muskat, McGraw 
Hill Book Co., Inc., New York, 1937, 62. 
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with flow the phenomenon “hydrodynamic dispersal”10 
understood and occurs the absence turbulence. 

Finally, note that recent mathematical work Tamada and 
the drag cylinders and other obstacles gives very graphic pic- 
ture the way which the non-linear regime develops the Reynolds num- 
ber increases and the inertia terms grow importance. 

These matters have been treated further Irmay discusses!4 
the non-linear laminar regime different, though closely connected way. 

Summarizing, then, this note advocates (andadvances reasons for) the recog- 
nition three regimes flow porous media. These are, order in- 
creasing Reynolds number: (a) Linear laminar (or Darcy) regime, (b) Non- 
linear laminar regime, and (c) Turbulent regime. The paper under discussion, 
andother references the literature, misleadingly label item (b) “turbulent.” 


Due diffusion across stream surfaces and disruption the stream system 
each stagnation point porous medium. 

“Dispersion Moving Salt-Water Boundary Advancing Through Saturated Sand,” 
Day, TAGU, 37, 1956, pp. 595-601. 

Steady Two-Dimensional Flow Viscous Fluid Low Reynolds Numbers 
Passing through Infinite Row Equal Parallel Circular Cylinders,” Tamada, 
Quart. Jour. Mech. Appl. Math., 10, 1957, pp. 425-432. 

«Physics Water Movement Porous Solids,” Philip, Highway Re- 
search Board, Washington, Special Report, 40, 1958, pp. 147-163. 


the theoretical derivations Darcy and Forchheimer formulas,” Irmay, 
TAGU, 39, 1958, pp. 702-707. 
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GENERALIZED DISTRIBUTION NETWORK HEAD LOSS 


ASCE.--While the method described has been used and 
found applicable gas distribution networks, felt more limit- 
its application for several reasons. The chief reasons being the differ 
ences the character the system load the gas industry compared 
the water industry, and the differences gas network design reasonof deal- 
ing with compressible fluid. 

The widespreaduse gas for space heating produces large swing 
not unusual find the maximum day the winter five six times 
the average summer base load. The load pattern gas system seldom 
smooth swing from maximum day summer base load. broken sever- 
places when blocks interruptible customers are permitted the 
line, thus changing the flow pattern the entire network. Hence, compari- 
son pressure losses the method outlined would limited comparing 
loads occurring the same system flow pattern. 

Due the compressible nature the fluid handled, gas distribution sys- 
tems are numerous points district regulators that, turn, 
are supplied system operating higher pressure level. This, again, 
limits the the method outlined smaller segments the distri- 
bution network. 

Even with the limits mentioned, the writer has found the method valu- 
able and time-saving tool making rapid calculations the load-carrying 
capacity sections gas distribution system. particularly valuable 
the load dispatcher who quite often faced with making quick estimate 


pressure drops without having enough time for detailed computation ana- 
lyzer set-up. 


Engr., Laclede Gas St. Louis, Mo. 
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THE FOURTH ROOT n-f 


Corrections 


64, INTRODUCTION, third paragraph. line omit 
comma after “chart.” line replace “at” “to.” Fourth paragraph, line 
omit semicolon. 

Page 66, first paragraph. line insert “he” between “but” and “did.” 
line replace “is” “we.” line replace “sociated” “sociate.” 
line replace “curvilinear” “logarithmic.” 


Page 68. 
Page 69. 
Page 72. 
Page 73. 


Number.” 


Page 72. 


line replace the references “3,4” “8,10.” 

section iv, line insert “it” before “will.” 

line after Eq. insert “n” before “should.” 

subsection line insert “square root of” before “Reynold’s 


the equation preceding Eq. insert missing bracket. 


January 1960, Blench. 
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BOUNDARY LAYER STIMULATION RECTANGULAR 


Corrections 


the February Journal the Hydraulics Division, 
page 29, please delete the superscript that appears the second line above 
Fig. add subscript thedelta thatsame line. Also, the second 
line below Fig. please delete the superscript and Eq. please adda 
subscript tothe deltathat precedes the equals sign. page please inter 
change the illustrative material that presently identified with that 
which presently identified Fig. 13. That is, the captions remain their 
present relative position the page but the drawings are moved. 
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TRAP EFFICIENCY RESERVOIRS, DEBRIS BASINS, AND 
DEBRIS 


Discussion Rolland Kaser 


voirs presented the authors will certainly helpful engineers who must 
plan similar reservoir projects. Some major dissimilarities between 
voirs cited the authors have not been mentioned, however, and the facts 
should available those considering use the data. The writer has not 
examined the circumstances peculiar all the facilities listed the au- 
thor’s Table but those concerning the Colorado River facilities, with which 
familiar, will serve illustrate the point. 

The Imperial Dam reservoir, Yuma, Ariz., listed “Normal Pond- 
Reservoir,” Table with drainage area 184,000 miles, and data 
are given for two 5-yr periods operation (1938-1942 and 1943-1947). 
illustrated the map Fig. Imperial Dam situated downstream from 
Hoover Dam (drainage area 167,800 miles; closure form Lake Mead 
February, 1935); Davis Dam (controlling 1,520 miles downstream from 
Hoover Dam;closure toform Lake Mohave 1950);and Parker Dam (control- 
ling 9,500 Davis Dam;closure toform Havasu Lake 
October, 1938). Imperial Dam was closed April, 1938. The drainage 
area which contributes sediment the reservoir formed Imperial Dam 
limited the 5,800 miles downstream from Parker Dam. 

Inflow Imperial Dam Reservoir, however, almost entirely up-river 
runoff released through Parker Dam. water, which has been 
released through the gates turbines Parker Dam, picks new sediment 
load from the bed and banks the Colorado River downstream from Headgate 
Rock was this coarse river-bed material, from nearly 100 miles 
degrading channel, that filled Imperial Dam Reservoir. 

Throughout the 10-yr period cited (1938-1947), the lower Colorado River 

was adjusting new regime brought about construction and operation 
new dams, reservoirs and diversion structures. Under these circumstances, 
the C/W and C/I ratios cannot expected bear the same relationships 
sediment production and trap efficiency they would for reservoir the 

same watershed which received unregulated inflows. noted that data for 
Imperial Dam Reservoir Fig. plot far the left the curve. would, 
course, only coincidence that data for reservoir would line 
with data for “normal” reservoirs. 

The authors have presented specific data for Imperial Dam Reservoir 

Fig. unfortunate that the peculiar conditions that apply were not empha- 

sized. The writer has called attention the source sediment the lower 


February, 1960, Charlie Moore, Walter Wood and Graham Renfro. 
Planning Engr., Study Commission-Texas, 980 and Bldg., Houston, 
Tex. 
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DISCUSSION 


Colorado River that explains the coarse nature the material transported into 
the reservoir. Such sediments would have different relationship between de- 
tention time and trap efficiency than would apply other particle sizes. An- 
other factor influencing the trap efficiency the practice weekly bi-weekly 
sluicing through the California sluiceway, Imperial Dam, remove sediment 
deposits, that might interfere with diversions the All-American Canal, and 
remove bars the sluiceway channel formed sludge pipe discharges 
from the All-American Canal Desilting Works. This practice, that was ef- 
fect during the period record cited, involved, sequence, (1) closure the 
canal headgates, (2) complete opening the sluicegates, (3) complete draw- 
down the reservoir, (4) closure the sluicegates until the normal reservoir 
pool level was restored, and (5) opening the canal headgates for normal di- 
versions. Extra water was from Parker Dam, synchronization with 
the sluicing operation, minimize the time required refill the reservoir 
that the canal shutdown period would short possible. The authors 
have not explained how the sediment data their Table were computed, but 
there nothing indicate that the above factors were taken into account. 
any event, quite unlikely that the relationship between trap efficiency and 
detention time, shown Fig. would apply any other reservoir. 

The foregoing comments and explanation also apply the All-American 
Canal Desilting Basins tosome extent. Certainly, the size the drainage area 
upstream from Imperial Dam (as listed column Table for record 
No. 38) had little influence the sediment entering the All-American Canal 
during the period record cited. 

The authors have presented, column Table figures for the average 
annual sediment production per square mile drainage areaand refer readers 
their reference paper4 Brune, for details. interesting note that 
rate 2,170 tons per mile givenfor Imperial Dam Reservoir, while only 
1,044 and 406 tons per mile are given for Lake Mead and the All-American 
Canal Desilting Basins, respectively. Obviously the rates given cannot all 
applied the drainage areas listed column and will most helpful 
anyone interested the data Table the authors will explain the deri- 
vation the figures and the limitations their use. 

appears from the data presented the authors that long the C/I 
ratio exceeds about 0.3, planners would justified assuming that all the 
sediment inflow would deposited and that trap efficiency matter requir- 
ing special study only ponds and relatively minor reservoirs. 
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SCOUR BRIDGE 
Discussion Blench 


ASCE.—The writer draws attention the collection data 
scour around bridge piers and other obstacles made Inglis, from 
1924 1942.10 The data severe scour around piers numbered and were 
for floods from 30,000 cfs 2,600,000 cfs. They appeared, with descriptions, 
Table 8-1 the reference and were plotted its Fig. 8-1 show the re- 
lation scour depth (measured from water surface bottom scour) the 
Lacey regime depth for the river the flood discharge; the line Fig. 8-1 
was excellent fit and indicated that, all other things being equal, scour depth 
varies the cube rootof discharge exactly the regime theory based con- 
trolled canals shows; the ratio regime depth averaged almost exactly 2.0. 
The writer has the Inglis data along with model data Inglis, 
the present author (some time ago), andof workers the University Alberta. 
The coordinates were chosen terms estimated discharge the 
sake uniformity and future applications; allowing for the discharge intensi- 
ties the rivers being averages over the whole river whereas for models they 
are discharge intensities attack, the river and model data were good ac- 
cord with each other and with the regime theory prediction variation the 
two thirds power discharge factors being equal. What 
when other factors vary has been 

Inglis has also, experiments model bridge piers, similar 
the recent Iowa ones, has plotted the results, and has given formula fit 
the graph. The writer has reduced this formula rather simpler form: 


where scour depth from water surface, regime depth the river, 
and the breadthof the pier projected right angles the direction at- 
tacking flow. 

interesting note that the mean depth flow between piers (as dis- 
tinct from the local scour round shown Leopold!2 follow the re- 
gime lawof variationas thetwo thirds power discharge intensity the bridge 
over the Powder River Arvada, Wyoming. 

February, 1960, Emmett Laursen. 

Prof. Civ. Engrg. and Cons. Engr. 


«The Behaviour and Control Rivers and Canals,” (with the aid models), Sir 
Inglis, Govt. India, Central Water Power Irrigation and Navigations Research 
Sta., Ponna, 1949. 


“Regime Behaviour Canals and Rivers,” Blench, Butterworths Scientific 
Publications, London and Toronto, 1957. 
Hydraulic Geometry Steam Channels and Some Physiographic Implica- 
tions,” U.S.G.S., Prof. Paper 252, 1953, Fig. 4.A. 
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The writer feels that Eqs. and should not have been derived from con- 
siderations total load since bed load that mainly effective 
mining depth and that notes that seems 
unrelated regime depth open objectionas standard for com- 
parison. Considering the evidence of*regime theory, plus the Inglis and Leo- 


pold observations rivers, the writer would replace Eqs. and their 
equivalents from 


and geometric coefficients various obstacles obtain scour 
discharge intensity and the appropriate bed factor. 
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